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‘Two-dimensional incompressible magneto- 
hydrodynamic flow across an elliptical solenoid 


By NELSON H. KEMP and HARRY E. PETSCHEK 
AVCO Research Laboratory, Everett, Massachusetts 


(Received 30 April 1958) 


SUMMARY 


An analysis has been made of the two-dimensional flow of an 
incompressible constant-conductivity fluid through an elliptically 
shaped solenoid containing a constant magnetic field directed 
normal to the flow plane. The effect of both Hall current and 
ion slip has been included in the generalized Ohm’s law used for 
the fluid. ‘The analysis is based on a perturbation procedure in 
two parameters, one being the magnetic Reynolds number &,, 
and the other the ratio S of magnetic force per unit area to dynamic 
pressure. Calculations have been carried to the first order in each 
parameter, and closed-form analytic expressions have been 
obtained for the force and moment on the solenoid, the current 
density, stream function, magnetic field and other pertinent 
physical quantities. 

It was found that, to the zeroth order, there is a force but no 
moment on the solenoid. ‘lo the first order in S, where the flow 
field is modified but the magnetic field is not, there is a moment 
and a force, the latter being anti-parallel to the zeroth order force. 
To the first order in R,,, where the magnetic field is modified but 
the flow field is not, there is a moment but no force. ‘Thus, to the 
first order the lift to drag ratio is the same as in the zeroth order. 
Graphs which illustrate some of the effects of angle of attack, 
fineness ratio of the ellipse, Hall current and ion slip, on the forces 


and moments are presented. 


SYMBOLS 
constantin Fourierexpansion  C, contour of elliptical solenoid, 
of (cyh;/C7),, see (4.7), c value of € on_ elliptical 
semi-major and semi-minor solenoid, 
axes of ellipse, D drag, 
magnetic field vector, d see (3.2), 
uniform magnetic field inside electric field vector, 
ellipse (dimensional), complete elliptic integral of 
constant in Fourier expansion the second kind, 
of (cx;/e7),, see (4.7), electronic charge, 
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F force, 

h fineness ratio b/a of ellipse 

I discontinuous factor which 
is unity inside ellipse, zero 
outside ellipse, 

j current density vector, 

K; constant in stream function 
bs’, see (4.28), 

K, K,, K,, K, constants defined 


below (4.32 b) and (5.9), 

k unit vector normal to flow 
plane, forming third vector 
of a right-handed system, 


k Boltzmann’s constant, 

L lift, 

AY, reference length, 

l coordinate along free-stream 


direction, 


M moment about centre of 
ellipse, 

m particle mass, 

n coordinate perpendicular to 


free stream direction, 
N wake half-width, see (4.19), 


N; number density of j species, 

Pr perimeter of ellipse, 

p fluid pressure, 

q velocity vector, 

q° free-stream speed (dimen- 
sional), 

O;, elastic collision cross-section 
for collision between species 
jand k, 

R, magnetic Reynolds number 
uLaq’, 

R moment arm, 

r position vector from centre 
of ellipse, 

S dimensionless parameter 
oL (B°)*/pq°, 

s arc length along a contour, 
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z temperature, 

U constant vector, defined 
below (5.9), 

is vector, see (4.32), 

v° constant vector, defined 


below (5.9), 

W., W, constants in vorticity jump, 
see (4.20c), 

x,y coordinate along and normal 
to the major axis of the 


ellipse, 

z complex position vector, 
Z=Xx+1y, 

4. angle of attack, measured 


between major axes of ellipse 
and free stream, see figure 1, 


z, 8, y, 8 constants in the solution 
of Laplace’s equation (with 
various sub- and super- 
scripts), 

A. see (4.7 a), 

€ eccentricity of ellipse, 

v (1—h’), 

Ej friction coefhicient between 
species j and k, see (A3), 

‘é €&+ in, 

1 angular elliptic coordinate, 

6 polar angle in (x, y)-plane, 

K Hall coefficient, see (A2b), 

n\ ion slip coefficient, see (A 2c), 

be permeability, 

E radial elliptic coordinate, 

p fluid density (dimensional), 

o electrical conductivity of fluid 
(dimensional), 

d electrostatic potential, 

E= —V4, 

yb stream function, 

= —kx Vu, 


Q vorticity vector, Q= V xq. 


Subscripts 
A applied, E electrons, 
e on the elliptical contour, t inside ellipse, 
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Subscripts—continued. 


I ions, w in wake, 
N neutral atoms or molecules, x,y components of vector in x-, 
ln components of vector in /-, y-directions, 
n-directions, 1 unit vector, 
0 outside ellipse, v outer normal to contour. 


Superscripts 
0 zero order, a first order in R,,, 
first order in S, . dimensional quantity. 


~ 


1. INTRODUCTION 

The possibility that magnetohydrodynamic forces may be_ used 
advantageously in high speed flight has recently been suggested 
(Kantrowitz 1953; Patrick 1956; Rosa 1956). At the gas temperatures 
encountered at speeds corresponding to a re-entering satellite, air becomes 
a reasonably good conductor of electricity (Lamb & Lin 1957). If this 
air is allowed to flow through a magnetic field, drag, lift, or control forces 
may be obtained. It can be shown, on the basis of very rough order-of- 
magnitude arguments, that it may be possible to obtain larger forces with 
less heat transfer by using magnetohydrodynamic forces than would be 
obtained by the conventional method of using gas pressure acting on solid 
surfaces. 

In order to make a more precise evaluation of the possibilities of 
magnetohydrodynamics in flight, it is necessary to develop methods for 
treating magnetohydrodynamic flow problems. The work reported here 
is intended as a step in this direction. Its aim was to analyse a physically 
realistic magnetohydrodynamic flow problem which maintains some of the 
features of hypersonic flight. However, some of the simplifications that 
have been introduced to obtain analytical solutions are such that the direct 
application of the results to flight is not possible. In particular, we will 
assume incompressible flow and uniform electrical conductivity in the 
entire flow field. While this is a case which is physically realizable in either 
subsonic flows of hot gas or liquid metal flows, it is not a good approximation 
to hypersonic flight. Nevertheless, it is hoped that the understanding 
gained will lead towards the solution of realistic flight problems. 

For magnetohydrodynamic calculations associated with astrophysical 
conditions or with the conditions obtained in highly pinched discharges, 
the magnetic Reynolds number (see (2.9b)) is usually assumed to be very 
large. This is frequently stated more crudely as the assumption of infinite 
conductivity. ‘This is justified in astrophysics because of the large lengths 
involved and in the pinched discharges because of the high electrical 
conductivity which is achieved at very high temperatures (~ 10°°K). 

2N2 
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For large magnetic Reynolds numbers the flow and the magnetic field are 
very closely coupled and the field lines may be considered almost rigidly 
attached to the fluid particles. ‘The value of the magnetic Reynolds number 
that can be expected in flight is, however, of the order of unity or less. 
The magnetic field is then only slightly distorted by the flow. The general 
analysis in this paper is oriented towards small magnetic Reynolds numbers 
and the final calculation is a perturbation calculation based, in part, on an 
expansion in powers of the magnetic Reynolds number. 

In the presence of a magnetic field, the relation between the current 
and applied electric field—Ohm’s law—in a gas is considerably more complex 
than it is in a solid or liquid conductor. At low densities and high magnetic 
field strengths the electrons perform many revolutions in their helical orbits 
in the magnetic field between collisions. Under these conditions the current 
may be reduced and is not parallel to the electric field, as seen in coordinates 
moving with the gas. At still lower densities, and in a partially ionized gas, 
an appreciable ditference between the mean velocity of the ions and the gas 
velocity arises. ‘This results in a reduction of the effective conductivity 
of the gas. Both of these effects will probably be important under the 
conditions of interest for flight applications. The generalized form of 
Ohm ’s law utilized in this paper is valid for general values of the parameters 
describing these two effects and is restricted only by the assumption that the 
degree of ionization is small. 

A type of configuration which might be of interest in hypersonic flight 
is an extended magnetic field such as that produced by a large ring current 
surrounding a relatively small body. Such a configuration would produce 
a drag of the order of magnitude of the dynamic pressure acting on an area 
comparable with the area enclosed by the ring. The area of solid surface 
exposed to high heat transfer rates is, however, comparatively small. In 
such a case the flow interacts principally with the magnetic field, and the 
the presence of the body and the field coils may be 


disturbance due te 
neglected to a first approximation. 

In the present paper this type of configuration is idealized by considering 
a two-dimensional flow perpendicular to the axis of an infinitely long solenoid 
(figure 1). ‘The solenoid is of elliptical cross-section with its major axis 
at an angle x to the flow direction. It has a uniform magnetic field normal 
to the flow plane inside the ellipse, and no field outside. We assume that 
the solenoid itself is completely transparent to the flow; that is, the current- 
carrying Wires are thin enough not to affect the flow. The only interaction 
with the fluid is that caused by the magnetic field inside the solenoid. To 
further simplify the problem, we take the fluid to be incompressible, with 
constant conductivity. 

The method of solution of the problem is based on a perturbation in 


two parameters. One, as mentioned above, is the magnetic Reynolds 


number which describes the perturbation of the magnetic field about the 
undisturbed field of the solenoid. The other is a parameter which describes 
the perturbation of the flow about the uniform flow which would exist 
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in the absence of a magnetic field. ‘This second perturbation procedure 
is similar to thin-aerofoil theory, in that it is a small disturbance of a uniform 
flow. The magnetohydrodynamic case is, however, not restricted to slender 
cross-sections of the solenoid since for weak magnetic fields the effect on 
the flow field will be small regardless of the shape. In connection with the 
remark about aerofoil theory it might be noted that, from an aeronautical 
point of view, the solenoid behaves as a ‘ magnetic aerofoil’, disturbing the 
flow and producing a force and moment on itself, as does a solid wing. 








Figure 1. Notation and coordinate systems. 


The first-order changes in both the flow and the magnetic field have 
been computed. The resulting forces and moments on the solenoid are also 
determined by considering the corresponding forces and moments due to 
the magnetic body forces on the fluid. 


2. BASIC EQUATIONS 

The basic equations for the steady motion (with velocity q*) of an 
incompressible fluid of density p and constant conductivity o in the presence 
of a magnetic field B* and electric field E* are a combination of the Euler, 
continuity, and Maxwell equations. In rationalized mks. units they are 
plqt.V)qt+Vtpt—jtxBr=0, V.qt=0, — (2.1a,b) 

V* x EY s @, Ve.» = 6, V*x BY = py, (2.2a, b,c) 

where y» is the magnetic permeability and j* the current density. 
Equations (2.1a,b) are the usual incompressible flow equations except 








ui 


58 Nelson H. Kemp and Harry E. Petschek 


for the presence of the term j* x B* which represents the body force 
exerted by the electromagnetic field on the fluid. (The remaining Maxwell 
equation which relates the divergence of the electric field to the net charge 
density may in all cases of magnetohydrodynamic interest be replaced by 
the quasi-neutrality condition which states that the electron and ion 
densities are equal at all points in the fluid*. ‘This condition enters into 
the derivation of the generalized Ohm’s law for the gas and is not required 
directly in the macroscopic equation.) For a partially ionized gas with a 
low degree of ionization, the current density is related to the electromagnetic 
and velocity fields by a generalized Ohm’s law of the form (see Appendix A) 
i = o[ E* + q* x B* — x*j* x B* + A*(j" x B*) x B*]. (2.3) 

We will work with a dimensionless form of the basic equations obtained 

by introducing the following dimensionless variables: 

ae ee ae p- &* _F i me 
q’ ? p(q’)" ‘ B® ; gq’ B° ‘ aq’ B° : 


Here g°® and B° are the free stream velocity and the constant field inside 
the solenoid, respectively. ‘The coordinates are also made dimensionless 
with a characteristic length “&. ‘Then (2.1), (2.2) and (2.3) become 


(q. V)q = 1Vg?—-qx Q= —Vp+SjxB, V.g=0, (2.5a,b) 
VxE=0, V.B=0, VxB=R,j, (2.6a, b, c) 
j= E+q~xB-—«jx B+A(jx B)xB (2.7) 


(2.4, b,c, d, e) 


where 2 = V xq and the dimensionless parameters « and A are defined by 
Kk = x*oB?, A= A*o( 30 )2, (2.8 a, b) 

The two dimensionless numbers S and R,, which appear are defined by 
af (BP 
pq” ? 
R,,, is the magnetic Reynolds number and S is a parameter which represents 
the ratio of electromagnetic body force per unit area to fluid dynamic 
pressure. Equation (2.6c) shows that the magnitude of R,, determines 


S = R,, = pod TZ. (2.9a, b) 


the magnitude of the changes in the magnetic field due to the currents in 
the fluid. 

The parameters that appear in the non-dimensional form of the equation 
depend, of course, on the method of non-dimensionalization. The non- 
dimensionalization of most of the parameters is fairly straightforward. 
There is, however, some freedom as far as the current is concerned. For 
small R,, the non-dimensional current defined by (2.4e) will be of the 
order of unity. However, for large R,, the flow and magnetic field adjust 
themselves so that the current is much less than og®B°. In this case the 
magnetic field can change by its own order of magnitude in the characteristic 
length. The current should, therefore, be non-dimensionalized with respect 


* For a discussion of this point see, for example, Cowling (1953). 
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to the magnetic field gradient B°/u/. If this is done the parameters entering 
into the non-dimensional equations are R,, and S/R,,. It would seem, 
therefore, that the appropriate parameter for measuring the electro- 
magnetic forces is S/R,, for large R,,, and S for small R,,. The parameter 
S/R,, is frequently found in the literature because of the preponderance of 
interest in the large-R,, case. In the present case of small R,,,, we will use S. 
Equation (2.7) can be solved for j by expanding the triple vector product, 
forming j x B, and then substituting the resulting expression back into (2.7). 
This leads to 
_ (1+AB?)(E+q x B)—«(E+q x B) x B+ [k?+A(1+AB?)][(E+q x B). B]B 
i (1 +AB?)? + «2B? 





(2.10) 


as the generalized Ohm’s Jaw with Hall current and ion slip. 

In the present two-dimensional case, where the applied magnetic field 
is normal to the plane of the flow we may also take the induced field to be 
normal, so that (2.6b) is automatically satisfied. ‘Then since q is in the 
plane, the last term in j vanishes. Furthermore, in view of (2.5b) and 
(2.6 a) we may introduce a stream function and an electrostatic potential 
by the relations 


q = V 4 kus = Vis x k, E ==_— V4, (2.11 a, b) 


where k is the unit vector normal to the flow plane. If we also introduce B, 
the magnitude of the magnetic field, by 





B= kB, (2.12) 
then the equation for j becomes 
(1+ AB*)(Vd+ BV) —KB(Vd+ BV) xk 2.13) 
j= (1+AB2P + eB - 
In view of (2.12) we also find from (2.6c) that 
V x Bk = VBxk = R,, j, (2.14) 


and of course, in general, from (2.6c)* 
V.j= 0. (2.15) 
The nature of the body force in the present two-dimensional case is 
of interest. Let us break up the magnetic field into two parts. Let B, 
be the applied field produced by applied currents, as for example the 
constant field inside the ellipse in the present problem. Let B, be the 
induced field, which may be jooked on as produced by the induced currents 
according to (2.14). ‘The body force on the fluid involves only the induced 
currents, since they are the ones flowing in the gas. Therefore, the body 
* This equation could, of course, also be derived directly from charge conserva- 
tion. Sources and sinks for gas current could exist in the flow if electrodes were 
introduced. If the current in both the gas and the electrodes is considered it must be 
divergence free. However, the gas current itself would appear to have sources or 
sinks at the electrode surfaces. In this paper we will not consider electrodes, although 
the work here could easily be extended to cover this case. 
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force may be written with the help of (2.12) and (2.14) as 
(Vx B,) x (B,+B,)/R,, = — VB?/2R,,—B, VB,/Rn- 


This shows that in the two-dimensional case the body force due to the 
induced field is the gradient of a scalar, and wherever the applied field 
is constant, the whole body force is such a gradient. 

The Euler equation (2.5a) becomes 


(q. V)q = 3V@?—qx Q = —V(p+SB?/2R,,)—SB,VB,/R,, (2.16a) 


Thus, in the regions where B, is a constant, we have the usual hydrodynamic 
case of body forces derivable from a potential SB,(B,+}B,)/R,,, 
Bernoulli constant along each streamline. Further by taking the curl of 
(2.16a), with due regard for the two-dimensional nature of the flow, we 
find that 


with a 


(q.V)&=0, where B, = constant. (2.16b) 


Here we have recognized that Q has only one component and put 
Q=Vxq= Ok, Q= —Vu. (2.17 a, b) 


The last relation, coming from (2.1la), is the usual one between the 
vorticity and stream function in plane incompressible flow. Equation (2.16b) 
expresses the familiar fact that the vorticity is constant along streamlines. 
In the present case, it holds everywhere except at the ellipse C,, across which 
B, is discontinuous, so the last term on the right of (2.16a) cannot be 
expressed as a gradient. ‘Therefore, vorticity is generated only on C,, 
where B, jumps, and nowhere else in the flow. 

To determine how much vorticity is generated on C., consider for the 
moment an applied field B, which varies rapidly but continuously near C,,. 
Rewrite the last term on the right of (2.16a) in its original form Sjx B,; 
where j is understood to be the induced current. ‘Then take the curl and 
so obtain 

V x [(qx 2) +S(jx B,)] = 0. 
Since the tangential component of a curl-free vector must be continuous 
across a surface, the jump in vorticity across the surface is defined by 
requiring that the tangential component of 
qx 2+ S(jx B,) (2.18) 
be continuous. 

The other boundary conditions on C, may be obtained by similar 
considerations. Since the electric field is also curl-free its tangential 
component must also be continuous. ‘The induced current density in the 
gas does not become infinite anywhere. The induced magnetic field is 
therefore continuous everywhere. In other words the jump in total magnetic 
field at C, is just equal to the jump in the applied field. The finite current 
also implies the absence of infinite forces. The velocity and pressure are, 
therefore, also continuous. In the absence of electrodes the gas current 
is everywhere divergence-free and, therefore. the normal component of 
current is continuous on C,,. 
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In summary, then, the problem of determining the hydrodynamic and 
electromagnetic fields consists of solving 


V.j=0, VBxk=R,j, (q. V)Q2 = 0, 


both inside and outside C,. On C,, where the applied field has a discontinuity, 
but there is no body, we require: 


q, j,, o¢/és, B, are continuous on C,, (2.19 a) 


where d/ds denotes differentiation along the contour and v the outer normal 
to C,. Far from C,, the induced fields must vanish, except possibly the 
velocity in the ‘wake’ of the contour. Let1, be a unit vector in the direction 
of the free stream flow. Then 


(q-1,, j, E, B,)— 0 at infinity, except in the wake. (2.19b) 


The remaining condition defines the vorticity jump on C, by continuity 
of the tangential component of (2.18), 
Q) -Q; = S(j,/¢,)B4; on C., (2.20) 


where the subscript 7 means inside C,. ‘To complete the formulation, one 
must use the relations (2.11), (2.13), and (2.17) between Q, y, ¢ and j. 

The set of equations (2.14)-(2.16) are non-linear in % and B. In 
order to obtain approximate analytical solutions a perturbation technique 
will be used, valid for small values of the dimensionless parameters S and R.,. 

The quantities of main interest are the electromagnetic force and 
moment on the solenoid, which are equal and opposite to those on the 
fluid. Let F* and M* be the force and moment on the solenoid, which 
arise solely from the body force j x B on the fluid. Then 





= = =— |] x B) dA Pp 
og’ L 2( Bo)? p(q°)?£ i | | (J )e ? ( ) 
M* 


M = [APB ~ ia rx (jx B)ddA = } (r.j)dA, (2.22) 
where the area integrals are taken over the whole plane and r is the 
dimensionless position vector of a point so that M is taken about the 
origin of coordinates. By analogy with wing-section theory, one might 
be interested in the drag amd lift of the solenoid, that is, the forces parallel 
and perpendicular to the stream direction, respectively. ‘These may be 
defined as 
D* 


Fee > . Za 
D og’ £( B®)? F 1, ( 3) 


L* 
L =——__,, = F.. (kx1,) = Fn, 2.24 
og’ f 2(B°)? ( 1) 1 ( ) 
where 1, and n, are unit vectors along and normal to the stream direction, 
as shown in figure 1. 
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3. PERTURBATION EXPANSION 
All dependent variables are expanded in a double power series in S 
and R,,. Zero-order terms represent the impressed fluid motion and 
magnetic field. In the present system of units, the zero-order velocity is 
a unit vector 1, in the direction of the free stream. ‘The impressed field 
is normal to the flow plane and has magnitude unity inside the elliptical 
contour C, and zero outside it. ‘The symbol J will be used for this 
discontinuous quantity. ‘Then we may write 
q= 1,+ qS+ q’R,,+... 
= f+ Wb’ S+ LR, secs 
= O'S+ OfR,, +... 
I+ B’S+ BR, +...; 
pt yS+ jy Rt... 
b= 69+ ¢'S+ O'R, 4+... 
F= F°o+ F’S+ FR, +... 
M = M°+M’S+M°’R,, +... . 
If these expansions are inserted into (2.13)-(2.16) and (2.18)-(2.22) 
and terms containing equal powers of S and R,, equated separately, we 


0 
Ww = 


kxq 
Q=Vxe 
B 

i 


obtain the following results. 


Zero order 


x. 3 d? 4 9 —". dh? be /, . 
a (1 +AI?)(Vd° + IV" 2 KI(Vd°+ IVy )x (3.1) 
¢ 


V.j° = (1 . AI?)(V7d° + IV*b°) d= 0. (3.1 b) 


Here d is a constant which has different values inside and outside C,. It is 





defined by 


d = (14+AI?)? + «7/?. (3.2) 

J’, &6°/es are continuous on C,, (3.1 c) 

Fo = — || (jk) dd =kx [| 5? dA, (3.1 d) 
M? = | { Zk(r.j°) da =k || (r.j)) a4. (3.1e) 


First order in S 
Here we first use (2.14) to obtain 
VB’ xk = 0. 
Since B’ is independent of the variable normal to the plane, this indicates 
that B’ is at most a constant. Since B’ is continuous on C, and vanishes 
at infinity, it must vanish everywhere. 
With B’ = 0 we then have 
= (1 +Al?)(Vd' + IVd') —K I(Vd' + IVE’) xk 
j= 


22 
= J.Ia 
d ( 
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Vij = —(1-+A02)(V24' +1V2y')/d = 0, (3.3b) 

(1, V)Q! = Vx (j?x kl) = I(k. V)j? — kV .j? = 0. (3.3) 

q’, j,» &d¢'/és are continuous on C,, (3.3 d) 

and 0-2, =p, (3.3 e) 


where the subscripts 0, 7 denote the regions outside and inside C,, respectively. 


F’ = — [| (xk) d4=kx [{ ja, (3.3) 
M’ = || Ik(r.j’) dA = k | (r.j') dA. (338) 


First order in R,, 
Here we first consider (2.16) and (2.20) which show that 


Q2—-Of=0 on C,. 

Thus {* is constant along zero-order stream lines with no jump at C,,. 
Since {7 is zero where the streamlines originate, it is zero everywhere, so 
Qt = V xq? = (0). 

Since we also have V.q@ = 0 from (2.5b), and since q is not discontinuous 

? 
q” is constant everywhere. ‘This constant must be zero to satisfy the 
condition (2.19b) at infinity. Therefore Vs* = 0. 
Now from (2.13) one finds j*: 


ak 2 dt —a q a x o's 0 Ra 
ie (1+Al*)Vd oat k+V°B (3.4) 
¢€ 





Here V° is a vector which is made up entirely of zero-order quantities, and 
which will be defined later. It is different inside and outside the ellipse. 


_ (L+A2)V26+ VO. VBt+ BV .VO 





V.j¢ = ; 0, (3.4b) 
VE xk = ff. (3.4) 

j’, 0b"/es are continuous on C,. (3.4d) 

Fe kx {{ jed4— [| (j°«B*) aA, (3.4e) 

M* =k {{ (r.jr) d4— [| rx (j?x BY) dA. (3.4f) 


Equations (3.1), (3.3) and (3.4) describe the present problem completely 
to the first order in the parameters S and R,,, and their solution will be 
presented in the next section. Readers who wish to do so may pass on to 
the sections entitled Results and Discussion. 


2, 


4. SOLUTION OF THE PERTURBATION EQUATIONS 
The solution of the perturbation equations is fairly straightforward 
but the algebraic details are quite lengthy. Here we shall omit most of the 
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algebra. Readers who are interested in the details may refer to Kemp & 
Petschek (19538). 

‘The basic equations to be solved are (3.1b), (3.3b) and (3.4b) for the 
electrostatic potentials 6°, ¢’ and ¢’. ‘They are actually Poisson equations 
because in each case the terms other than the Laplacian of the potential 
are known. Either they are given, as in the case of ys", or they are determined 
from zero-order quantities; ys’ is found from j’ through Q’ by way of (3.3 e) 
and B* from j’ according to (3.4c). Once the various ¢ are found the 
currents can be obtained from (3.1a), (3.34) and (3.4a) and then the forces 
and moments from (3.1d, e), (3.3f, g) and (3.4e, f). 

[In actual practice it turns out to be easier to proceed slightly differently. 
The zero order and first order in S cases satisfy identical differential 
equations and boundary conditions, except that ¥° is given while ¥’ must 
be determined, as pointed out above, from j® through 2’. Therefore, the 
solution of (3.1 b) and (3.3 b) as Laplace equations in ¢ + Is, can be carried 
out simultaneously. ‘Then the appropriate values of % can be inserted to 
Lb’ is found from the zero-order 
solution. On the other hand, to find ¢%, equation (3.4b) is treated as a 
Poisson equation with V° and Bb“ determined from zero-order quantities. 


complete the solution in each case, after y 


The natural coordinate system to use is an elliptic one so we introduce 
elliptic coordinates €, 7 by the transformation 


z=x+iy=re'’ = eacosh(€+i) = eacoshl. (4.1) 


Here x, y are Cartesian coordinates parallel respectively to the major axis 2a 
and the minor axis 24 of the ellipse C, (€ = c), and « is the eccentricity 
given by 

= ]—f’, h = dja. (4.2) 


Each value of & defines a member of a confocal family of ellipses with foci 
at x = +ae. On each ellipse 7 is an angle running from 0 on the positive 
x-axis to 27. 

Since we deal with Laplace and Poisson equations we must have available 
solutions of Laplace’s equation in elliptic coordinates. ‘The regions inside 
and outside the elliptical contour C,, denoted respectively by subscripts 7 
and 0, must be treated separately. Furthermore, since the gradients of ¢ 
and ys are physical quantities, they must be finite everywhere and vanish 
at infinity. Suitable solutions of Laplace’s equation, which can be found 
by separation of variable in the elliptic coordinates, are 


f,=%,+ > %,, cosh mé cos my + B,, sinh mé sin mn, (4.3 a) 
fo = YotVeEtygntVingnt DL Com (Ym COs mn + 8, sin mn) (+.3b) 
1 


where z, 8, y, 5 are constants. 

With these solutions available it is only a matter of algebra to find ¢ 
in terms of ¢ for the zero order and first order in S cases, while for 6“ some 
simple particular solutions must also be constructed. 








en 

















Magneto-hydrodynamic flow across an elliptical solenoid 565 


Expressions for the zero order and first order in S potentials, forces and moments 


As pointed out above, the zero order and first order in S equations for ¢ 
are the same and may be handled together. According to (3.1 b) and (3.3 b) 
they both satisfy 

V7 +1V*s = 0, (4.4) 
and the boundary conditions that, on C,, j, and cd/cs are continuous. 
In both cases, j is also given by the same expression, (3.la) and (3.3 a). 

A solution to (4.4) is obtained by using (4.3a) for ¢;+%,; and (4.3 b) 
for ¢,. ‘The constants are found from the boundary conditions, and can 
5,6n ON C,. If we write 


j 


be written in terms of the Fourier coefficients of ey 


(cy,/en), = A,+ > A,, cos my + C,, sinmy, (4.5) 


m 


the constants are 


Ye=Ve, =%, 7, =A,, 

bm = (A, «—C,,[(1 +A)coth me + d,]}/mA,, cosh me, | 

Bm = (A,,[(1 +A)tanh me + d;] + C,, «}/mA,, sinh me, | (4.6) 
Yme™ = {A,,«+C,,[1+ (1 +A)tanh mc]}}/mA,,, ; 
6,6” = {—A,, [1+ (1+A)coth mc] +C,, «}/mA,,, | 
where 

A,, = 1+d;+ (1+A)(tanh mc + coth mc), (4.7 a) 
d,= (1+A)? +x. (4.7 b) 


The constants «, and y, are not significant and do not enter the equations. 

The current vector, both inside and outside C,, can be found from (3.1 a) 
or (3.3a) by using (4.5) and (4.3), but is too lengthy to write here. When 
written in Cartesian vector components it has the form of a series, each term 
of which is the product of an exponential function of € and a trigonometric 
function of 7. 

The forces and moments are obtained from (3.1d, e) or (3.3f, zg) and 
require the integration of j, dd and r.j,dA over the inside of the ellipse. 
Because of the form of j,, both these integrands also involve only terms 
which are the product of a function of € and a function of 4, so the integration 
is very simple. All but one term in each component vanishes because of 
the periodicity in 7, and the results are 

= —7a(x,hd,-y,y,)e~, (4.8 a) 
M = tkza*e*y,e-”. (4.8b) 

Now, in order to express F and M in terms of given parameters, it is 
only necessary to determine the Fourier coefficients A,, and C,, in the 
expansion of (Cci/7),, as indicated in (4.5). 


é 


Zero-order current, force and moment 
For zero order the stream function is, of course, given as that of a uniform 
stream at angle « to the major axis of the ellipse: 


Y° = ycosa—xsing = ae(sinhésiny cos%—cosh€cosysina). (4.9) 
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From (4.5) we see immediately that 

A? =CO=0, m#1; A! = bcosz, C?=asinz (4.10) 
When these are inserted into (4.6), the potential ¢° is completely determined 
(except for an additive constant which we shall ignore) by (4.3a) for 


d) +? and by (4.3b) for ¢. 
fp) = ye~x/at+d°ey/b, (4.11 a) 
dh? = e-*(y} cos 7 +4! sin 7). (4.11b) 


The current is likewise found from (3.1a) to be 


j; = X, ye C/b + y,o,é Cla =) x, +Ji Yu (4.12 a) 





si Xv — Y1 0!) (cos n — e~?*) + (x, 69+ y, y?)sin 27 
0 _ ( 1 yi : a ont 2 yi71) le (4.12 b) 
2a (sinh? € + sin* 7) 
This shows that the zero-order current inside the ellipse is a constant with 
Cartesian components 
Je, = ye*/b = [xcosa+sina(1+A+h-")]/A,, (4.13 a) 
Ji, = Sfe-*/a = [«sina—cosa(1+A+h)]/A), (4.13 b) 
where A, is found from (4.7). 
Since j? is a constant, the calculation of the zero-order force and moment 
directly from (3.1d, e) is simple: 
F° = (k x j)zab, M? = 0. (4.14) 
These, of course, agree with (4.8). The lift and drag, normal and parallel 
to the stream velocity respectively, are easily found by resolving j) along 
Cartesian coordinates m and / in those directions (figure 1). If we write 


jv = Lyi +My], (4.15) 
we find from (4.13) and (4.14) that 
L® = F°.n, = zaby') = zab[«k +sinacosa(h-!—h)]/A,, (4.16a) 


D® = F°.1, = —zabj’, = wab[1+A+hcos?~4+h-sin®«]/A,. (4.16 b) 


Discussion of these results is postponed to a later section. 


First order in S velocity and stream function 
To find the first order in S quantities one must know the stream function 
ys’. ‘This is determined by the vorticity Q’ according to (2.17b). The 
vorticity is in turn determined from zero order quantities by (3.3¢c, e). 
These equations are 
V4’ = -Q’, (1,. V)Q’ = 0, (4.17 a, b) 


Q’ — QO = j°/l,, on C,. (4.18) 


Ui U 


Equation (4.17 b) indicates that Q’ is constant along zero-order streamlines, 
i.e. in the /-direction, and is thus a function of only. Equation (4.18) 
specifies the jump in {’ at the ellipse caused by the discontinuity there in 
the applied field. The zero-order (free stream) streamlines which do not 
pass through the ellipse suffer no jump in vorticity, and since they come from 
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a region with no vorticity, they never have any. In other words, the vorticity 
is confined to the inside of the ellipse and to the wake. ‘The wake is bounded 
by the zero-order streamlines just tangent to the ellipse. In the (/,7) 
coordinate system these streamlines are defined by 

n= +N, N? = a®(h? cos* x + sin? x). (4.19) 


From (4.17) and (4.18), using (4.12) for j?, we find the vorticity inside 


the ellipse (/ and in the wake Q', to be 
Q. = W.— W,n[/(N?—n?), Q) = —2W,n//(N?2-n?), (4.20, b) 
where the constants are 
W, = —a*(h*j, cos % +); sin x) N=. (4.20c) 
W,, = @h(j,2sin x—j)) cosa)/N? = —a*hy?/N*. (4.20 d) 


With the vorticity known, the stream function follows by a straight- 
forward superposition of the stream functions of elementary vortices of 
strength 02’ dndl: 


1 -- sia = 
I,! = | y(n y—z 2 
om ee | _Q’(a)log's — 2| dndl. (4.21) 
om TJ tre 
Here z = x+iy and |s—2| is the distance between the vortex and the field 
point. 
One convenient way to find ¥%’ from this formula is to first obtain 
q, = —(Vy’), The integral for (Vy’), can be converted by integration 


by parts and use of the divergence theorem into an integral on C, with 
(Q'—Q)) in the integrand. Use of (4.20) then leads to 


In 


(V¥'), == | (hjgcosa +! 


, Sin H)log|s —z,| dip. 


The log term can be written in elliptical coordinates and expanded in a 
power series in exp[—(¢+¢,)] and exp[+(¢—Z,)], the plus sign being 
used inside the ellipse and the minus sign outside. ‘The integration on 7 
is then easily performed and the result, in the (/,m) coordinates, is 





, Ifj;—-UA—h)j2cos a] +n[j2 + 1 —h)j,2 sin « 
mn [hin— ( 2)J 2 COS =! Sus ( b)j;~ sin «| (4.22 a) 
G,, = —ahe-*[j9 cos(n — «) +7), sin(y — «)]/e. (4.22 b) 
To find #%’ we can now integrate the relation g, = —0oj’/el. First we 


work outside the ellipse. g,’, can be integrated with respect to / if we notice 
from (4.1) that 

de _ i ae(2t + e-2 )+ie mae 4.23 
ied bae(20+e-**), +in = ze“, (4.23) 
This integration determines y’ except for an undetermined function of n, 
which must be found from the vorticity relation (4.174), and from the 
boundary conditions on q,, = oy,/dn. This function turns out to be zero 
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outside the wake, but not inside it. The results for g,, and #’ are 

q,, = ahe~*[jsin(n — x) —j;) cos(y — x) ]/e+ 
0 outside the wake, (4.24) 
2 2W/ (N?2—2?) inside the wake, 


uf tha" y' (e-*§ cos 2n + 2£) yi (e-? sin 2y - 2n)] : 
-(~) outside the wake. 


= le forn>O\]. . 4.25 
-W,, N?| sintcost+t— {5 il inside the wake. ‘ ) 
a7 for n < U0 


Here ¢ ts an angle defined by 


n/N=sint, 0<t< 4m for n>0, 8a <t<27 forn<0. (4.26) 


An entirely similar procedure, beginning with g!., and keeping in mind 
5 I Y in £ 
the continuity of q’ on C,, leads to the velocity and stream function inside 


the ellipse : 


q; U oe . sin a] n [W.+ — cos @] — W./(N?-2n?), 
(4.27) 

wb. = 3(ng,,—1q;,)—4W, N*t+K, (4.28) 

where K;, is a constant chosen so that 4’ =u on C,. It differs for n S$ 0. 


As a check on the first order in S flow field, the velocities far from the 
ellipse can be found. Because Bernoulli’s equation holds everywhere but 
in the wake the pressure can also be found far from the ellipse except in 
the wake. An application of the integral momentum theorem to a large 
circle then yields the force and moment on the fluid due to the first order 
in S flow field. It is found that this calculation exactly checks the zero-order 
force and moment. ‘This shows the advantage of obtaining forces and 
moments from the body force, since to obtain them from the flow field 


requires a higher order calculation. 


First order in S force and moment 

Knowing 4s’, the force and moment are found by returning to (4.5), 
where they are expressed in terms of y,, 6, and y,. ‘To determine these 
constants according to (4.6), the coefficients 4,,, C,, in the Fourier 
expansion of ¢ys’/en on C, must be found for m = 1, 2. 

Since vs’ is continuous on C,, either yf’ or us’ may be used. When the 


expansion is written in the form of (4.5) we find that 


A’ = — — W,Nhacosa = {- at 2 VA 


3a 

, ee ae sige Nees 
C,= —- — W,Nasin« = & = WN) ¢ D. 
IT IT ‘ 


where (4.10) for the zero-order coefficients have been recalled. This 
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result shows that F’ is parallel to F’, and proportional to it: 


F’ ( = WeN)F 2 aA (4.29) 





32 4 (A? cos? x + sin? x) ” 
The second-order coefficients, which determine the moment, are 
ha*W., rae ) lene h{(i—h) .. 
= = “ (hcos* x—sin* x), C, = =| —— W, sin 2 oe W.. |. 
. L+h - ee l+h 
When these are inserted into (4.6) for y, e-*’, and that into (4.8 b), we find that 


a'a(1—h?) ( 2kh 








M’ k —_— :_—__ W , (h cos? a — sin 2x) 4 
SA, L+h 
{1 +(1+A)tanh 2c] E = W,.sin 22 oe W, | ». (4.30) 
2 th 


he quantities IV, and IW. are defined in (4.20c, d) and A, in (4.7a). 
For a circle # = 1, and we see that M’ = 0, as required by symmetry. 


First order in R., fie ld and pote ntial 
As has been pointed out before, to the first order in R,, there is a perturbed 
magnetic field, rather than a perturbed flow field as found to the first order 
in S. ‘This field is related to j’ by (3.4), from which, since Jj? is a constant, 
it is easy to see that inside C 
BY Je, +], 9. (4.31 a) 
Outside C,, j) = —V¢" according to (3.1a) and (3.2). Using (4.11b) 
for d” and (4.13) for y! and 5, we find that 
B! = e-*(y' sin yn — 8" cos 7) = ae’—*(Aye. sin y J) cosy). (4.31 b) 
To find the potential 6" we use (3.4b). ‘The vector V°, which is obtained 
from the expansion of (2.13), is given by 


V°d = (14 3A) Vib? + 2AIV 9 — 2k IV? x k— «V9 x k + 297 [2A(1 + Al?) + &?). 


(4.32) 

This is constant inside C, since all zero-order terms have been found to 
be constant there. Outside C, it is variable, but V.V° = 0, as can be seen 
from the definition (4.32). ‘Therefore, 6” satisfies the Poisson equation 
V24 Vo. VB(14A22). (4.33) 


‘This equation is solved both inside and outside, using (4.3) as the comple- 
mentary solutions. ‘The particular solutions are simple exponentials of 2€ 
and trigonometric functions of 27, obtained using (4.31) and (4.32) and the 
zero-order quantities already found. ‘The constants in the complementary 
solution are determined as before by matching j% and ¢e¢"/cy on C,. This 
results in a set of simultaneous equations whose inhomogeneous parts 
(which come from the right side of (4.33)), have only constant terms, and 


\\ 
terms second harmonic in 7. ‘The determinant of the coefhcients is d;A,,, 
which does not in general vanish and we conclude that only «‘, y’, y* and 
the coefficients of the second harmonic terms are non-vanishing. 

20 


F.M. 
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The resulting expressions for $” are 


¢! = Ka’e*(cosh 2é + cos 27) + 43 cosh 2€ cos2n + P38 sinh2€sin2y, (4.34) 
db! late2(K ,e-*§ — K,.cos 2n t K, sin 27) + e—*8(y" cos 2n + 62 sin 2n) ry. 
(4.34 b) 


The constants K, K. and A, come from the right-hand side of (4.33) and 
are detined below (5.9). ‘The constant Kp, arising from the same source, is 
K, =2(ae) 2(; +6") + (ae) (y'' cos x%—6" sin). 
y3, # and f§ are related to A", AJ and C) by (4.6); the latter two constants 

are defined below (5.9), and 


(1+A)abK,—ab(j?,U,—-ji,U,) , @e(1-hA)K, 


]” 
2d, 2(1+h) 





’ 


with U also defined below (5.9). 74 and 03 are given by 
vie? x" cosh 2c + La*e?(K + K,,), 


ote 34 sinh 2c — 1a*e?K.. 


First order in R,, current, force and moment 

The current is found by inserting (4+.3+) into (3.4a) and using (4.31) 
for B’. Inside C. the result is 
ji = K[(1+A)(x,x+y,¥)+«(—%,9+y,x)]/2d,+ V}(—-j) x +y).v)/d, 

+ (4/a*e?d;)\ [(1 + A)a3 — «B3](x,¥—y,y¥)+ 
\\ Qa ‘in ’ P va 
+ [(1 +A)B% + xaS](x,yv+y,x)}. (4.35) 
The current outside could be found also, but we will now show that the 
force and moment do not depend on it. 

The force and moment are given by (3.4e,f). The second integral 
in each formula can be simplified by using the relation j® = V x B’. The 
integrands become 

jx B LV (B2), rx (j° x Bt) = $V x 2(B*). 
Use of vector integral theorems for the divergence and curl then reduce 
the area integrals to line integrals which can easily be shown to vanish 
because 8B" is continuous across C, and vanishes as €-- 2 (see (4.31b)). 
The formulae for force and moment therefore reduce to 


F'=kx {| jtdd, M*=k~ || (r.j*) dd. 


They depend only on j’, not on the outside current j?. 


The integrations are very simple. Since j? is linear in x and y, the 


symmetry of the ellipse indicates that 
F’ = 0. (4.36) 


The only contributions to the moment come from x? and y? terms which 
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arise when r.j? is formed from (4.33). After some reduction, the result is 
M’ — atn(1—h?) (h(V oj +P yin) , Ade+C$[1+(1+A)2h/(1+h2))) 
7 + \2 d; A, . 








(4.37) 
V’, which is closely related to V°, is defined below (5.9). Again the moment 
vanishes for a circle (2 = 1), which is required by symmetry. 


5. RESULTS 
The results for force and moment to the zeroth and first order in both 
S and R,, are collected here for convenience of reference and discussion. 


Zero order 


Fe =] D +n, 2. (5.1) 
From (4.2), (4.7) and (4.16) the lift and drag are 
L° = nabj) = ma*h[« + 3 sin2a(h-! —h)]/A,, (5.2 a) 
D® = —zabj}, = za*h[1 +A+hcos? «+h sin? «}/A,, (5.2 b) 
where 
A, = (14+A) (1+h)° L 24 2. (5.3) 
h 


Also from (4.14) we recall that 


First order in S 
From (4.29) and (5.2b) the force is related to FY by 








F — 8 ahF°;!), _ _ 8 ah{1+A+hcos?a+h sin? x] FE 
3x / (h? cos? x + sin? x) 3a \/ (h? cos” x + sin? x) : 
8S F°D°® 
= — (5.5) 





32a \ (h? cos? x + sin? x)” 


From (4.30), (4.20c, d) and (4.13) the moment can be put in the form 








M’ a'x(1—hA)(1 +h)? (sin 2x i 2h? 
i = 9 7 (1 Ls h*) - 9 9 a zd 
k 8A,A,(1+h?) ( 2 h? cos? x + sin? x 

h 


| Tame g sin2%+«cos2x) x 


2h?(1 4 +h)? sin 2x = 
; (a .h) + 2h?(1+A) ) Ah(1 +h)? sin =|} (5.6) 


h? cos? «+ sin? h* cos? x + sin? x 








where 





First order in R,, 
From (4.36) we may recall that 
F2 = 0), (5.8) 
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‘The moment is too lengthy to write out in one expression. From (4.37) 





it 1S 
M wai(1—h?) (hV Yj? +0 jo Aw +C$[1+(14+A)2h/(1+A?)] 
k 4+ 2 (14+A)?+«? A, is 
(5.9) 
where 
ji nc Se ee ee 
2(1 +A)? + 2k? 4 ? 
ps a* h( jiy U, +49 at ae?(K + K..) 
2(1 +A)? + 2k? + : 
U= [(1+A)?+«?]V!— V?, 
V = [(1+Al?)? + 4?) (1 + SAD) Vib? + 2AIV 9 — 2k ITV? x kk - 
KV x k+25°7[2A(1 + Al?) + J}, 
K (ae) '(d'sin x—y')cos x), K, = (ae) 1(8) cosa + y'!'sin x), 
K Vek xj (1 +A), 
e° x«cosaz+sing[1l+A+h"] dfe"  «sinz—cosa[1+A+AhA] 
h A, a A, 
Vd" x, (y° e/a) + y,(6°e-*/b), Vu X, SIN % + Y, COs 4, 


x (ye “/b)+y,(d/e/a). 


j Xie tt YiJi, 


For the special case of « = A= 0, M” reduces to the following very 


simple expression 
! — 1J,2 
M’ 7 ath*(1 h) a (5.10) 


Si 


& lth 


6. DISCUSSION 


The problem that has been considered basically contains six parameters ; 
the interaction parameter S, the magnetic Reynolds number &,,, the Hall 
coefhicient «, the ion slip coefficient A, the ratio of minor to major axis of the 
ellipse A, and the angle of attack x. Since the analysis has been carried 
out as a perturbation in S and R&,,, the results for each order, as given in the 
previous section, are functions of the four parameters «, A, A and x. 
Discussion of these four-parameter functions is, necessarily, limited to 


consideration of certain limiting cases. 


Zero ordei 

‘To the zero order the gas current distribution and the resulting forces 
ire computed from the undisturbed flow and magnetic field. Since to 
this order the hydrodynamic properties of the fluid do not enter, the zero- 
order solutions are not restricted to incompressible fli WS. 

As may be seen from (4.12a) both the magnitude and the direction 
the ellipse are independent of position. While 


of the gas current inside 


this is true for all values of the parameters x, A, # and a, it cannot be 
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generalized to all closed contours, but seems to be a property peculiar to 
ellipses. It can easily be demonstrated, for example, that it would not 
hold for a rectangular cross-section. 

The fact that the zero-order moment vanishes for all ellipses follows 
directly from the fact that the current and, therefore, the body forces on 
the fluid are uniform throughout the ellipse. 

Let us consider first the special case of high densities where « and A 
are both zero and discuss the force on ellipses one of whose axes is parallel 
to the flow (z = 0 or }z). Equation (5.2a) shows that the lift is zero, and 
equation (5.2b) may be written in the dimensional form 

\ 
D®* = og? B’ra*b* (=n) (6.1) 


where H is the ratio of the lengths of the axes perpendicular and parallel 
to the stream; 1.e., 


H=h for «<=0, H=kh* for «= tn. 


The factor preceding the bracket is the drag the ellipse would have if there 
were no electric field to impede the current flow, i.e. if the current in the 
ellipse were oq’ B?. The factor in brackets is then the reduction in drag 
due to the fact that an electric field is required to close the current loops 
outside the ellipse. In other words, it is a measure of the electrical 
impedance seen by the e.m.f. which is generated inside the ellipse as a result 
of the motion of the gas through the magnetic field. For ellipses of fixed 
area the drag increases monotonically as the length of the axis perpendicular 
to the stream increases. 

The case of general angle of attack for « = A = 0 may be considered 
as the superposition of flows parallel and perpendicular to the major axis. 
Since the problem has been linearized, the total force is the vector sum of 
the forces resulting from the component flows. Because the drag per unit 
velocity is larger for the flow perpendicular to the major axis than it is for 
the flow parallel to it, the total force will not be a pure drag force but will 
also have a lift component. 

In figures 2 and 3 the drag and lift for ellipses of different fineness 
ratios / have been plotted as a function of angle of attack from (5.2). The 
perimeter P has been used as the representative length Y in this plot, 
so that the forces are compared for ellipses of the same perimeter. This 
choice was motivated by the fact that for a given magnetic field and a given 
mass of copper in the solenoid the joule dissipation required in the copper 
is a function only of the perimeter. On this basis, in addition to the factor 
due to the change in impedance with shape discussed above, the drag 
decreases as the ellipse becomes thinner because the area decreases. 

The lift is zero for a circle because of symmetry, as mentioned above, 
and also zero for a very fine ellipse because the area goes to zero. It has 
a maximum fora fineness ratio of about 0:3. The ratio of lift to drag becomes 
infinite for infinitely thin ellipses at small angles of attack. At this point, 
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of course, both the lift and the drag themselves approach zero. For a 
reasonable fineness ratio the lift-drag ratio is not large. For example, for 
a fineness ratio greater than 0-2 the lift is less than the drag at all angles of 
attack. 
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Figure 2. Non-dimensional zero-order drag D® = D®*/aq°P?(B°)? (based on peri- 
meter, Y P) vs angle of attack « for varving fineness ratio h; plotted from 
equation (5.2 b) with Hall and ion slip coefficients « and A equal to zero. E is 
the complete elliptic integral of the second kind. 


In order to describe the effect of the Hall and ion slip coefficients « and A 
let us consider first a circular solenoid, h = 1. The Hall coefficient tends 
to produce a current in the stream direction. The cross product of this 
current with the magnetic field produces a lift even for a circle. Since the 
sign of the Hall current is independent of the sign ot the magnetic field, 
the sign of the lift will depend on the sign of the magnetic field. Lift and 
drag have been plotted as a function of « in figure 4 from (5.2). In drawing 
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this graph we have assumed that 
A = «7/500, (6.2) 
which is a limiting case of (A5). 
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Figure 3. Non-dimensional zero-order lift L° = L°*/aq°P?B° (based on perimeter, 
gY P) vs angle of attack « for varying fineness ratio h; plotted from equation 
(5.2 a) with Hall and ion slip coefficients k and A equal to zero. F is the com- 


plete elliptic integral of the second kind. 


As « increases from zero the lift initially increases and then decreases. 
However, both the drag and the total force decrease monotonically. 
Lift-drag ratios up to 7-9 can be obtained, but these large ratios occur when 
the total force has been reduced by almost a factor of ten compared with 
the case « = 0. 
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As can be seen by looking at the general formula for zero-order drag 
(5.2b) the reduction in drag associated with x depends on the geometry. 
It becomes less significant for flat ellipses. ‘Vhis reduction can probably 
also be made less significant by the use of electrodes whose potentials are 
adjusted to inhibit the flow of the Hall current. Non-uniform conductivity 
in the flow field, which would be encountered in an actual flight case, might 
also tend to produce a similar effect. 





Figure 4. Non-dimensional zero-order lift L° and D® zs Hall coefhicient « for a 


circle, / 1; plotted from (5.2.a, b) with ion slip coefficient A K7/500. 


It is instructive to examine two limiting values of « and A for general 
g 
ellipses. Equation (6.2) shows that « can be large compared with unity 
while A is still small. ‘This case will be approximated by «—- «x, A = 0. 
In the second case examined, both «-- «2, A> x. 
In the first case, the lift is much larger than the drag. ‘The dimensional 


form of the lift may be obtained from (5.2a) as 


a. 


aq® BY xa*bh* ‘ 
! (NV, eg") B°za*b*. (6.3) 


K 


In order to obtain (6.3 


equations (A2b) and (2.8a) have been used to 
define x. Since the current inside the ellipse was shown to be constant, 


the force is the product of this current, the (constant) magnetic field, and 


the area. ‘Thus, the term in parentheses in (6.3) must be the current, 


In order to produce this current the difference between the ion and electron 
velocities must be g’. For A = 0, (A 1b) shows that the ion velocity is equal 


to the gas velocity g°. The mean electron velocity must therefore be zero 


inside the ellipse. This is phvsically reasonable since, for « large, the 
I pry ; g 


electrons make tight spirals around the magnetic field lines. When the 
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electron suffers a collision the centre of the spiral is moved by something 
less than the radius of the spiral. Since in this region the radius of the 
spiral is small and collisions are relatively infrequent, the electrons are 
practically stationary in the magnetic field. We may conclude that in the 
limit of « large but A still small the current may be approximated by assuming 
that in the magnetic field the electrons are stationary and the ions move 
with the gas velocity. 

In the limit of both « and A large the drag becomes the dominant force. 
Using equations (A 2c) and (2.8b), we find that the expression for zero- 
order drag becomes 
og’ B’za*b* 


From the definition of ¢€,,, the term in parentheses is the frictional force 


D* 


(N, Ny €,y 9° )7a*b*. (6.4) 


due to collisions between ions and neutrals when their relative velocity is q’. 
In this case, then, the ions are held almost stationary in the magnetic field. 
It should be noted that in this limit the drag becomes independent of the 
magnetic field strength. 

When the ion slip term is significant, the ion velocity is appreciably 
less inside than outside the ellipse. ‘The assumption of uniform conductivity, 
which we have used throughout, implies that the ion concentration is uniform 
in space. In order to satisfy the continuity equation for the ions under 
these two conditions, a large portion of the ions entering the ellipse must 
recombine with electrons near the boundary. Also, where the streamlines 


leave the ellipse a corresponding number of ions must be produced. If 


the rates of ionization and recombination are sufficiently rapid, this require- 
ment can be satisfied. ‘The rapid rate of ionization would maintain thermal 
equilibrium everywhere. ‘The uniform temperature which is implied by 
the assumption of incompressible flow then assures uniform degree of 
ionization. However, since ion slip occurs at low densities, the rate of 
recombination may not be sufficiently rapid to give thermal equilibrium 
everywhere. In this case the degree of ionization inside the ellipse would 
be greater than outside. Equation (6.4) would then suggest that the drag 
would be higher than one might expect assuming a uniform degree of 
ionization. ‘The reduction in drag due to ion slip may, therefore, not be 
as large as the present calculations indicate, if the rates of ion recombination 


are not fast enough. 


First order in S 

To the first order in S we have caiculated first the modified flow field 
and then the resulting forces and moments. ‘The perturbation velocity 
perpendicular to the free stream velocity is given in (4.22a). The parallel 
component inside and outside the ellipse is given in (4.27) and (4.24). 
As was mentioned previously, for a circle the external flow outside the wake 
may be represented by a point source and vortex at the origin 

The gas pressure may be computed by substituting the first-order 
velocity and zero-order current and magnetic field into the momentum 
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eguation (2.5a). It is interesting to note that in general the pressure 
gradient inside the ellipse is a constant, and its component in the free 
stream direction is negative. For a circle this pressure gradient is equal 
to one-half the zero-order magnetic force. For an ellipse, it is more 
complicated and is not in general parallel to the magnetic force. 

The first-order forces (5.5) are of opposite sign to the zero-order forces. 
‘The fact that the forces must be reduced by the first-order terms is obvious 
since the flow velocity in the ellipse is reduced by the magnetic field. It 
is interesting to note, however, that the calculation shows that the first-order 
force is exactly anti-parallel to the zero-order force. 

‘The magnitude of the first-order force may be used as a crude indication 
of the range of validity of the perturbation procedure. ‘The ratio of first 
to zero-order forces is, from (5.5), 











'S 8 Ys 
SFr =~ Fpi (iPeosta sina) 12 = 
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‘The term in parentheses may be interpreted as a drag coefficient based 
on the major axis. If we take the square root as about unity, the first order 
., correction when the zero-order drag coefficient 
is unity. One would therefore expect that the terms which have been 
calculated are insufficient beyond this point. It should be remembered 
that when « and A are not zero the drag coefficient may be small for fairly 


in S produces about a 25° 


large values of S. The first-order terms become comparable with the zero- 
order terms only when the drag coefficient becomes comparable with unity, 
not when S becomes comparable with unity. 

‘The moment, which vanished to the zero order for ellipses, appears in 
the first-order approximation. In figure 5 the moment for « = A= 0, 
has been plotted from (5.6). The moment has been plotted in terms of the 
non-dimensional moment arm 1 F°*Y divided by the zero-order drag 
coethcient 2SD°. ‘The plot therefore compares the moment arm for different 
ellipses with the same drag. It will be noted that the moment is zero at 
angles of attack of both 0° (major axis parallel to free stream) and 90° (minor 
axis parallel to free stream). At 0° the neutral point is unstable, so that the 
ellipse would orient itself at the other neutral point with the major axis 
perpendicular to the free stream. ‘lhe sign of this moment may be explained 
crudely by noting that the flow through the downstream parts of the ellipse 
has already been decelerated somewhat in the magnetic field. ‘The flow 
velocity is then higher in the upstream part of the ellipse and therefore 
the forces are larger there. A small angle between the major axis and the 
free stream direction therefore produces a moment which tends to increase 
the angle of attack. 

For general « and A the moments become quite complex. One may to 


some extent regard the Hall effect as introducing another angle into the 
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problem, the angle between the current and the electric field in a coordinate 
system moving with the gas. The introduction of this angle tends to 
change the angle of attack for a given moment. In order to illustrate this 
effect the stable angle of attack for an ellipse with 4 = 4 has been plotted 
in figure 6 as a function of «x. This was calculated numerically from (5.6), 
and (6.2) was used to relate « and A. It will be seen that for « both very 
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Figure 5. Non-dimensional first-order in S moment arm R’ = M’/YF° divided by 
zero-order drag coefficient C® = D®/}p%q@@Y = 2SD° ws amen of attack « for 


varying fineness ratio h. Plotted from equations (5.6) and (5.2 b) with Hall 
and ion slip coefficients k and A equal to zero. 


small and very large the stable angle is 90°. For intermediate values the 
stable angle is reduced to almost 45°. In the limit «— co and A = 0 the 
stable angle would be 45°. For small values of « the stable angle depends 
on the ellipticity. However, for large « the curves for different ellipses 


approach one another asymptotically. 


First order in R,, 


The first-order calculations in R,, principally determine the changes in 
the magnetic field caused by the flow. The change in field is given by 
(4.31a) and (4.31b). Asa result of the uniform zero order current (4.12 a) 
inside the ellipse the field inside varies linearly with position, the field 
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strength increasing toward the downstream side. For « and A zero the 
magnitude of the field change across the ellipse is roughly R,, times the 
zero-order field. In this case the perturbed field is small compared with 
the applied field for R,, less than unity. For « and A not zero the current 
and therefore the field change is reduced and the ratio of perturbed to 
applied fields remains small to larger values of R,,,. 
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Figure 6. Angle of attack x, at which ellipse has zero moment to the first order in S, 
ind is stable with respect to that moment, zs Hall coefficient « for two fineness 
ratios / L, Plotted from equation (5.6) with ion slip coefficient 
K- S500 


‘The fact that the field variation is linear immediately implies that the 
a the downstream 


~ 


first-order force is zero (5.8). The increase in force o1 
side of the ellipse due to the increased magnetic field is exactly balanced 
to this order by the reduced force on the upstream side. 

‘he moments to this order for « = A = 0 have been plotted in figure 7, 
from (5.9). ‘The shape and magnitude are very similar to that for the first 
order terms in S. ‘The sign of this moment is, however, opposite to that 
of M’. ‘Therefore, for R, somewhat larger than S the stable angle of 
attack for the ellipse is 0°. ‘This may be seen qualitatively from the argument 


above that the force on the downstream end is larger than on the upstream 


end. 
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Figure 7. Non-dimensional first-order in R,, moment arm R4# = \/**/AYF® divided 
by zero-order drag coefficient C° D®/4 pq” F 2SD®° vs angle of attack « for 


varying fineness ratio /. Plotted from equations (5.10) and (5.2 b) with Hall 


and ion slip coefficients « and A equal to zero. 


Concluding remarks 

This analysis has been carried out as a step toward the solution of the 
more complicated situation which exists in hypersonic flight. Further 
work in this direction is presently under way at the AVCO Research 
Laboratory. Dr Frank Fishman is considering theoretically the corre- 
sponding perturbation problem for supersonic flow and Mr John Lothrop 
has begun an experimental investigation of the supersonic case utilizing a 


shock tube. 
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APPENDIX A. GENERALIZED OHM’S LAW 

The electrical conductivity of a partially ionized gas in the presence 
of a magnetic field has been calculated by Schluter (1950, 1951). A review 
of this work with some extensions and some algebraic corrections is given 
by Finkelnburg & Maecker (1956). Very briefly Schluter’s approach 
consists of considering the separate momentum equations of the component 
gases, ie. neutrals, ions, and electrons. ‘hese equations include, in 
addition to the acceleration, pressure gradient, and electromagnetic force 
terms, a frictional type of force arising from collisions between different 
components of the gas. This frictional force is proportional to the difference 
in mean velocity of the components. ‘These three momentum equations 
may be combined to give the overall gas momentum equation (2.1a) and 
in addition the mean velocities of the ions and electrons may be determined 
in terms of the local flow properties and their gradients. Instead of using 
the electron velocity explicitly it is more convenient to use the difference 
between electron and ion velocity, which is essentially the current. 

‘The result (equations (51.17) and (51.18) in Finkelnburg & Maecker 
(1956) for the case of arbitrary degree of ionization of the gas is quite 
complex. However, for the special case of very small degrees of ionization 
where the partial pressures of the ions and electrons may be neglected 
the equations for current and ion velocity may be reduced to the form 
jt = o(E* + qj x B* — «*j* x B*), (Ala) 
q; = q* +A*j* x B*, (A1b) 


where q7 is the mean ion velocity, and o, «* and A* are coefhcients which 
describe the gas conductivity. In terms of kinetic theory parameters 





N, e 
oo = __ — . (A 2a) 
Vy egy t+ Nrere/ Ny 
K* Ne, (A 2b) 
A* 1/N, Ny€zy; (A2c) 


where the subscripts \V, J, and FE refer to neutrals, ions, and electrons, 
respectively, VV is the number density of particles, e is the absolute value 
of the electronic charge and ¢,,. is a coefficient which describes the frictional 


force between species and is given by 
: 
2 in MM . 
(2x1 ume Oy (A3) 
NJ 7 
where k is Boltzmann's constant, 7 is the temperature, m is the particle 


mass and O,, is the elastic collision cross-section. 
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The reduction of the general equations to the form given in (A1) is 
not completely straightforward since it involves neglecting terms with 
gradients of the dynamic pressure of the ions. In the case being considered 
the discontinuous change in magnetic field produces a discontinuous change 
in ion velocity according to (Alb). In the immediate vicinity of this 
discontinuity the above equations are therefore not valid. The width of 
this region would, however, be of the order of one mean free path and, 
therefore, very small compared with the characteristic dimensions of the 
solenoid. Furthermore the total force due to the current resulting from the 
terms which have been dropped cannot be larger than the dynamic pressure 
of the ions. The equations (A1) are, therefore, valid everywhere except 
in the immediate vicinity of the discontinuity in magnetic field and the 
discrepancies there do not produce significant effects on the flow field as 
a whole. 

The generalized Ohm’s law given in (2.3) may be obtained by direct 
substitution of (A1b)in(Ala). The physical significance of the coefficients 
is, however, more easily described in terms of the equations (Al). By 
considering the case of zero magnetic field strength it is clear that o is the 
electrical conductivity in the absence of a magnetic field. From (A1b) 
it may be seen that the coefficient, A*, determines the magnitude of the slip 
velocity of the ions relative to the neutrals. The term involving «* gives 
rise to a current which is not parallel to the electric field when observed 
from a coordinate system moving with the ion velocity. Using terminology 
carried over from solid state physics this current is usually referred to 
as a Hall current. 

Using (A2) and (A3) it can easily be shown that the non-dimensional 
Hall coefficient «x, given by (2.8a), is the product of the angular frequency 
at which the electrons spiral around the magnetic field lines and the mean 
free time between collisions of electrons with ions or neutrals. In other 
words the Hall coefficient is 27 times the average number of revolutions 
an electron experiences between collisions. ‘The magnitude of the Hall 
coefficient depends, of course, on the physical conditions being considered. 
However, to obtain an estimate of the magnitude we will consider the case 
where the magnetic pressure B” 2u is equal to the gas pressure p*. This 
condition is equivalent to taking both S’R,, and p*/pg® of order unity. 
Using this condition and further assuming that the degree of ionization 
is so low that the second term in the denominator of equation (A2a) is 
negligible and taking O,, = 1 x 10- cm? one finds 

k = O-2y (ps,/p) (A4) 
where p,, is the sea-level density of air. From the above formula it is clear 
that the Hall currents must be considered in the density regions that are 
of interest for high speed flight applications of magneto-hydrodynamics. 

The ion slip coefficient A may be expressed in terms of the Hall coefficient 


as follows: 
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where the masses of the ions and neutrals are taken as the molecular weight 
of Ne and it is assumed that 0, 0, \ Si As pointed out by Patrick 
(1956), under typical flight conditions €,, €,, is of the order of 500. For 
degrees of ionization less than 2 x 10°? the term in brackets may therefore 
be taken as unity. 
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Influence of diaphragm opening time on shock-tube 
flows 


By DONALD R. WHITE 


General Electric Research Laboratory, Schenectady, New York 
(Received 13 May 1958) 


SUMMARY 

Shock waves generated in a shock tube by use of hydrogen 
or helium as a driver gas and air, nitrogen, oxygen or argon as 
a driven gas have higher velocities than predicted by simple theory 
when sufficiently large diaphragm pressure ratios are used. 
Expected shock-tube performance curves have been constructed 
using the equilibrium Hugoniot for the driven gas, both for the 
usual model of shock tube flow, which assumes instantaneous 
diaphragm removal, and for a suggested model based on a finite 
rupture time for the diaphragm. Agreement between experiment 
and the latter model is in general good, and the differences are 
qualitatively accounted for by the pressure waves expected to 
result from mixing between driver and driven gases at the contact 
surface. These waves may be either compression or expansion 
waves, depending on the relative heat capacities of the two gases. 
The maximum shock strength observed as a shock goes down the 
tube was found to occur at a distance from the diaphragm which 
increases with the shock strength, and the strongest shocks were 
found to be still accelerating at the end of a 42 ft. long shock tube 
of 3} in. square cross-section. Diaphragm breaking time has 
been measured and found to be consistent with the observations 
on the shock formation distance. 


1. INTRODUCTION 


A number of investigators have shown that calculations of expected 
shock-tube performance based upon the usual ideal model of shock-tube 
flow (i.e. no attenuation, instantaneous diaphragm removal, and one- 
dimensional flow) agree satisfactorily with experiment, the observed shock 
being slightly weaker than predicted for the given initial conditions. Even 
the etfect of a change in area at the diaphragm section has been shown to 
be amenable to a one-dimensional treatment (Alpher & White 1953). 
Most of the reported comparisons between expected and observed shock- 
tube performance, however, have concerned relatively weak shocks generated 
by diaphragm pressure ratios less than z= ),/p, = 10° (for example, 
Bleakney, Weimer & Fletcher 1949). With a z of the order of 10*, shocks 
stronger than predicted were observed and reported without comment by 
Yoler (1954) and Rabinowicz (1957). One might expect that the greater 
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velocity, temperature and pressure ratios encountered as stronger shock 
waves are generated would cause greater departures from the expected 
performance, and that study of these departures might illuminate the 
nature of real shock-tube flows. 

Previous investigators have demonstrated photographically (Geiger & 
Mautz 1949;Glass, Martin & Patterson 1953) that the process by which a 
diaphragm actually ruptures involves formation of a small opening near the 
centre, followed by the outward tearing of the material and the folding back 
of the leaves so formed. ‘The effective cross-sectional area of the shock tube 
at the diaphragm location therefore increases gradually with time instead 
of attaining its maximum value instantaneously as assumed in the usual 
ideal analysis. ‘lhe effective piston velocity pushing on the driven gas 
increases with the degree of opening of the diaphragm. Since each 
positive increment in this piston velocity sends out a compression wave, 
the gradual opening will generate a train of compression waves which must 
overtake each other to form a shock wave. Further, the high velocity jet 
which issues from the first hole in the diaphragm generates a high degree 
of turbulence as it penetrates the driven gas; and as the contact zone 
between the driver and driven gases moves down the tube, the initially 
gross mixing proceeds toward mixing on a microscopic scale. These two 
processes resulting from the gradual diaphragm removal, generation of 
a large amplitude steepening compression wave and the mixing of the 
driver and driven gases, will be examined more closely to determine what 


observable etiects might be expected to result. 


2. "THEORY 


Formation of a shock from a compression wave 

‘To estimate the result of a finite but unspecified opening time for the 
diaphragm, consider that an unsteady isentropic compression wave is first 
formed in the driven gas, and that the pressure profile of this compression 
wave then steepens as it progresses down the tube so that a shock wave is 
formed. Since the gas velocity is not the same behind a shock wave as 
behind a compression wave having the same pressure ratio, a third reflected 
wave is required to match the boundary conditions as the shock is formed. 
It is readily shown that this third wave must be an upstream-facing 
expansion wave, and that the transmitted shock wave has a lower pressure 
ratio than does the original compression. ‘To determine the expected 
shock-tube performance as a function of the diaphragm pressure ratio, 
the pressure and velocity of the expanded driver gas (ps, u,) will be equated 
to the pressure and velocity behind the compression wave (p,, u,.); and then 
the strength shock (py fp, or us a,) resulting from the steepening of this 
compression wave will be determined. ‘This model of the shock-tube flow 
is displayed in the (p,x,¢)-diagram of figure 1. Since further wave 
reflections and interactions are not taken into account, the strength computed 
on this basis must represent a transient condition, ‘The shock strength 
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should approach that predicted from the usual model (neglecting losses) 
as the effects of the starting process disappear. 

For simplicity in the analysis required to determine the shock formed 
from a given compression wave, consider that the individual elements 
of the compression all overtake the leading edge at the same point in the 
tube, resulting in a transmitted shock and a reflected expansion fan as 
shown in figure 1 and the insert of figure 2. Let the state of the gas in 
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Figure 1. Suggested model of shock-tube flows as depicted by a 
(p, x, t)-diagram. 


the various regions be indicated by the suffixes 1, 2, c, and e. In figure 2, 


for y, = 1-4, the curve labelled ‘isentropic compression’ represents the 
locus of states which may be reached through an isentropic unsteady 
compression wave from the initial state (p~,, “4, = 0) assuming ideal gas 
properties. Similarly, the curves labelled ‘shock states’ represent the 
locus of states (p., vu.) which may be reached from the initial state through 
a shock wave, assuming ideal gas properties or equilibrium air, as indicated. 
From any selected state (p,/p,) on the compression curve, the locus of 
states which may be reached through a backward-facing unsteady isentropic 
expansion wave (again assuming ideal gas properties) may be constructed 
as has been done to obtain the dashed lines labelled ‘expansions’. ‘The 
intersection of such an expansion curve with the shock curve determines 
that gas state which can satisfy the boundary conditions p, = p, and u, = Us, 
and hence defines the shock which will result from the steepening of the 
initial compression. 

When the driver and driven gases have been selected, and the sound 
velocities a, and a, and specific heat ratios y, and y, are known, then from 
the velocity ratio ua, behind the selected compression wave p,/p,, the 
value of uv, a, is determined by the requirement that the velocities be matched 
across the contact surface, u,=u,. From us/a, the expansion pressure 
ratio py ps for the driver gas is found, assuming as usual that u, = 0. The 
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diaphragm pressure ratio is then determined from (p,/p3)(p,/p,) =z. The 
procedure for obtaining the expected shock-tube performance curve 
based on the ideal model with instantaneous diaphragm removal is well 
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A graphical method for determination of the strength of the shock wave 


y, = 1-4. 


Figure 2. 
formed by steepening of an isentropic compression wave ; 


For the large diaphragm pressure ratios considered here, the gases 
cannot in general be considered ideal for either model of the shock-tube 


flow. Equilibrium shock Hugoniots for the shocked gases were used and 


were obtained from the following sources: for air, Gilmore (1955); for 
nitrogen and oxygen, Alpher & Greyber (1958); and for argon, Bond 
(1954) and unpublished calculations supplied by AVCO Research 
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Laboratory. In addition, hydrogen used as a driver gas expands to such 
low temperatures that changes in the degree of excitation of its internal 
degrees of freedom must be considered (Huber 1958). It may be noted 
that use of the equilibrium isentrope (Hilsenrath 1955) for the hydrogen 
expansion results in a decrease by about one Mach number in the predicted 
shock strengths for s ~ 10°, this decrease being greater for the ideal model 
because of the larger expansion ratio. 

Expected performance curves for the various gas combinations are given 
in figures 3 to 7, both for the ideal model and for the ‘formation-from- 
compression ’ model for a shock tube with an area reduction of 1-51 across 
the diaphragm section. ‘The method of Alpher & White (1958) is useful 
when equilibrium gas properties are to be incorporated. 

It is characteristic of the flow processes associated with this formation- 
from-compression model that the initial rate of extraction of energy from 
the driver gas is lower than that rate obtained using the ideal model. ‘That 
is, for a given diaphragm pressure ratio the driver gas expands through a 
smaller pressure ratio and acquires a lower velocity. Yet in some instances 
a stronger shock is predicted. ‘This seeming contradiction is associated 
with the transient nature of the flow being analysed. We have considered 
in this model an isentropic compression wave whose leading edge propagates 
at the velocity of sound a, in the low pressure gas. Although energy has 
been extracted from the driver at a lower rate, up to the point where the 
shock forms this energy has gone isentropically into a quantity of gas 
proportional to a,, whereas in the ideal model the energy goes anisentropically 
into a quantity of gas proportional to the product of the shock Mach number 
and a,. When the transmitted shock forms from this compression, it is 
maintained both by the energy still being taken from the driver gas through 
the first expansion wave (state + to state 3) and by the additional energy 
extracted from the isentropically compressed driven gas through the second 
upstream-facing expansion wave (state c to state e). This region of driven 
gas which has undergone an isentropic compression acts in a sense as an 
energy accumulator, feeding energy back to generate transiently a shock 
wave stronger than can be maintained by the energy being taken from the 
driver gas alone. 

This formation-from-compression model can be only an approximation, 
of course, because moderately strong shocks are found relatively near the 
diaphragm even when the maximum strength occurs much further down 
the tube. Since the entropy rise across a shock wave is proportional to 
the cube of the shock strength, however, this isentropic model should be 
a reasonable approximation even if, for example, a series of weaker shocks 
overtake each other to form a strong shock. 

Mixing of gases in the contact region 

It is clear that during every shock-tube experiment the driver and driven 
gases mix to some extent. While the magnitude of such mixing must be 
influenced by many factors, including the gases used and the shock strength, 
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a jagged opening projecting into the tube (the breaking diaphragm) would 
in particular promote mixing. While a quantitative prediction of the effect 
of such mixing on the shock strength cannot be made, nevertheless the 
direction of the etfect may be determined. 

In those experiments where the driver and driven gases have different 
heat capacities, this mixing if considered to be carried forward at constant 
pressure will be accompanied by a change in volume, since the gases will 
generally be at ditkerent temperatures. While this change in volume is 
taking place, there ts etfectively a velocity gradient across the region of 
mixing resulting in unequal velocities for the regions of uniform flow on 
either side. ‘l’o indicate the direction and possible magnitude of this change 
in volume, consider on either side of an ideal contact surface two equal 
volumes |} containing 7, moles of gas a at temperature 7, and n, moles 
of gas 6 at 7), respectively. Let \ T JT, =n,/n.. The change in 
volume after mixing has been completed at constant pressure is obtained 
from 

14 Al +V (ee) 
V 2 N+Cp, an 
and for V > 1, this becomes $(1+ Cp, Cp,). 

Considering only shock-tube experiments where the driver and driven 
gases are initially at the same temperature (not combustion driving), a 
region of the flow which includes the contact zone would be expected to 
undergo a volume increase trom the mixing of a monatomic driver gas 
and a diatomic driven gas. Conversely, a volume decrease would be 
expected when the heat capacity of the driven gas is smaller than that of 
the driver gas. ‘his mixing can be considered to affect the shock strength 
either through generation of pressure waves or through a change in the 
effective piston velocity. ‘To give a numerical example, in a shock tube 
of constant bore, a diaphragm pressure ratio of 104 between helium and air 
would produce on an ideal basis a Mach 7 shock. ‘The temperature of the 
expanded helium would be 38 K, and that of the shocked air 3140° kK. 
This yields V = 83, and Al) V = 0-194. The limiting value for large 
diaphragm pressure ratios for these gases is AV J = 0-20, 

The effect on the shock wave cannot be determined quantitatively, 
since one does not know either the extent of the mixing region or the rapidity 
with which mixing approaches completion ; and hence neither the magnitude 
of the compression waves nor the difference in flow velocity across the 
contact zone is known. It may further be pointed out that with the stronger 
shock waves the excitation of further internal degrees of freedom of the 
shocked gas results ina still greater heat capacity, and hence a greater value 
for the strong shock limit of AT}. Wray (1956) has reported observations 
on strong shock waves (., ~ 8) generated in fluorine by use of helium 
as a shock-tube driver gas. Fluorine is largely dissociated under these 
conditions and hence has a large heat capacity. He observed the shock 
velocity first to decrease with increasing distance from the diaphragm, 


reach a minimum, and then increase to the end of the tube. 
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Effect of diaphragm opening time on formation distance 

The expected effect of the diaphragm opening time upon the formation 
distance may also be estimated by use of a simple model. Assume first 
that in a shock tube of constant bore the flow at any time is that predicted 
by use of the ideal model, i.e. with instantaneous diaphragm removal. 
Let the gas combination and the diaphragm pressure ratio be given, and 
further assume ideal-gas behaviour for both driver and driven gases. 
‘Then one determines the distance travelled by the shock wave at the time 
it is overtaken by a sound wave originating from the diaphragm location 
at some given time after the flow was initiated. ‘This distance becomes 
relevant if one assumes that the maximum shock strength occurs as a result 
of the last compression sent out as the diaphragm is fully removed. 

Proceeding downstream from the diaphragm location, the flow may be 
divided into three regimes, an upstream-facing, isentropic unsteady 
expansion wave, a region of uniform flow of expanded driver gas, and 
a region of uniform flow of shocked gas. ‘The velocity of propagation of a 
sound wave will be the sum of the fluid velocity and the local velocity of 
sound. Assume for simplicity that the average velocity of propagation 
through the unsteady expansion wave is the mean of the velocities at the 
diaphragm location and at the foot of the expansion wave. All flow velocities, 
local sound velocities and propagation velocities may be readily calculated 
on this basis. A consequence of the highly idealized model utilized is that 
the distance travelled by the shock will be directly proportional to the time 
after flow initiation at which the signal is sent out. Representative results 


are given in table 1. 














| | 
[Diaphragm | Shock | Formation distance per 
Gases | pressure Mach | 100 microseconds effective 
| ratio number | diaphragm opening time 
a —|- | 
He/N, | 2500 6 | 6-5 ft. 
H./N, | 2500 8-8 | 8-3 ft. 
H./N, | 33 6 | 3-3 ft. | 


| 








‘Table 1. Estimation of formation distance for finite opening time of diaphragm. 


It may be noted from table 1 that, for the same shock strength, the shock 
should form in only half the distance with 2 hydrogen driver as compared 
with a helium driver. For the same diaphragm pressure ratio, the formation 
distance with the hydrogen driver is only 30°,, greater. ‘The use of a 
combustion driver would further shorten the formation distance for a 
given shock strength, because of the higher temperature and hence higher 
propagation velocity through the expanded driver gases. 
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3. EXPERIMENTS 
Apparatus and procedure 

These data were obtained in the 3-25 in. square shock tube at the 
General Electric Research Laboratory. The driver section of this tube is 
of 4:5 in. inside diameter and 6-5 ft. long, and the driven section is 3-25 in. 
square and 42 ft. long. Most of the diaphragms used in these experiments 
were made of type 302 stainless steel, 0-018 in. thick, scored in the form 
of an X by a milling cutter with a 0-010 in. radius tip. The groove depth 
was generally 0-008 or 0-009 in., and the breaking pressures varied from 
15 to 30atm. ‘The unsupported diameter was 4:5 in. 

For most of these experiments the initial pressure in the driver section 
was between 0-040 and 10 mm of mercury, and was recorded as the average 
of the readings of two McLeod gauges or, for the higher pressures, a 
precision Bourdon gauge. A calibrated Bourdon gauge was used to read 
the driver gas pressure as it was slowly increased to the point of rupture 
of the diaphragm. From these readings the pressure ratio across the 
diaphragm was obtained as the independent variable. 

The arrival of the shock wave at each of up to twelve stations along the 
tube caused a signal to be displayed on a cathode-ray oscilloscope which 
Was tracing a raster pattern at the rate of 20 microseconds/cm. <A photo- 
graphic record of the trace permitted determination of the variation of 
shock velocity with distance from the diaphragm. ‘The transducers 
employed in the tube were thin-film heat gauges, ionization gauges, or 
pressure gauges of our own design, depending upon the conditions of the 
experiment; and they were typically 32-6 in. apart. 

The average shock Mach number was computed for the interval between 
each pair of shock detection stations, and the maximum of this set of values 
plotted against the diaphragm pressure ratio. 


Shock tube performance curves 

Experimental shock-tube performance data together with predicted 
curves are shown in figure 3 for helium/nitrogen and hydrogen/nitrogen, 
and in figure 4 for helium/argon and hydrogen/argon. ‘The tendency for 
the experimental points to follow the formation-from-compression curve 1s 
clear, and in fact shocks stronger than predicted by the ideal model are 
noted for all gas combinations. For the case of helium/nitrogen in figure 3, 
as well as helium /oxygen in figure 5 and helium/air in figures 6 and 7, the 
shock strengths exceed even the stronger of the two predicted curves. 
These gas combinations are such that one would expect the contact zone 
mixing to enhance the shock strength. For hydrogen/argon, where the 
mixing would be expected to attenuate the shock, the experimental points 
are somewhat weaker relative to the compression-formation curve than are 
the hydrogen nitrogen data. For helium/argon and hydrogen/nitrogen, 
mixing should have no effect to the extent that further excitation of internal 


degrees of freedom may be neglected. 
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Figure 3. 
Compression waves due to contact zone mixing 
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Figure 4. Observed and predicted shock-tube performance for helium/argon 
Expansion waves due to contact zone mixing would weaken 


hydrogen/argon. 
the shock in hydrogen/argon experiments. 
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A set of experiments was run with helium oxygen to determine whether 
significant differences could be found relative to helium/nitrogen. For the 
strongest shocks (.V7, ~ 10), the dissociation of the oxygen is significant 
whereas it is still negligible for nitrogen; and the heat capacity, and hence 
the effect of mixing with the helium, would be greater for oxygen than for 
nitrogen. ‘The data are given in figure 5. Comparison with the 
helium/nitrogen results given in figure 3 shows that the helium/oxygen 


shocks appear to be somewhat stronger. 
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Figure 5. Observed and predicted shock-tube performance for helium/oxygen. 


The data shown in figure 6 were obtained in connection with work on 
the effect of an area ratio across the diaphragm section (Alpher & White 
1958). ‘These experiments show that the incremental shock strength due 
to a change in the area ratio across the diaphragm section is adequately 
predicted by a theory based on the ideal model even though the actual 
shock strength may be appreciably greater than predicted. 

Use of the equilibrium rather than the ideal-gas Hugoniot for the shocked 
gas results in the prediction of a weaker shock for a given diaphragm pressure 
ratio. Expected performance curves have been computed both on the basis 
of an ideal gas and of a real gas, and the comparison with experiment is 
shown in figure 7 for helium air. ‘The rather striking agreement with the 


real-gas compression-formation curve is probably partially fortuitous, in 
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Figure 6. The effect of an area ratio across the diaphragm section upon the 
maximum shock strength observed. 
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Figure 7. Comparison of experimental shock-tube performance with predicted 


performances based upon various assumptions. 
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the sense that apparently the enhancement effects of mixing approximately 
balanced the attenuation mechanisms which must always be present. 

In each of figures 3 to 7, the shock strengths observed for the largest 
diaphragm pressure ratios represent a lower bound on the maximum 
which could have occurred in the experiments, since the shock was still 
accelerating at the end of the tube where the recorded strength was measured. 


Figure 8. 'l'ypical records of shock strength at various distances from the diaphragm 
in the 3}-in. square shock tube. The horizontal lines indicate the average 
velocity of the shock wave over the corresponding distance. For the largest 
diaphragm pressure ratios the shock is still accelerating at the end of the tube. 

Further indication of the relative magnitudes of the mixing and 
formation-from-compression effects may be noted from the shock velocity 
records. In the hydrogen/argon experiments, for example, in addition 
to the losses due to the interaction between the shocked gas and the wall, 
the mixing at the contact surface would, for the weaker shocks, cause the 
shock to decelerate. ‘The attenuating effect of mixing would decrease, 
however, for the stronger shocks when the hydrogen is cooled to the point 
that its heat capacity approaches that of a monatomic gas, and when the 
ionization of the argon greatly increases its heat capacity. Representative 
data showing the variation of shock strength with distance along the tube 
are given in figure 8. The horizontal lines indicate the average velocity 
of the shock over the corresponding distance, and do not imply that the 
actual shock velocity varied in a stepwise fashion. For the larger diaphragm 
pressure ratios, the shock wave is observed to accelerate for the full length 
of the tube, emphasizing the dominant role played by these non-ideal 


processes. 
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From such velocity records the distance along the tube required for 
realization of the maximum shock strength has been determined. ‘This 
distance is plotted vs shock Mach number in figure 9, and vs diaphragm 
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Figure 9. Shock-formation distance as a function of shock Mach number in the 

34-in. square shock tube. Points for helium/argon and hydrogen/nitrogen 
generally lie between the helium/nitrogen and the hydrogen/argon values. 
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Figure 10. Shock-formation distance as a function of diaphragm pressure ratio in 
the 3}-in. square shock tube, Points for helium/argon and hydrogen/nitrogen 
are similarly distributed. The solid point was obtained by loading each of 
the four leaves of the scored diaphragm with 6-2 gm of solder. 


pressure ratio in figure 10. The curves probably apply only to this particular 
tube and to these diaphragms, and should be used only qualitatively in 
estimating the expected performance of other shock tubes, 
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Diaphragm opening time 

To determine experimentally the time required for the leaves of the 
diaphragm to fold back to the full open position, the transmission of light 
through the diaphragm section was measured as a function of time. A light 
source was mounted near a 3 in. window in a blind flange at the end of the 
driven section, 42 ft. from the diaphragm. A barium titanate pressure 
transducer was mounted in the driver chamber wall near the diaphragm 
section and responded to the metal vibration as the diaphragm first began 
to break under the gradually increasing pressure, producing a signal to 
trigger a cathode-ray oscillograph. A photomultiplier viewing the side 
wall of the driven chamber observed the variation with time of the intensity 
of the scattered light and displayed it on the oscillograph. 

‘The diaphragms used in these experiments took about 600 microseconds 
to open fully. ‘The initial rate of opening was low however, since the light 
intensity recorded had reached less than 10°,, of its final value after 
200 microseconds. ‘his slow initial rate of opening was also characteristic 
of both 0-010 in. and 0-035 in. thick diaphragms, which took somewhat 
shorter and longer times respectively to open fully. 

These measurements of formation distance and diaphragm opening 
times are considered to be consistent with the highly simplified analysis 
previously described, since the many approximations involved preclude 
other than qualitative interpretation of the results. For example, a 
convergence at the diaphragm section would tend to shorten the formation 
distance because of the higher temperature of the expanded driver gas, 
and attenuating effects would do likewise by inhibiting attainment of the 
maximum possible shock strength. 

However, one may conclude that an opening time of several hundred 
microseconds should result in an appreciable formation distance for the 
larger diaphragm pressure ratios, and that this distance should be more 
nearly independent of diaphragm pressure ratio than of shock Mach 
number as various gases are used. This is observed in figures 9 


and 10. 


Shock attenuation 

Since these stronger-than-predicted shock waves (as in helium/ nitrogen) 
exist only transiently, successive wave interactions may be expected to 
reduce the shock strength. Mixing in the contact zone may either have 
no effect or act to enhance or reduce the shock strength, depending upon 
the gases employed. ‘Therefore, these two effects associated with the 
starting of the shock tube flow may both strengthen and ‘attenuate’ the 
shock at different times in the same experiment. Interpretation of 
attenuation data from strong-shock experiments in terms of energy loss 
mechanisms must be done carefully, since these other adiabatic processes 


may also cause the shock strength to vary with time, 
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4. CONCLUSIONS 

The finite time required for removal of the diaphragm can result in 
an appreciable modification of the shock-tube flow for the case of the 
larger diaphragm pressure ratios. For a monatomic driver/diatomic driven 
gas combination, this opening time can cause generation of a shock signifi- 
cantly stronger than predicted by the usual model. An analysis based 
upon the formation first of a compression wave, then of a shock wave, 
agrees with experiment. ‘lhe mixing of the driver and driven gas in the 
contact zone may also modify the shock strength, either strengthening or 
weakening it according to the relative heat capacities of the gases used. 
Observed shock deceleration may result partially from these adiabatic 
processes of mixing and wave interactions, as well as from attenuating 
effects of energy losses to the wall. 

The shock-tube starting process can occupy a significant portion of 
the length of the tube. In these experiments, for the largest diaphragm 
pressure ratios (s ~ 10%) the shock was still accelerating 40 ft. from the 
diaphragm. 

The one-dimensional analysis of Alpher & White (1958) successfully 
predicts the shock enhancement due to area ratio at the diaphragm, even 
for the larger diaphragm pressure ratios where the actual shock strength 
may not be predicted by use of the ideal model. 

In experiments employing combustion driving, approximately equivalent 
performance curves are computed from the ideal model, whether hydrogen 
or helium is employed as a diluent. ‘The processes discussed here, however, 
might cause one or the other to be advantageous in a particular experiment. 

This research was supported by the Ballistic Missile Division of the 
United States Air Force under contract AF-04-(645)-24. A preliminary 
account of the work was given at the annual meeting of the American 
Physical Society in New York, January 1957. 
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two-dimensional steady flow 
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SUMMARY 
A set of two-dimensional subsonic flows past certain cylinders 
is obtained using hodograph methods, in which the true 
pressure-volume relationship is replaced by various straight-line 
approximations. It is found that the approximation obtained by 
a least-squares method possibly gives best results. Comparison 
is made with values obtained by using the von Karman-Tsien 
approximation and also with results obtained by the variational 
approach of Lush & Cherry (1956). 


1. INTRODUCTION 

There are a number of ways available for finding the approximate 
steady two-dimensional flow past a fixed cylinder. One such method is 
the tangent-gas method proposed by Chaplygin and later extended by 
Tsien, which uses a tangent to the curve, pressure (p) vs volume (zv), as a 
representation of the true relationship. ‘The significance of this device 
is that the hodograph equations can be reduced to the Cauchy-Riemann 
differential equations. Demtchenko (1932) considered a tangent at the 
point on the curve corresponding to the stagnation conditions, and as a 
result his theory can only be applied, with any accuracy, to flows with 
speeds up to about one-half sonic speed. ‘Tsien (1939), on von Karman’s 
suggestion, took the tangent at the point corresponding to conditions at 
infinity, assuming uniform flow there, and thus succeeded in extending 
the range of accuracy. ‘These cases are only two examples of a set of flows 
obtainable by taking a tangent at each point of the (p, v)-curve, or, more 
precisely, a series of straight-line approximations to the curve within the 
significant limits for the particular problem. 

A method is suggested for finding the particular line which possibly will 
yield the best results, and comparison is made with the accurate values 
obtained by Lush & Cherry (1956), using a variational method, for the 
adiabatic flow past a circular cylinder. 


2. EQUATIONS OF MOTION 
A 


ume 


DQ 
nD 


p = a-bv (1) 
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as the straight-line approximation to the (pf, v)-curve, which for illustration 
may be taken to be the adiabatic relation 
po’ = 1. 
(For convenience the stagnation pressure and density are taken as unity.) 
Since v= 1/p (p representing fluid density), Bernoulli’s theorem 
immediately yields 
g? = b%(e2— 1), (2) 

where q represents fluid speed such that v = 1 when g = 0. 

If ¢ represents the velocity potential and % the stream function, the 
hodograph equations are 


TC 


where @ is the angle the fluid velocity makes with a fixed direction. 
Equations (3) can be reduced immediately to the Cauchy-Riemann 
differential equations 


d., alae Po, db, = Pas (4) 
by using (2), and setting 
w= | (p/q)dq 
= log[ Bq {b+ (0+ 4)'*)], (5) 


where B is an arbitrary constant. 
From (4), 6+7 must be an analytic function of w—16, that is to say, 
W = +i = f(4), 
with o = w +76 and the over-bar representing a complex conjugate. 
If x,v represent coordinates in the plane of flow, then quite simply 
dz = q-'e!°(df+ip— dy), z= x+1y, 
and, in particular, on a profile given by % = constant, 
dz = q"e"® dW. (6) 
Now the function f(¢) will represent the flow of an incompressible 
fluid in a certain ¢-plane, so that 
dW , ; 
—— = Oe-*, O = Ce”, 
dl 7 
where C is an arbitrary constant and QO, @ are respectively the magnitude 
and angle of inclination of the velocity vector of the incompressible fluid 
in the ¢-plane. Equation (5) shows that 


O = 2bq/{b+ (b? + q?)*}, (7) 
if the condition O - g as g — 0 is assumed, so that 
q = 4b*O/(4b? — O?). (8) 


\ssume the flow in the ¢-plane to be uniform and parallel at infinity, 
and write W = O,,G(¢), O, being the fluid speed at infinity; then 


F.M. 2Q 
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and from (6) and (8), 
. O08, :/dG\? |, 
eat ee (3) df, (9) 
since dG = dG along the profile % = constant. 


3. PRESSURE AND DENSITY 
Let the straight line (1) be parallel to the tangent to the curve pv’ = | 
at the point (p,,7,). Then 


dp ; (? dp Per 


where c, is the local sonic speed. This linear assumption yields 


 -  - 
c=: —_- = —. 
dp p* 
so that b* = cp? = c pi’. (10) 
Equation (2) now gives 
p? = 1— M2, (11) 


where VW = q/c is the Mach number. ‘This means that the density and 
pressure are given very simply in terms of .V/, and as p = 0 when M = 1, 
it appears that (11) will probably not be a good approximation near the 
critical velocity. 

I . . 7 : 3 . , : ~ 

It seems a reasonable assumption therefore to find only the fluid speed 
oy this approximation method, and to use the correct adiabatic relations 
by tl ) t thod 1 t tl rrect adiabatic relat 


y b \ 7-2 y—1 ,\Mr- ‘ 
a (1 aa ) ; ei = #) (12) 


for pressure and density once q has been determined. 





4. BEST APPROXIMATION WHEN y = 1-4 
Obviously a set of flows can be constructed by taking values of 6 to 
correspond to tangents at various points on the (p,v)-curve between the 
stagnation and critical conditions. Such a set has been calculated by taking 
tangents at points on the curve corresponding to Mach numbers 0, 0-1, 
()-2, ..., 1-0, and applied in each case to uniform flow past a circular cylinder 
of unit radius in the ¢-plane, so that 


G(t) = C4 1/6. 


4 


‘The Mach number at infinity in the z-plane is taken to be WM, = 0-35, 
and since 


g = yM?/{1 4+ 3(y—1)M*, 
g.. = 0-4091, when y has the value 1-4. 

Table 1 shows the thickness ratio of the corresponding cylinder in the 
s-plane, also the ratio qmax/¢:, Ymax being the maximum speed on the 
profile, which occurs at the points where thickness is greatest. An inspection 
of this table shows that although thickness ratio increases with MM, 











Tangent-gas approximation for steady flow 603 


Ymax/Y increases even more quickly, so that it is reasonable to suppose 
that for any given thickness ratio, the value of gmax/qg. is too small if the 
tangent corresponding to M = 0 is taken and too large if the tangent 
corresponding to M = 1 is used. The problem is now to select the ‘best’ 
tangent, or in effect the best value of 6*, since only the slope of the 

















| | 

‘ Thickness 

| M o* 0. ratio Fmax/ Ts | 

0 1-400 0:3975 1-041 27191 | 
0-1 1-383 0-3974 1-041 2-194 
0:2 1-334 0-3970 1-042 2-201 
0-3 1-257 0:3964 1-045 2-215 
0-4 1-144 0-3953 1-050 2-239 
0-5 1-044 0:3941 1-055 2-262 
| (6 (0)-9224 0:3921 1-062 2-300 

| 0-7 0-7989 0-3894 1-072 2-35 

r (08 0:6797 0:3865 1-086 2:423 
| 0-9 0-5798 0-3832 1:101 | 2-508 
| 1-0 0-4687 03779 | 1-127. | 2-663 








Table 1. Thickness ratio and g,,,,/¢, for various tangent approximations with y = 1-4. 











Figure 1. Approximation to (p, v)-curve. 


2Q2 
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straight-line approximation occurs in the results. Tsien (1939) chose the 
tangent at the point on the (p,v)-curve corresponding to the conditions at 
infinity, and it appears that the speeds obtained by this assumption are too 
small. It is suggested here that the value of 5? should be taken to 
correspond to the best straight-line fit to the (p, v)-curve, found, say, by the 
method of least squares, between the significant limits for the particular 
problem under consideration. For example, it is found by using the 
least-squares method over 11 points between the stagnation and critical 
points, that the best straight-line fit over this range corresponds to 
a = 1:7544, b? = 0:7982. This is shown in figure 1. 


5. CIRCULAR CYLINDER 
In order to test the accuracy of this method, the compressible flow, 
uniform at infinity, past a circular cylinder of thickness ratio unity is 
considered, and comparison made with the accurate results obtained by 


Lush & Cherry (1956). 























| t | x y q/45 (x?+ y?)! 2 
| eset = is 
| 
0 | 4-325 0 0 4-325 
10 | 4-247 0-811 0:3424 4-324 
| 20 4-022 1-590 0-6846 4-325 
30 3-661 2-298 1-0231 4-322 
40 3°191 2915. - || 13506 4-322 
50 2-643 3-423 | = 1-6589 4-325 
60 2-028 33820 | 19275 4-325 
70 1-370 4-103 | 237i 4-326 
80 | 0692 | 4:270 gew2-2 4-326 
90 0 | 4325 " 2-3981 4-325 
| | | 
Table 2. Coordinates and fluid speed for a nearly circular cylinder at 7, = 0-4. 


‘lo obtain such a flow in the s-plane, the uniform flow past an elliptic 
cylinder in the incompressible plane is considered, represented by 
G(L) = C +020, 
where € = ¢’+1/¢’; the profile is given by ¢’ = xe“. This is substituted 
in (9), from which it is found on integrating and separating into real and 
imaginary parts, 


v (- re Jeos A{ 20 +2?)cos t+ }(a?— 1)? log 


es 


(a?-—1)? eet || 


x? — 1) 
(a? + 2xcost+1)? | 





} 
| 
| 
{ 
S 


ae 
y= (2- -)sin , Af a( —g?)sint + (22-1)? tan-! era | | 
- = | 

(13) 





with A = Q? /45". 
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The semi-axes of the resulting cylinder can be found by setting ¢ equal 
to 0 and 37 in (13), and the ratio of these gives the thickness ratio § as 


9 9 2% } 

_— 1+a) —28+ Ja(o? — 1)tan( > ) | | 

o= 9 | é r fr. 
i a 








—_— a ee ! | tat! % 


C 3O oe) 9c 





oO 








Figure 2. qg/q, plotted against tan-'(y/x) for circular profile at M,, = 0-4. 
Comparison is shown with the von Karman-Tsien approximation and 
Lush & Cherry’s S, result. 


If now M,, = 0-4, then q,, = 0-4659, O,, = 0-4380, and, by employing 
the best straight-line fit previously found with 5? = 0-7982, A takes the 
value 0-0601. Using this value of A, it is found that when «? = 24-00, 
5 = 1 to four significant figures. ‘Table 2 shows coordinates, calculated 
from (13), and fluid speeds on the profile for various values of t. This table 
shows that the profile is almost circular, the error in radius being nowhere 
greater than 0-1°,. 








606 G. Power and P. Smith 


Figure 2 shows the graph of g/q, plotted against tan-'(y/x). From 
this an interesting comparison can be made with the results obtained by 
Lush & Cherry, and in particular with their S, approximation. Comparison 
is also made with a similar approach using the von Karman-'T’sien method. 
Agreement with Lush & Cherry is extremely good, although it must be 
observed that they use y = 1-405. 

It would therefore appear that by this means of approximation a very 
accurate result can be obtained, at least for flow past cylinders with 
elliptic-type cross-section. Unfortunately there are very few solutions 
that are known to be sufficiently accurate for further comparison to be made. 


One of the authors (P.S.) wishes to acknowledge the receipt of a 
maintenance grant from the Department of Scientific and Industrial 
Research. 
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Surface waves on water of non-uniform depth 


By JOSEPH B. KELLER 
Institute of Mathematical Sciences, New York University 


(Received 20 March 1958) 


1. INTRODUCTION 

Gravity waves occur on the surface of a liquid such as water, and the 
manner in which they propagate depends upon its depth. Although this 
dependence is described in principle by the equations of the ‘exact linear 
theory’ of surface waves, these equations have not been solved except in 
some special cases. ‘Therefore, oceanographers have been unable to use 
the theory to describe surface wave propagation in water whose depth 
varies ina general way. Instead they have employed a simplified geometrical 
optics theory for this purpose (see, for example, Sverdrup & Munk (1944)). 
It has been used very successfully, and consequently various attempts, 
only partially successful, have been made to deduce it from the exact linear 
theory. Itis the purpose of this article to present a derivation which appears 
to be satisfactory and which also yields corrections to the geometrical 
optics theory. 

The present derivation is based upon an asymptotic solution of the 
equations of the exact linear theory for periodic waves in water of arbitrary 
non-uniform depth. ‘The solution is asymptotic in the sense that the depth 
and wavelength must be small compared to the horizontal scale length of 
the bottom contour. Our solution is formal since we do not prove its 
asymptotic character. ‘Therefore, it is significant that the first term agrees 
exactly with the asymptotic form of the solution for waves in water with a 
uniformly sloping bottom, as the bottom slope tends to zero. ‘This 
asymptotic form was deduced by Friedrichs (1948). The accuracy of our 
result is indicated by the fact that for a bottom slope angle of 6°, the 
asymptotic form agrees very well with the exact solution at all points beyond 
one wavelength from the shore (see Stoker 1947). 

Previously Lowell (1949) has shown that in shallow water the geometrical 
optics theory follows from the linear shallow-water theory. The latter is 
a simplified form of the exact linear theory which applies when the depth 
is small compared to the wavelength. ‘The significance of this restriction 
is not clear, in view of the fact that Lowell’s derivation depends upon letting 
the wavelength tend to zero. Our derivation seems to clarify this point 
since we require the depth to tend to zero in a fixed ratio to the wavelength, 
while the horizontal scale of the bottom contour is constant. We could 
have obtained the same results by keeping the depth and wavelength constant 
while permitting the horizontal scale length of the bottom contour to become 
infinite. 








608 Joseph B. Keller 


A previous attempt by the present author to derive the geometrical 
optics theory from the exact linear theory was based upon an asymptotic 
solution obtained by requiring the wavelength alone to tend to zero. The 
resulting geometrical optics theory was independent of the depth. That 
result can be explained by noting that waves disturb the water only within 
a few wavelengths of the surface. As the wavelength tends to zero, the wave 
motion is confined to a vanishingly thin layer near the surface, so the waves 
do not ‘feel’ the bottom. This derivation, which may be useful for other 
purposes, is contained in Appendix I. 

We will now describe the geometrical optics theory to which we have 
already referred. In this theory a propagation velocity is defined at each 
point P of the surface. This velocity is just that which waves of the given 
period would have in water having uniform depth equal to that at P. By 
using this propagation velocity and employing Fermat’s principle of optics, 
rays can be defined. It is assumed that surface waves propagate along these 
rays. ‘The variation of the amplitude of these waves along a ray is then 
determined by using the principle of conservation of energy. In its optical 
form this principle states that the flux of energy is the same at all cross- 
sections of a tube of rays. In the present case a tube is bounded by two 
neighbouring rays. ‘The energy flux is proportional to the square of the 
amplitude of the waves and to the distance between the rays. It is assumed 
that at each point P the proportionality factor, which depends upon the 
depth, is the same as for plane waves in water of uniform depth equal to 
the depth at P. From these assumptions the variations of amplitude along 
a ray can be determined. ‘This amplitude variation has been studied by 
Miche (1944). 

2. FORMULATION OF THE PROBLEM 

We consider the periodic irrotational motion of an inviscid incompressible 
liquid bounded below by the rigid surface Y = —H(x,z) and bounded 
above by an unknown free surface. In the absence of motion the free 
surface is the plane Y = 0; when motion of angular frequency w occurs, 
itis Y = A#fe'n(x,z)}. Here t denotes the time, 7 is the complex amplitude 
of the surface wave motion, and the positive Y-axis points vertically upward. 
In the exact linear theory of surface waves, the wave amplitude 7 is given 
by 

n(x, 2) = (tw/g)@(x, 0,2). (1) 


Here (x, Y,z) is the velocity potential of the fluid motion, and g is the 
acceleration of gravity. 
The velocity potential satisfies the following conditions (see Stoker 
1947) 
A®=0 in 02 Y2 -A(x,2), 


= om Y= (B = w?'g), ) 
®,.+H,9,+H,®,=0 on Y= —H(x,2). (4) 


Furthermore, © satisfies other conditions which we will consider later. 
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We now introduce the new quantities y, h and ¢ defined by 
y = BY, h = BH, $(x, y,2) = O(x, Y,2). (5) 
In terms of these quantities, (2), (3) and (4) become 
Bo yyt+brr+%:2=90 in ODy> —A(x,2), ( 
d,=¢ on y=), ( 
B'6,+h,o,+h.¢,=0 on y= —h(x,2). (8 
( 


for large values of f. 

As the first step towards finding solutions, we introduce three unknown 
functions k(x,z), S(x,z) and A(x,y,z,8). Without loss of generality we 
may express ¢ in terms of these functions by the relation 

¢ = Acosh[k(y +h) ]e’*. (9) 

Next we insert (9) into (6), (7) and (8). These equations then become 
the following, in which « =k(y +h) and V = (0/éx, 3/dz): 

B?[{k? — (VS)*}4 cosha+A,,, 

+1B[(V2S)A cosha+2VS.V(A cosha)]+V?(A cosha) =0; (10) 

A,coshkh+kA sinhkh = A coshkh y= 0; (11) 

B?A,+1BAVh.VS+Vh.VA =0 y= —h. (12) 

We now want solutions of (10), (11) and (12) for large 8 which represent 

waves travelling on the surface. 


cosh « + 2k4 , sinh «] + 


3. SOLUTION OF THE PROBLEM 

To solve (10), (11) and (12), we assume that, for large 8, A has an 

asymptotic expansion of the form 
A ~ A,(x,z)+ > A, (x,y, 2)/(iB)”. (13) 
it l 

Then we insert (13) into (10), assuming that termwise differentiation 
is correct. When we equate coefficients of each power of 8, we obtain the 
following equations: 
Ps (14) 
S.V(A,cosha)+ AycoshaV2S; (15) 
S.V(A,,_,coshx)+A,,_,coshaV2S + 


(VS)? =k 
(A,),,cosha+2k(A,), sinha = 2 
(A,,),,, cosh « + 2k(A,,),, sinh x = 2 


+ V*(4,,_,.cosha«) (n>2). (16) 
In a similar way, we obtain from (11) the equations 
Rtanh kh = 1, (17) 
(A4,), = 0 y= (mn > 1). (18) 
Finally, from (12), we obtain 
(A,), = Ag VA.VS y= —-A, (19) 


(A,), = 4A,-1 Vk. VS+Vh.VA, 2 y=-—h (n>2). (20) 
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Let us now examine the equations (14) to (20) which must be satisfied 
by k, S and the 4,. Equation (17) is a transcendental equation which 
determines k(x, 2) in terms of the known depth A(x,z). ‘This 1s the same 
equation that occurs in the case of uniform depth and its solution is well 
known. Once k is known, (14) is the eiconal equation of geometrical optics 
for S. If the value of S is prescribed at each point of any given curve, the 
solution S can be found by means of the characteristics of (14). These 
characteristics are just the optical rays in a medium of refractive index k, 
and the surfaces S = const. are the corresponding wave-fronts. In terms 
of these rays the directional derivative AS . Vf, which appears in (15) and (16) 
is just k(df/d7) where + measures arc-length along a ray. 

Next let us consider the amplitudes 4, which satisfy (15) and (16). 
These equations can be simplified by means of the identity 


(4,,),, cosh x + 2k(4,,)ysinh« = (coshz)~'[(4,,),, cosh? «],,. (21) 
Now (21) is inserted into (15), which becomes 
[((A,), cosh? «],, = (2VS.VA, +A, V2S)cosh?« +A) VS.Vcosh?«. (22) 
Now we integrate (22) from zero to y, noting that, by (18), (4,), = 0 at 
y=0. Thus (22) yields 
(A,),cosh?a = (2VS.VA,+ A, V?2S+ 4, VS.V) | cosh? dy. (23) 
/0 
The integral in (23) is given by 
| cosh?k(y +h)dy = jk" 
< fsinhk(y +h) coshk(y +h)—sinhkhcoshkh]+}y. (24) 
We now use (24) in (23) and obtain 
(A,), cosh?« = 1(2VS.VA,+A, V2S +A, VS.V) x 
< {k-"[sinhacoshz—sinhkhcoshkh]+y}. (25) 
From (25) we can calculate A, up to an additive function of x and z, once 
Ay 


determine this additive function. For that purpose the equation (16) for 
n = 2 would have to be examined. 
To find A, we set y = —A in (25). Then we can eliminate (4,) 


means of (19). In this way we obtain the following equation for A): 


24,VA.VS = —(2VS.VA, + AyV2S)(sinh? kh +h) + AypVS.(VA— Vsinh? kh). 
(26) 


is known. The method which we will now use to find A, will also 


, by 


This can be written as 
VS.V{[A?(sinh? kh + h)] + [A?(sinh? kh +h)]V2S = 0. (27) 


Since V.S.V = k(d/dr), (27) is an ordinary differential equation along a 
ray. Its solution is 
RAV2S “ (28) 


A?(sinh? kh +h) = [A?(sinh* kh +h)], exp< — | 
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Luneberg (1944) has shown that the exponential in (28) is given by 


k(t) do(7 9) 

k(t) do(r) 
In (19), do(7,))/do(r) is the ratio of the width of a narrow strip of rays at 7, 
to its width at 7, or rather the limit of this ratio as the width tends to zero. 
When (29) is used in (28), the equation for the variation of A, along a ray 
finally becomes 








exp - is k2AV2S ar} - (29) 


A?(sinh? kh + h)k do = const. (30) 

Equation (30) has a simple physical interpretation. Except for a constant 

factor, the left side of (30) is the energy flux per length do of a plane wave 

in water of uniform depth h. ‘Thus (30) expresses the fact that the energy 
flux is constant along a tube of rays. 


4. CONCLUSION 

Let us now summarize our results. ‘The wave amplitude 7 is given by (1), 
the potential is given by (9), & is given by (17), S is determined by the 
eiconal equation (14), and A, is determined by (30). This yields the leading 
term in the solution for large 8. It conforms with all the principles of the 
geometrical optics theory and thus provides a derivation of that theory. 
Additional terms are given by the 4,,, which satisfy (16) and the boundary 
conditions (18) and (20). 

In Appendix II our result is compared with the asymptotic form of the 
exact solution for plane waves normally incident on a uniformly sloping 
beach. ‘That asymptotic form is obtained from the work of Friedrichs 
(1948). It agrees exactly with our result. 

Our result is not valid at caustics of the ray system since do is zero on 
them, and then (30) shows that A, is infinite. A method for obtaining 
the correct finite result at and near caustics was given by Kay & Keller (1954). 
At a shore line where / = 0, our result also fails because it yields an infinite 
value for A,. In this case the solution for the uniformly sloping beach 
given by Peters (1952) can be used near the shoreline to yield the correct 
result there. Both of the exceptions just described are instances of 
‘boundary layers’ in which the asymptotic behaviour of the solution is 
different from that outside these layers. 


APPENDIX I. ANOTHER ASYMPTOTIC SOLUTION 
We will try to solve (2), (3) and (4) for ® by writing 


® = ue'4S, (A1) 
In terms of u and S, (2) and (3) become 
82(VS)2u—iB[2VS. Vu+uV2S] + V2u = 0, (A2) 


BS ac — 
iBS,,ut+u, = Bu y= 0. (A3) 
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Here V is the three-dimensional gradient and y = Y is used for notational 
convenience. We now assume that u has the asymptotic expansion 


x 


u~ y u;(x,y,2)/(¢8). (A4) 
j=1 
Upon inserting (A4) into (A2) and equating coefficients of 8, we obtain 
(VS)? = 0, (A 5) 
2VS.Vu; +u,V?S = —Vu;_, j29, u_,=0. (A6) 
Similarly from (A 3) we obtain 
is, = 4 y = 0, (A7) 
(u;), = 0 y =0, J 29%. (A8) 
From (A5) we conclude that S is complex, so we write 
S= Rd: (A9) 
where R and / are real. Now (A5) yields the two equations 
VR.VI =0, (VRP = (VIFF. (A 10), (A11) 


Equation (A 10) shows that the wave-fronts R = const. are orthogonal to 
the surfaces J = const., which are essentially surfaces of constant amplitude. 
This is common in evanescent waves. 

When (A9) is used in (A7), there results 


R, = 0, i= —1, y= ©. (A 12) 
To satisfy (A 10) and (A 12), we assume 
I y, R = R(x, 2). (A 13), (A14) 
Then (A11) yields the eiconal equation for R 
(VRP = 1. (A 15) 


To solve (A6) with j = 0 for uy, we assume that uw, is real. Then (A6) 
implies that uy = uy(x, 2), and that 
2VR. Vig + uy V?2R = 0. (A 16) 
As before, (A16) is an ordinary differential equation along a ray. Its 
solution is 


Uy (7) = Up(7)[da(7) do(r)]' ? (A 17) 


From (A 15) it follows that the rays are straight lines, and along any ray 
R= R(t))+7. Then do(r,)/do(r) = (p+79)/(p +7), where p denotes the 
radius of curvature of the wave-front at 7,. 

Upon collecting our results, we find 


1) ac Ar 125 Br+By. (A 18) 


In writing (A 18) we have collected the factors independent of 7 into the 
factor A and shifted the origin of 7 so that p+ 7 becomes 7. In other words, 
we now measure 7 from the centre of curvature of the wave-front—i.e. 
from the caustic—on each ray. ‘The result (A 18) does indeed conform 
to the principles of geometrical optics, but the index of refraction to which 
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it corresponds is unity. The waves travel with the velocity w/8 = g/w, 
which is the velocity of waves of frequency w in water of infinite depth. 
This was explained in the Introduction. 

In constructing the solution (A 18) the condition (4) at the bottom was 
ignored and the result in (A 18) does not satisfy it. It is also impossible to 
satisfy this condition by including additional terms of the sum (A4). 
Therefore, the present expansion applies in water of infinite depth, when 
condition (4) is absent. When the depth is finite we must modify the 
expression (A1) by including another term of the same form, which we 
may think of as arising by reflection of the original term from the bottom. 
This reflected term must be so chosen that the two terms together satisfy 
condition (4) at the bottom. But then the reflected term violates the 
condition (3) at the top, which the first term alone satisfies. We therefore 
add another term, reflected from the top, to remedy this, then another 
reflected from the bottom, etc. 

On the basis of these considerations, we assume that ® has the following 
asymptotic expansion 


© ~ S exp(igs,)pe™ ¥ u,,/(iB). (A19) 


Here [n/2] denotes the greatest integer not exceeding n/2, and the S,, and u,, 
are functions of x, y and z. The terms with n = 0 are those considered 
before. By proceeding as before we can obtain equations for all the new 
terms. It appears that at y = 0 all the S, with m > 1 will have positive 
real parts which depend upon the depth. ‘Therefore they will all be of 
exponentially lower order in f than the terms with »=0 which we 
considered before. ‘Therefore, our previous result will still determine the 
surface shape or wave pattern. 


APPENDIX II. COMPARISON WITH AN EXACT SOLUTION 
Friedrichs (1948) has obtained the asymptotic form of the exact solution 
for waves on a sloping beach as the slope angle tends to zero. He considered 
only points on the top surface y = 0. However, it can be seen that his 
equation (16) also holds for y < 0 provided that xis replaced by x + zy in (12). 
If we then treat y as small compared to x, we obtain 
A(x +1v) = A(x) +1wydA'(x)+..., 
and 
R[A(x+1v)] = R[A(x)] +iwy/A(x) +... 
When these results are inserted into his (16) and a typographical error in 
the second exponent is corrected, (16) becomes 
XP (x,y, t) = —(1—A)A(A)exp[i(w 1k +t + $7) —y/A] — 
—~(1+A)A(A)exp[—i(wtk+t+}7)4+v/A]+.... 
Then the real potential function ¢ is given by 
d = AXP) = — A(A)[(1—A)e—4* + (1 +A)e”*]cos(w tk +t +37) 
= —2A(X)[cosh(y/A) +Asinh(y/A)]cos(wk +t +47). 
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Now we reintroduce dimensional variables by replacing x by Bx, t by ot 
(o is the angular frequency and 8 = o/g) and we define the depth A = Bwx 
(w is the slope of the bottom). We also introduce k’ = A. Then (12) 
becomes k’ tanhk’h = 1. Now ¢ is given by 
cosh[k'(y+h)] 

cosh k’'h 
From the definition of A(A), it follows that 4(A) = A, coshk’h, where A, 
is given by our equation (30). Thus, ¢ finally becomes 

d - 2A, cosh[k'(y +h) ]cos(wk + ot +47). 


This is exactly the sum of two expressions of the ferm (9) given by our 


@ = —2A(d) cos(w1k+ot +47). 





asymptotic solution. 


‘This paper represents results obtained under the sponsorship of the 
Office of Naval Research, Contract Nonr—285(06). 
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SUMMARY 

The effects of compresstbility on a radial laminar wall jet are 
investigated. On the assumption that the coefficient of viscosity 
is proportional to the temperature, it is shown that a similarity 
solution for the velocity distribution exists, which is expressible 
directly in terms of the corresponding solution for an incom- 
pressible wall jet. For arbitrary Prandtl number the energy 
equation is studied in detail and solutions are obtained for a variety 
of temperature conditions. 


1. INTRODUCTION 

‘The term ‘wall jet’ was introduced by Glauert (1956) to describe the 
flow due to a jet of air spreading out over a plane surface, either radially or 
in two dimensions. With compressibility neglected, Glauert studied the 
velocity distribution in a wall jet for both laminar and fully turbulent flow, 
and for laminar flow found an exact solution of the boundary-layer equations 
in the form of a similarity solution. This solution should be asymptotically 
approached whatever the initial form of the jet near the axis. 

The present paper considers the influence of compressibility on a 
laminar radial wall jet. ‘The effects of viscous dissipation, wall temperature 
conditions and of the initial temperature of the fluid in the wall jet are 
analysed. Some of the results of this study have been briefly referred to 
by Glauert (1957). 

Glauert (1956) proved that, in an incompressible wall jet, the ‘flux of 
exterior momentum flux’ is constant; this forms the basic starting point 
of his theory. It is here shown how this result can be generalized to 
compressible flow, the flux of exterior momentum flux being constant 
when the viscosity is proportional to the temperature. 

In the boundary-layer equations of momentum, continuity and energy, 
the stream function may be conveniently chosen as independent variable 
in place of the coordinate normal to the wall. This transformation was 
first used for incompressible flow by von Mises (1927) and later by 
von Karman & Tsien (1938) for compressible flow. Under the simplifying 
assumption that the viscosity is proportional to the temperature, it is possible 
to solve the momentum and energy equations independently, the momentum 
equation becoming identical with the corresponding equation for 
incompressible flow, and the same solution being applicable. The 
interpretation of the results in terms of geometrical coordinates differs, 
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and is fully investigated. In particular it is shown that the skin friction 
has the same value as in the corresponding incompressible flow. 

Still retaining the simple viscosity-temperature law, but imposing no 
restrictions on the Prandtl number, solutions of the energy equation are 
obtained. These solutions describe the effects of viscous heating, wall 
temperature and initial jet temperature on the velocity and temperature 
profiles, both when there is no heat transfer across the wall and when the 
wall is maintained at a constant temperature. If the assumption that the 
Prandtl number is unity is made, the well-known Crocco relation between 
the temperature and the velocity provides a particular solution of the energy 
equation, being a special form of the solutions obtained in the more general 
analysis. When a more accurate viscosity-temperature law is assumed it 
is possible to develop solutions of the momentum and energy equations 
in series form, though the details are not set out here. 

In practice, wall jets will almost always be turbulent. No attempt is 
made in this paper to extend Glauert’s analysis for turbulent flow, since 
it is not easy to predict how the eddy viscosity and diffusivity will vary with 
temperature. However, the general nature of the changes in the velocity 
profiles in the laminar case may perhaps serve as an indication of the sort 
of modifications to be expected for turbulent compressible flow. 

Glauert (1956) shows that there is a close analogy between wall jets in 
two and three dimensions, the same velocity profile occurring in each case. 
Similarly all the results obtained here, with simple changes, become 
applicable to compressible plane wall jets. 


EQUATIONS OF MOTION 
On the boundary-layer approximation, the momentum, continuity and 
energy equations governing a compressible, laminar, radial wall jet flowing 
over a plane wall are 


e 0 Of a1 
pits t pv = = ( =) (2.1) 
ex dy dy\" dy 
e 0 
— (pxu) + — (pxv) = 0, (2.2) 
Cx CV 


wal yp OF _ @/ aT\. up bs 2 
iS it —_ = = Cee 8 sees: Fe . 
Ox 7 ey o oy oy & =) 


The boundary conditions areu = v = Oaty =0; u>0, T>T, asy> a, 
together with a condition on T at y = 0, x = 0. 

Here x and y denote distances along and normal to the wall, x being 
measured from the jet axis, w and v the corresponding velocity components, 
p the density, 7 the temperature and y the coefficient of viscosity. The 
specific heat at constant pressure C), and the Prandtl number o are assumed 
constant. ‘The subscript «x is used to denote values at v= «x. On the 
boundary-layer approximation the pressure is uniform everywhere and 


hence the equation of state implies that 


pT = const. (2.4) 
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For an incompressible wall jet, Glauert established an integral relation 
which he interpreted as saying that the flux of exterior momentum flux is 
constant. We can follow the same procedure as Glauert to obtain a 
corresponding result for compressible flow. Multiply equation (2.1) by x 
and integrate with respect to y between the limits y and ©, using the 
condition that u > 0 as y > «©; then, since 

co Cu oe .- 
| pxv dy = [pxvu]? — = (pxv) dy 
Jy cy Jy coy 
-2 fa] 
= —pxcu+ | u— (pxu) dy 


“yy 


using (2.2), we have 


wn 
~_— 


Z4 (pxu?) dy — pxt ree. = @, (2. 
Jy 7 oy 


Multiplying (2.5) by pxu and integrating with respect to y between the 
limits 0 and «, we have 


4 ee “x ) “cm 7 ( ~@ ) 
| pxus | pxu? dy>dy— | — (pxu) | pxu? dy> dy — 
-_ Sp | Os by 


\ 


¢ 





-e 


Ley 
:< ngs ou 
— | p?x®ou® dy + | pux*u—dy=0. (2.6) 
Jo : Jo oy ~ 
From the continuity equation, the second term of (2.6) is 


*2 2 - co 
pxt | pxu® dy + | p?x?vu? dy. 
- 0 “ O 
Now at vy = 0, u =v = 0, therefore equation (2.6) reduces to 
a , . . A 
g (° wi i 9. Ou 
=~ | pxu< | pxu?dy>dy+ | ppx*u— dy =0. (2.7) 
OX Ja Jy F i /0 oy : 
For incompressible flow, where p and pu are constants, the second term 
of (2.7) is zero. For compressible flow it will vanish only in certain 
circumstances. In particular it will do so if pu is a function of the velocity u, 
or is constant. In these cases (2.7) reduces to 
pxui | pxu® dy) dy = F, (2.8) 
y J 
where the constant F is the flux of exterior momentum flux. It may be 
noted, that as defined here, the constant F differs by a factor p? from the 
definition used by Glauert (1956). 
The continuity equation (2.2) implies the existence of a stream function & 
such that 


Ous 
pxu = ae ’ | 
‘ | 2. 
a (2.9) 
Ou | 
xv = —-—. 
f Ox 


Following von Mises, we now take x and & as the independent variables, 


F.M. 2R 
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and hence write 


3 3° nH) 
7 =a TPC Sos 
(=), ~ (a), (a). 
7] y | (2.10) 
(5). - (3) | 
— = pux —_— e ] 
Oy) z Cus}, J 


Substituting (2.10) into the momentum and energy equations (2.1) and 


Cu Pas, Ou 
= x2 puu— ; (2.11) 
Cx C ub ’ ¢ us 


(2.3) we have 





of <x a oT ppx*u ou\? 719 
— =— =| peu — )+ - —). (2.12) 
Ox Co x us r4 us ( p C yb 
On introducing dimensionless variables by writing 
u= Ud, x = v,, &/U, b = px v2, £/U, | (2.13) 
it : es \ 243 
P = PxPs B= BoB, T=T,f, J 


where U is an arbitrary constant velocity and v the kinematic viscosity, 
we find that equations (2.11) and (2.12) become 





Since pT is constant by (2.4) ,we have, if uw « T, pu = const. = px, or 
pa = 1. (2.16) 


Using (2.16) we see that equations (2.14) and (2.15) become 








3. VELOCITY DISTRIBUTION 

The density f, and hence the temperature 7, does not appear explicitly 

in the momentum equation (2.17), and so (2.17) may be solved independently 

of the energy equation (2.18). Moreover a particular solution of equation 

(2.17) must be the similarity solution, obtained by Glauert (1956), of the 
corresponding incompressible equation, that is 

i= 3 f’(n)x-32, (3.1 

where f(y) = bz" (3.2) 


and f(7) satisfies the equation 


"+ ff" +22 =0,) 


f(0)=f'0)=0, f(x) =1, 
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The boundary condition at the outer edge of the wall jet only requires 
f'(«) =0, but writing f(0) = 1 involves no loss of generality, since the 
velocity U in (2.13) remains arbitrary. Glauert obtained the solution of 


5 


equation (3.3) in the form 





2 Un? as @ tein nt Ve 2 
fxs’, n = log ae + V3 tan7! (+e) (3.4) 
More recently, Glauert (1957) has shown that 
a= if (n)(8 +P), (3.5) 
where f(y) = b(#8 +B), (3.6) 


is also a solution of the momentum iad (2.17), 7 being an arbitrary 
constant length. Since the velocity given by equation (3.5) remains finite 
as ¥ -> 0, this solution may have some significance even near the axis. We 
shall make use of this latter solution of the momentum equation. 
Although the velocity function f’(7) is the same in the incompressible 
and compressible problems, the interpretation of the variable » differs. 
Connection of 7 with the physical variable y can be made as follows. Using 
(2.9) and (2.13) 
pe | a ee df. (3.7) 
: lo pxu Us sJ 9 pux 
Substituting p../p = T and hemodiiahis a and » according to equations (3.5) 
and (3.6) we obtain, from equation (3.7) 
Tf ak 3\ -” 
a (SS) 3 |" Td. (3.8) 
Vi X J0 
Thus, in order to obtain the velocity in terms of the physical variable y, 
we see that the temperature distribution throughout the wall jet must be 
known. The temperature distribution will be the subject of §4. The 
skin friction is independent of the temperature distribution when p « T, 
since it is given by 


‘ 5), 07 ahr oc) (3.9) 


which has the same value as for a oie flow. 

The arbitrary velocity U may be expressed in terms of the flux of exterior 
momentum flux F defined by equation (2.8). Substituting the values 
of (3.5) and (3.8), we obtain 

_. & whee 
F=———. 3.10 
40 UL aaah 
Proceeding exactly as in the corresponding incompressible case, if the wall 
jet is the result of a free jet impinging on a flat plate, a rough estimate of 
the magnitude of F may be obtained from the conditions in the free jet as 
F = } x (typical velocity) x (mass flow/radian)?. (3.11) 


The constant length / can be estimated by equating the maximum velocity 
at x = 0, given by (3.5), to the maximum velocity in the impinging jet. 
2R2 
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Alternatively, if the wall jet is formed by fluid flowing out radially, 
with uniform velocity U, and temperature 7), from beneath a circular disc 
of radius a and at a distance h from the wall, then equation (2.8) may be 
evaluated at x = a to give 
U? h?a? ‘ 
a (3.12) 

216 
In this case / can be estimated by requiring the maximum velocity at x = a, 
given by (3.5), to be Up. 


F= 


4, "TEMPERATURE DISTRIBUTION 
In order to determine the temperature distribution, the energy equation 
(2.3) must be solved for 7, When the Prandtl number a is unity, a particular 
se of (2.3) is the Crocco relation 


at 
T- 2G, A+ Bu, (4.1) 
where A and B are constants. —_ u->( as y-> 0, A must have the 
value 7,,. The value of B depends on the degree of heating in the initial 
jet. In particular the solution with B = 0 will be appropriate for a jet 
issuing from a reservoir of fluid maintained at the same temperature, T,,, 
as the fluid outside the wall jet. In this case the heat transfer across the 
wall is seen to be zero. If we suppose the wall jet to be formed by fluid 
flowing out uniformly from beneath a circular disc as in § 3, then equation 
(4.1) can be evaluated at x = a, giving 
+ 3(y—1)MR-1 

Uy 
where y is the ratio of the specific heats and My = U,/a,, a, being the 
velocity of sound outside the wall jet. 

We shall now seek solutions of the energy equation for arbitrary Prandtl 
number. The simple solutions, for Prandtl number unity, given by (4.1) 
occur as special cases of these more general solutions. 

When the assumption » «< T is made, the energy equation takes the 
form given in equation (2.18), a linear equation for T. Since a(#,7) is 
known, it is convenient to change the independent variables again from 
(%, 75) to (%,7), where 7 is defined by equation (3.6). To do this we have 
the transformation formulae 





B on (4.2) 























(si), 7 “a7 Ba). 
dx/5  \Ox 4 f’ (®+B)\dn] =’ | 
; alla (4.3 
O\ _ 7] 
(a5); 77 of (B+P ii) 
Substituting for @ frcini equation (3.5) the final form of the energy equation 
(2.3) i 1S 
02 T 4(#8+3),0T oT ¢ j2 "2 
at (5 ze pe 7). 9 UL j 4.4) 
07 > OX ~*~ On "6 C726 (+P) 
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This equation (4.4) is the form in which we shall study the energy equation. 
A boundary condition at the wall, 7 = 0, must be specified. We shall 
consider only the two cases of a perfectly conducting wall, over which the 
temperature is maintained at a constant value, and of a thermally insulated 
wall, at which the heat transfer is zero. ‘The outer boundary condition 
is T=1. A condition at * =0 must also be specified. This is related 
to the temperature of the initial jet. When the solution for T has been 
found, the correspondence between 7 and the physical variable y may be 
obtained from equation (3.8). 
A particular integral of equation (4.4) may be found in the form 


T = 14+C,(#+P)-10,(n). (4.5) 
Equation (4.4) also has an infinity of complementary functions of the form 
T,, = C,(#+B)-*6,(n), (4.6) 

where 8, satisfies the equation 
0” + o( f0), +4f’6,) = 0. (4.7) 


By adding to the particular integral (4.5) suitable complementary functions 
of the type (4.6), the temperature conditions we wish to be able to prescribe 
can be satisfied. 

The effects of viscous heating, wall temperature and initial heating 
will be discussed separately, and the solutions appropriate for both the 
required conditions at 7 = 0 will be determined. 


Viscous heating 

Let us first examine the temperature distribution due to viscous heating 
in the wall jet. As we have seen, we require a solution of equation (4.4) in 
the form (4.5), where, if we choose 


9 U? 
“ae (4.8) 
the equation to be satisfied by 4)(7) is 
6” +0( f0, +4f'0,) = of”, (4.9) 
where 
6,(0) =0 } 
and either 6,(0) =0 for constant wall temperatureT,,, | (4.10) 
or 6°(0) =0 for a thermally insulated wall. | 


When o = 1, 6(7) = }f'7(m) satisfies all the conditions (4.10), a result 
indicated by (4.1). For arbitrary o, no solution of equation (4.9) has been 
found in closed terms. The equation has been integrated numerically for 
o = 0-72, the value appropriate for air. 

When the wall jet is flowing over a thermally insulated wall, the wall 
temperature is, from equations (4.5) and (4.8) 

9 U?  4,(0) 

16 C, Ta (@ +P) 
The numerical integration gave 4,(0) = 0-0060. 


at (4.11) 





uw 
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When the wall is maintained at a constant temperature 7\,, the rate at 
which heat is transferred across the wall, per unit area, is 


er as) (4.12) 


dy 


ZT \ pl3 yt \14 1 
eset ees Se A’ (0), 43 
64+ — Ge ~ PYG) M9) s ‘. 


using equations (3.8), (4.5) and (4.8), where A is the thermal conductivity. 
The numerical integration gave 0/(0) = 0-0036. 

If (1 —T,,)/T,, > (Mach number)? at all points in the wall jet, then 
the effects of viscous heating may be neglected and the particular integral 
(4.5) may be taken as T = 1. 


Wall heating 

To describe the effects on the temperature distribution of maintaining 
the wall at a constant temperature 7, 4 T,,, we require a complementary 
function (4.6) with « = 0. Thus 


T, = (T,,—1)6,(n), (4.14) 

where 
¢,(0) = 1, ” 
@,( ne (4.15) 


The required solution of equation (4.7) is immediately obtained as 


- { “7 1 { +x -” 7 
6.(7) = | exp: a| f dy > dy | eXp- G | f dy > dy. (4.16) 
As shown by Glauert (1956), f = g?, where g’ = }(1—g°), and hence 


— | f dy = log(1—g’*), (4.17) 


“0 
and equation (4.16) may be written 


-] / +] 
0,= | (1-g*)?"dg/} (1—g*)’ dg. (4.18) 


Jg [| Jo 
When o = 1 this reduces to 
6, = (1-g). (4.19) 
When o # 1, 6, may be represented as an incomplete beta function by 
writing (1 —g*) = ¢ in equation (4.18), which then becomes 

A, | t? = ¢ 7 t) 3/3 dt | | 7° 1(] —t) 23 dt 


/ 0 / / 0 
= Ba, 4)/B(o, 4). (4.20) 
The contribution to the rate at which heat is transferred across the wall, 
per unit area, is given by 


, ™ a i \l4 
tem —DPe— TI TeseRR) WO) —_ 


from equations (3.8), (4.12) and (4.14). Expressing the beta function 
in terms of [’ functions, we find that (0) = —0-2861. 
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Initial heating 
When the jet is heated initially, we require another complementary 
function of the form (4.6), 


T, = C,(%5 + P)-*6,(7) (4.22) 

where 6@,(7) satisfies equation (4.7). The boundary conditions are 
6,(00) = 0 (4.23) 
and either 6,(0) = 0 for constant wall temperature, (4.24) 
or 65(0) = 0 for a thermally insulated wall. (4.25) 


The determination of x is an eigenvalue problem. We shall consider 
first the case of a thermally insulated wall. Integrating equation (4.7) 
with respect to 7 between the limits 0 and y, and using (3.3) and (4.25) 
we have 


-” 
4 


6, +040, f+(4x—-1) | f'0,dy> = 0. (4.26) 
L /0 J 
On letting 7 — ~, equation (4.26) becomes, from (4.23), 
(42-1) | f'0, dy =0. (4.27) 
/0 


Unless the integral in (4.27) is zero, which on physical grounds is unlikely, 
this shows that x=}. With this value of x, equation (4.7) may be 
integrated in closed terms, using the result (4.17), to give 


6, = (1-g)’. (4.28) 
It follows that the contribution to the wall temperature is 
To, = Co(#? + )-"4. (4.29) 


The value of the constant C,, which depends on the initial conditions in the 
jet, is discussed later. 

For a wall maintained at constant temperature, it may be verified that 
the solution of (4.7) satisfying the boundary conditions (4.23) and (4.24) 


is given by 


Pa) 9 (4.30) 
dso 
6, = g(1—g?)’. (4.31) 


When o = 1, these equations become x = }, 6, = 3f’, in accord with the 
term Bu in equation (4.1). The contribution to the rate at which heat 
is transferred across the wall, per unit area, is given by 


re p,,\ (119 +3) 20 x4 ae 


from equations (3.8), (4.12) and (4.22), where, from (4.31), 65(0) = 4. 
An estimate of the constant C, may be obtained as follows. The 
2 


when o = 0:72. Hence, by requiring the maximum temperature at * = 
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(or ¥ = a when the jet flows from beneath a circular disc) to be the excess 
of total temperature at this point in the physical flow, the value of C, may 
be roughly determined. 

Since the equation (4.4) is linear, we may construct a solution which 
is any linear combination of the solutions obtained above, and thus we are 
able to deal with a wall jet of arbitrary Mach number, wall temperature 
and initial heating. When the wall is maintained at constant temperature 
we require (4.5), (4.14) and (4.22), the total rate at which heat is transferred 
across the wall, per unit area, being given by a combination of (4.13), 
(4.21) and (4.32). For a thermally insulated wall (4.5) and (4.22) are the 
appropriate solutions, the wall temperature being given by a combination 
of (4.11) and (4.29). 


5. RESULTS 
The results obtained in §3 and § 4 are presented in figures 1-5. Figure 1 
shows the velocity function f'(7) whilst figures 2-5 show the temperature 
distribution across the jet and the resulting dependence of the geometrical 
distance y on 7 for all the cases considered above. This dependence of y 
on 7 is given directly by equation (3.8) when the temperature distribution 
is known. Since the velocity u is a function of the mass flux ys, a decrease 
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Figure 1. The velocity function f’(7). 


of p (corresponding to an increase of 7) causes a physical broadening of 
the wall jet, and an increase of p causes a narrowing. The broadening or 
narrowing occurs most where |7—T,,| is greatest. Since, in general, the 
temperature distribution across the wall jet changes as the jet develops, 
it follows that even for our similarity solution the geometrical form of the 
velocity profile is not independent of x, as it is in the corresponding 
incompressible flow. 
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The effect of viscous heating is shown in figure 2 where the temperature 
distributions across the wall jet, for a thermally insulated wall and for a 
wall maintained at constant temperature 7.., are given from equation (4.5). 
The simple Crocco relation for o = 1, given by equation (4.1) with B = 0, 
is included for comparison. In this case the total energy in the wall jet 
is constant, so the effect of a decrease in velocity is exactly balanced by an 
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Figure 2. T'he temperature distribution due to viscous heating; 
thermally insulated wall, 
— — — — wall maintained at constant temperature T,,, 
a Crocco solution for o = 1. 


increase in temperature. When the wall is maintained at a temperature T,, 
there is no heat transfer across the wall and hence both conditions at the 
wall are satisfied together. For o=0-72 this energy balance is not 
maintained, and a solution for each of the wall conditions is required; 
qualitatively, however, the results are similar to the case o = 1. Figure 3 
shows how the wall jet is narrowed in the region of higher velocity, in 
comparison with the profile at low Mach number. The Mach number M 
is here defined as the ratio of the maximum velocity across the wall jet 
to the velocity of sound outside the wall jet. 
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When the wall is maintained at a constant temperature T,,, 4 T,, the 
contributions to the temperature distribution and to the variation of y 
with 7 are as shown in figure 4+. These results were obtained by numerically 
integrating equation (4.20) for the value o = 0-72. As we might expect, 
the temperature increment is greatest near the wall, gradually falling to 
zero with increasing distance from the wall. 
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Figure 3. The variation of » with 7 due to viscous heating; 
——————————. thermally insulated wall, 
wall maintained at constant temperature T,,, 
—+—-+—-+— Crocco solution for o = 1. 


Figure 5 depicts the contribution to the temperature distribution and 
to the variation of vy with » when heat is added initially to the jet, for the two 
cases of a thermally insulated wall and a perfectly conducting wall. Fora 
thermally insulated wall the effects of the added heat are similar to those 
shown in figure 4, though the variation with distance from the axis is 
ditferent. For a perfectly conducting wall, the effects of the added heat 
are mainly confined to the faster moving regions of the wall jet. 


As we have seen in §4, when a hot wall jet flows at high speeds over a 
heated wall, the temperature distribution can be obtained by taking a suitable 
linear combination of the distributions that have been considered above. 
Similarly, we see from equation (3.8) that to obtain the total variation of y 
with 7 we need the same combination of the appropriate distributions shown 
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in figures 3—5, the total increment in y being the sum of the increments of y 
in the separate distributions. As the distance from the axis of the jet 
increases, the Mach number decreases and the added heat is diffused away 
until their effects on the temperature eventually become negligible. We 
see from equation (4.5) that the effects of viscous heating, which fall off like 
(x? + /8)-1, decrease faster than the effects of added heat which, from (4.22), 
behave like (x?+/*)~*, where « = 0-25 for a thermally insulated wall, and 
2 = 0-548 for a perfectly conducting wall when o = 0-72. However, the 
effect of a wall temperature 7, 4 T,, does not decay with increasing x. 
As a result the form of the velocity profile changes as the wall jet develops, 
and only when the effects of viscous heating and added heat have become 
negligible is geometrical similarity achieved, the final velocity profile being 
the incompressible one unless the wall is maintained at a temperature 
different from T.,. 


The author is indebted to Mr M. B. Glauert for suggesting this problem 
and for his encouragement at all times, also to the Department of Scientific 
and Industrial Research for a maintenance grant whilst this work was being 
carried out. 
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SUMMARY 

In this paper transonic small-disturbance theory is applied to 
a simplified model of the flow near the front of a ducted body. 
The body is assumed to consist simply of two parallel flat plates 
which extend from the inlet station to infinity downstream. The 
velocity far upstream is sonic, and the velocity far downstream 
in the duct, which is assumed to be known, is slightly subsonic. 
Air is therefore ‘spilled’ around the intake edges. An analytic 
solution is found for the resulting flow field up to the ‘limiting 
Mach wave ’, and asymptotic solutions are found for the supersonic 
flow and for the shock wave far from, and near, the intake edges. 
The pressure distribution along the outside walls is then known 
at both ends, and its computation is completed by an empirical 
procedure. Distributions of pressure along the centre-line and 
along the inside and outside walls are shown. ‘These results may 
be used to compute the drag of sharp-edged intakes with a very 
small frontal area. 


1, INTRODUCTION 

Under most flight conditions, the intakes of high-speed aircraft ‘ spill’ 
air: that is, the mass-flow requirement of the engine is such that the entering 
stream-tube of air has a cross-sectional area, far upstream, which is smaller 
than the area of the hole at the front of the body. As a result, air is diverted 
around the intake edge, and at high subsonic and at all supersonic flight 
speeds the flow field is of mixed type (see figure 1). In the subsonic case 
there are regions of supersonic flow near the edge and, possibly, further 
back along the external wall of the intake. In the supersonic case there is 
a subsonic region which is bounded by the bow shock wave, by a sonic 
surface from this shock to the intake edge, and by the walls of the duct. 
As far as the external flow is concerned, the outer wall of the intake acts 
somewhat like the upper surface of an aerofoil, increasing spillage corre- 
sponding roughly to increasing incidence. ‘These mixed intake flows have 
been the subject of considerable experimental work, and various empirical 
theories have been developed to predict the effect of spillage on the drag 
of aircraft, but (as far as the writer is aware) nothing approaching a 
theoretical description of the flow field has been developed. 

In the present paper we study the simplest possible example of such 
a flow. The intake is assumed to be two-dimensional and to consist simply 
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of two parallel, semi-infinite, flat plates. Alternatively, the upper half of 
this configuration may be interpreted as a simplified model of one of the 
side intakes found on some aircraft. The final velocity in the duct (which 
in practice is known from the air mass flow of the engine) is assumed to 
be only slightly less than that in the free stream, so that perturbations from 
a uniform stream are small, except in the neighbourhood of the intake edge. 
If the free-stream Mach number is near one, transonic small-disturbance 
theory may then be applied, and the central problem, in which the free 
stream is sonic, is considered here. 








\\ / \ 
\/ 
0 SUBSONIC \ \ 
\ \ 
—>— = STREAMLINE (b) SUPERSONIC 
—-—-- SONIC LINE 
—— SHOCK WAVE 
Figure 1. 


It is usual in intake work to use a definition of drag which for plane flow 
becomes ; 
D =2 | (P—P,) ay, 
where P, denotes the pressure in the free stream, and the integral is taken 
over the entire length of the dividing streamline, ASEB in figure 1. This 
definition is consistent with the usual definition of engine thrust. It is 
natural to consider the drag as the sum of two parts: (i) the ‘pre-entry 
drag’, which is the contribution of the upstream part, 4S, of the dividing 
streamline, and can be calculated from the known change in momentum 
and pressure of the internal flow; and (ii) the ‘cowl drag’, which in the 
case of a flat-plate intake is simply the force on the edge E. In incompressible 
flow theory these two drags precisely cancel each other, as might be expected, 
regardless of whether exact or small-disturbance theory is used. In the 
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present transonic problem, however, the singularity in velocity at the edge 
is rather weak, and the edge force is found to be zero. Both in incompressible 
and in transonic flow, the edge flow has, of course, the same character as 
that at the leading edge of a lifting flat plate, for which the transonic solution 
has been found by Guderley (1954). Hence the result of zero edge force, 
which was not discussed by Guderley, applies to that problem also. 

Some further information about cowl drag can also be obtained from 
the present model. Consider a wedge intake (that is, a two-dimensional 
intake of the type shown in figure 1) in a sonic stream. For moderate wedge 
angles the flow up to the ‘limiting Mach wave’ (figure 2 (a)) is the same 
as that past a flat-plate intake. (The significance of the limiting Mach 
wave will become clear from what follows.) Further, if the wedge angle 
is of even smaller order than the slope of the dividing streamline ahead 
of the intake, then the dominant part of the pressure acting on the wedge 
face is that on the outer wall of the flat-plate intake. Hence for such 
slender-wedge intakes a first approximation to the (negative) cowl drag can 
be found from the present solution. 


2. 'THE EQUATIONS OF TRANSONIC FLOW 
Let X, Y be Cartesian coordinates such that a free stream of near-sonic 
velocity U,) and of Mach number M, flows in the direction of increasing X. 
Let the velocity at a point in the field be (U)+U, U’, U,V’). The fluid 
is assumed to be a perfect gas whose specific-heat ratio is y. Then for 
small disturbances the equations of transonic flow may be written 


(1—M2—TU’)U;,+Vj;, =0, or 


2.1 
a oh | sot 


} Xx 
where subscripts denote partial derivatives, and where we have adopted 
Spreiter’s (1954) version of these equations. In order to work with 
variables whose magnitude is O(1), and to express the equations in 
canonical form, we introduce a small parameter 6, such that V’ ~ O(8), 
and a length L, representative of the lateral extent of the field, and make 
the transformation 


X = L(sl)'3x, Y = Ly, ] 
(2.2) 
U’ = (1— M2) 1+ 642-18 u(x, y), V’ = 5v(x,y). 
Then equations (2.1) become 
—uu,+v, = 0, | 
' (2.3) 
u, = G. J 


with the boundary condition 
u->(M2-1)(8T)-2° = K, 050, as%y> oo. 


Here K is the transonic similarity parameter. 








632 L. E. Fraenkel 


If the roles of the independent and dependent variables are now 
interchanged, equations (2.3) become 
x,—uUy, = 0, 
x,—y,, = 0, 
so that the Legendre contact potential, y = { (x du+y dv), satisfies the 
Tricomi equation 
Xuu— "Xo = 9. (2.4) 
The characteristic curves of this equation are v = 3u°?+ constant, (u > 0): 
hence we introduce the variables ¢ = 3u? in the hyperbolic half-plane, 


and w= —u, zx = 3w>? in the elliptic half-plane. Then the equation 
becomes 

1 2. 7 

Xeor 30 Xe— Xov = 0, XzzT 3, X27 Xvv = 0. (2.5) 


Two types of elementary solution of this equation, and the relations between 
them, are discussed in Appendix I. 

To the lowest order, the stream function is ‘Y’ = [py U,) dY, so that 
y (= x,) may be taken constant on streamlines. ‘The variable y also satisfies 
the Tricomi equation. 


3. SOLUTION OF THE INTAKE PROBLEM UP TO THE LIMITING MACH WAVE 
Figure 2 shows qualitatively the physical and hodograph planes of the 
present problem. The origin of coordinates is taken at the intake edge, 
and the length Z in (2.2) is taken as half the distance between the two 
flat plates. The final velocity in the duct being U,, we define 6 by 
s23p-3 — Uy Uy (3.1) 


0 
By continuity, the overall change in width of the entering streamtube is 
then O(6*3), and this occurs essentially over an X-length of O(6!*), so 
that the slope of the dividing streamline is O(6). 

The physical and hodograph fields may be described as follows. ‘The 
origin in the physical plane maps to infinity in the hodograph plane, and 
vice versa. ‘The most important streamline is the dividing streamline, 
y=0. As this approaches the stagnation point, S, it has positive slope 
(v > 0) and the velocity is reduced (w>0). At S (x=0,y=0,w= ~, 
v = 0 on the approximate theory) the dividing streamline branches; one 
part, vy = 0—, running along the inner wall SJ, while the other, y = 0+, 
expands at infinite velocity about the edge FE. In fact it overexpands, 
as do all the external streamlines, and ultimately attains negative slope 
(v <0). Recompression and the return to zero slope occur first through 
a shock, and then continuously. The shock corresponds te a jump in 
position in the hodograph plane, along the appropriate shock polar, and the 
subsequent flow maps on to a second sheet of the hodograph, which is here 
shown superposed on the first. 
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sonic line ——>} ri 


limiting Mach wave ~,! 
































O shock ! 
ies Yo 
y= 0 
stnitiarinanea acta y =0- 
yey, 
0 aki I 
‘ yee 
(a) PHYSICAL PLANE 
V 
/ 
E / 
Dae O* Py 
(at infinity) E = _s ‘imiting Mach wave 
a 
w < u 
(b) HODOGRAPH PLANE 
Figure 2. 
The central streamline of the intake, y= —1, proceeds from O to I 
at zero slope. The internal flow is contained between y=0, y =0-, 
and y = —1; the external flow between y = 0, y = 0+, and the shock. 


The limiting Mach wave is the last of the expansion family emanating from 
the edge to meet the sonic line; hence disturbances downstream of this 
line do not affect the subsonic field. It maps on to v = ¢. 


F.M. 25S 
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In this section we consider only the flow up to the limiting Mach wave, 
that is, the problem in A(v—¢) > 0. This is essentially a boundary value 
problem of the Tricomi type. ‘The solution is constructed in two parts 
as follows (y = x; + x2). 

(1) The potential y, is to represent the far field of the physical plane, 
that is, the flow near O. We require: 

V1 = Xiy > © as u,v, 

¥,70 as u,v> o, 

Vy =) On 2 = 0, wW> 0, 

y, finite onv =f, v > 0, |} 
and a one-to-one mapping from the hodograph to the physical plane. 
This ‘free stream singularity’ for any symmetrical body has been found by 
Frankl (1947) and Guderley (1948). It is one of the similarity solutions 
x = ¢-"fn(x), where « = v?/{? (Appendix I). In the present case m = 3, 
and the fn(x), which is a hypergeometric function, reduces to elementary 


form. In fact 


* 





) Gia sho ChS | e-4°K 4 19(AC) aA, (3.3. 
‘osh!! cos 1 > Se 
— pea = =a rN ) (v20, «a< 0), (3.3b) 
Vv (1-2) 
Ss] Sl oy l eee ) 
= —p2¢ 23 inhi sinh™ + (2 a)5 (v>0, #21), (3.3c) 
\/(x—1) 


where yz 1s a constant to be determined in any particular problem. Although 
it is irrelevant at this stage, we also note that 
y--93 Sint} sin~! +/(1— 2), 


X1 Meg a 


‘ << i) 53 
v (1-«) iat. saath 





where sin-!4/(1—«) varies from 0 to 7 as v varies from ¢ to —€. 

To formulate precisely the problem for y., we observe that y, gives 
velocities near the edge FE which are altogether too large, and must be 
removed by y,. Indeed if y, were to be the dominant term of the complete 
solution near F, this would result in an infinite drag, as will be shown in § 6. 

(ii) ‘he potential y, is introduced to satisfy the boundary conditions 
on the body: 

Ve —lonv=0, 0O<w<l, 
=0 on v=0 w> i, 
VYoty, ~ 0(y,) as w,v> 0, 


¥y is finite on v = €. 


Although p is as yet unknown, the second of these conditions does fix the 
order of magnitude of y, at infinity in the hodograph plane. It is shown 
in Appendix II that (3.4) determine a unique boundary value problem in 
the sense of Tricomi. 
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We seek a solution of form 


/ - 00 
xe= — Ea [eK QQ /0)>, 
where the J,/5-solution has been omitted because it leads, in general, to 
a singularity on the limiting Mach wave. (‘This step is justified rigorously 
by the fact that (3.5) will satisfy all the boundary conditions, and is therefore 
the unique solution.) f(A) is then determined from the integral equation 


y,(2,0) =3"8 | [Jyg(Az) +7 aa(As)]{(A) dA = -1, 3 >2>0, 


v0 


by means of the Mellin transform 
+00 


B(p) = | 2°14 (e) de 


(see Appendix I for further details). ‘The integral equation transforms to 
Yis(p+ 3) +I as(P + 3)1f(—p + 3) = — (3)? /p, 
where (Erdélyi 1954) 


> sin }zpcos ba Zp 


J 3(p) +9 _1(P) a cae tune ye tp 7 ae p+ ty’ (3 < Rp << 3). 





Hence, with —p+2 = q, we find after a little reduction that 
27-3 cosmq+cos jz ,, : 24 
f(q) = (3)?8-9@ —— 2 dg a (49+ 3), (\Aal <1), (3.7a) 


La Eh Ss 
7 COS 47qQCOS ¢7 


) 


and 
fit) =| r-*F(q) dq, (le < 1). (3.7b) 
471 J ein 
The only singularities of f(q) are simple poles at 
g= +(2n+1), w= OE A cs 


the poles of the ['-functions being removed by zeros of cos7q+ cos 47 
The evaluation of the contour integral in (3.7 b) is therefore straightforward : 
we close the contour with a large semi-circle in Zq < 0, whose contribution 
vanishes in the limit for all finite A, because of the strong convergence 


induced by the I’-functions in (3.7a). Then 


, 23 4 2 — fR\2n41 1 - 
ss aa (5) v3 2 (-) (5) T(n+%)0(n+4)’ ay 








V n 
2\232,/3,  /2\23 84/3 
\ D) \V 2Q 
— = = A+ — eee) ee eee «( 
(3) 5a (5) i | ll aie 


Now consider yx. If the substitution g= —(4/3/7)f is made in 
equation (I.8), and the resulting integral is evaluated by contour integration, 
X2 is given by various series of hypergeometric functions. However, we 
are primarily interested in x(w,0), which gives the pressure distribution 
along the centre-line and the inner wall of our intake. Now on v= 0 

282 
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the hypergeometric functions take on simple values, or alternatively (3.5), 


(3.6) and (3.7a) may be used directly. 


Xe(—q— 4,0) = - (Ji 3(-—9g) +d vs(—9)]F(9) 


In fact (3.5) may be written 


(-—1< #q< -}), 


so that, substituting from (3.6) and (3.7a) and inverting, we have 


X2(2,0) = — 


} 1 press. yer 1. cos7q+cos 47 
- —— | = tan }7q—_—__—_ x 
2m 2 COS 7g — COS $77 





l(3q- 4) (3q+ 4 
P(3q+e(3q+% 
COS 77 


X(w, 0) = - | w>72-12 tan arg 


PGq— 3)PGq+3) 
D(og+ oP (q+ 8) 

The last integrand has simple poles at 
g = +(2n+1), 





2n+3, 2n+3, 


dq, 


where n = 0, 1, 2, .. 


1 


) dq, 


+cos iz 
: x 


cos 7q—cos iz 


") 


and for w 2 1 we may use large semi-circles in #g < 0, respectively. Then 
< 5 > ’ 


for w > l, 


where 











x,(w, 0) = bd y (6, +C,, w+d, w)zey3”, 
a=. 
where 

, -Pe-pre+d) 1, 1 | 
"  D(n+4)P(n+1) nn 12n? 
Pe ee |) 
"" T(n+3)P(n+$) on 42? ‘ 

Pint 4)P(n+ 9) 1 7 





i = ow 
ie P(n+1)0(n+2) nn 


for n> &, | 


| (3.10a) 
| 


for n> ©, 


| (3.10b) 


The asymptotic forms of the coefficients show that these series (which 
could also be written as generalized hypergeometric functions) converge 


rapidly except near w = 1, and that they diverge at w = 1. 


However, the 


singular part of the function is easily recognized and removed from the 


series. Equation (3.10a) may be written 


mw*x,(w,0) = —log(1—w-*) — 3[(1-—w-)log(1 — 





w 3) ‘ w-3] By 


a 


3 


= 1 
zs ~*~ Se cuen — 3 
2 («, 5 Bech 5) ~ (34%) 
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where it is to be understood that the terms —(1/n) and 3/{4n(m—1)} are 
to be included only for »>1and for n> 2, respectively. Equation 
(3.11) displays explicitly the logarithmic singularity of x (which corresponds 
to an exponential decay of the excess velocity in the duct), and the new series 
converges at w = 1 like ¥ n-%. Equation (3.10b) can be put into a similar 
form. 

It remains to find the constant », which determines the scale of the 
far-field solution y,: to this end we consider the flow at E (w,v— o). 
A suitable series for x, as v-—> 0, is obtained by expanding f(A) in 
ascending powers of its argument: 


Pel” see gas in 
x2 = vgs eK, a(00)| (5 


JO ated 
-(3) att | a (3.12) 





Applying the condition 
Xo+x1~ (x1) for v>o, 
which follows trom the requirement ot finite drag, and using (3.3 a) for y,, 


we obtain 
2\23 24/3 = 
ex (5) tal (3.13) 
3 


The asymptotic value of the complete solution near the edge F is now 
given by the second term of (3.12), and may be written 


2\23 81/3 a /y, 
ee i eee: CA | 3.1 
XB (5) 5.7.11.7 =(*) lala 


This is an elementary function which can be found from equations 
22 72 ¥ 
(3.3 b) to (3.3 d). 
A more revealing (but less rigorous) method of evaluating p is the 
following. ‘The condition at the stagnation point S may be written 





=O af wae GBA, o= 6. 
Since w is large, the leading term of an asymptotic series for x,(w, 0) suffices, 
and this is given by the first term of (3.10a). ‘Thus we require 


3\ 2/3 44/3 
Ny TX ~ - 12 ) w* 4 or=O 26 w= 0-75F 4, 


Zz 7 





In other words, as w-> ~ on vw = 0, x should — 0 as rapidly as possible, 
in order that the stagnation point of the exact hodograph be simulated as 
closely as possible. ‘This requirement leads to the same value of as before. 


4. THE FAR FIELD AND THE EDGE FIELD IN THE SUPERSONIC REGION 

The far-field potential y, and the edge-field potential y, are both 
similarity solutions up to the limiting Mach wave. Barish & Guderley 
(1953) have demonstrated that when the initial solution is a similarity one, 
the flow continues analytically to the shock wave, which is itself a similarity 
curve, and that the solution behind the shock, while a different function, 
has the same similarity as that before the shock. However, Barish & 
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Guderley, who studied plane and axi-symmetric flows simultaneously, 
used a non-linear differential equation associated with the flow in the 
physical plane, and had to integrate this equation numerically to obtain 
their results. 

In this paper we prefer to work in the hodograph plane, where our 
asymptotic solutions are elementary functions, both before and after the 
shock. The only numerical work required will be to find the roots of 
certain algebraic equations, which provide the constants in the solutions. 

In place of « we introduce the variable 


— 
* 9 ° | a a 
€ = sin-(1—«)!? = ce ia 


un 


which increases from 0 to 7 as v decreases from ¢ to —¢. The upstream 
solutions are both of form 


x = u-"fn(E), (n = 3m), | 
a nm oe (4.1) 
x =u" TAL}, y=u 32k gg (E), | 
and the solutions downstream of the shock will be written 
v= w-"—"k, f,(n), y =u" 32k, 2,(n), (4.2) 
where 
Ws ee a 
1) _ S = Tee 


f, and g, are determined by the similarity index , and by the boundary 
condition y= 0 on 7 =0: they are, in fact, hypergeometric functions 
which again reduce to elementary form. 

For the far field, by (3.3 d): 























24/3 
n ile oo Sa 
7 
2 sin 3 
f,(é) = S88 _ coség,(€), 
siné 
See On Be | 
sin 4€+2sin 3é | 
g,(é) ae (4.3) 
sin? € 
eri 
2cos 17 
f,() : st 4 sin 21(), 
COS ”) } 
Seema ee 
- ) sin 7) Z2s1n 37) | 
2,,(y 3 ’ | 
COs‘ 7) 
and for the edge field, by (3.14): 
3 
n 4, k, faa caecdht 
so 1 ] uae 
ici 8(5 sin Zé — 14sin 2€ + 35 sin 1) 
Filey s a -°" —cos&9 (&), 
sin’ € dae ites 
10(—2sin #é + 7 sin $€ — 30 sin 4€ + 42 sin 3€) 
¢ — > . 3° state 3 = 35 
£,(§) see ae , (+.4) 
sin’ é 
8(—5 cos fn — 14cos 2n + 35 cos $n) 
Jol) = “—— cos? : a gyn ); 


10(—2 sin yn —7 sin =) + 30 sin 4n + 42 sin $n) 





cos‘ n 
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Let suffices (), and (), denote, for the moment, conditions immediately 
before and after the shock, let u,/u, = o, and let k,/k, = «x. The following 
conditions apply. 


(i) x and y must be continuous across the shock: 











fol) FF Salts) _ (4.5, b) 
Fal&1) K Sa($1) K 


(ii) The jump in (u,v) must be along the transonic shock polar: 


1/2 
u,+U,\'* 
Vg—Ty = (ty -4(45"*) , 








bt. 
2 1+o\12 
: ee ‘ 
5 (—o¥*sin y,— cos.) — (1-0) 5 = 0. (4.6) 
3 2 
(iii) The shock slope must be normal to the vectorial velocity jump: 
dx, dx,  U,-, u,+uy\!? 
dy, dy,  Uy—Uy 2 : 
ie. 











n+l f,(&) _ (est), n+1 filma) _ (=) (4.7a,b) 


n+ £a(§1) 2 n+ 3 $r(72) 20 
These are four, rather than five, independent equations, since (4.5a) and 
(4.5 b) combine to 


and so do (4.7a) and (4.7b). The unknowns are €,, 7, o, and «; and 
the solution of these equations is not difficult. From (4.7) one can 
tabulate or plot €, and 7, as functions of o. ‘The left-hand side of (4.6) 
is then a known function of o, and one seeks its zero. ‘There is only one 
such zero within the physically possible range of the variables. ‘The results 
of this calculation are shown in the following table: here c is the constant 


2) 


in the equation of the shock wave vy = cx@"*9 @n+®), 








| | n o é, Neo K r 

| — ene 

| Farfield | ! 0-345 154-0? 27-2 1:124 0-616 
| Edge field | 4 0-735 126-7? 21-2 1-004 1-430 
| | 





5. "THE PRESSURE DISTRIBUTION ALONG THE OUTER WALL 
The edge field and the far field are now completely determined, but 
the flow in the intermediate region is known only up to the shock, whose 
position is also not fully known. Strictly, the solution should be completed 
by means of the numerical method of characteristics, but in the present case 
the following empirical procedure was though to be sufficient. 
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There are two lines u > 0, v = 0 in the problem: the first lies between 
the limiting Mach wave and the shock, and on the first sheet of the hodograph 
plane; the second is the outer wall of the intake‘and lies on the second sheet 
of the hodograph. The first of these lines would probably be used as a 
starting point of the characteristics computation (cf. Guderley 1954). 
Let the values of x on these lines be denoted by x,(u,0) and x,(u, 0), 
respectively. Both these functions are monotonic decreasing and might 
be expected to be qualitatively similar, since both represent a transition 
from the edge field to the far field. Now the foregoing results show that 


, ) 
*i(t0) _ 2.248[1 + O(u2)] for u—>0, 








x,(u, 0) 

and f (5.1) 
Xp(u, 0) = 2-008[1+O(u-)] for u> o. 
x,,(u, 0) = J 


Accordingly it has been assumed here that 


4 0 A 
*o(u, 0) ~ 2 for all u. ( 
x,,(u, 0) 


x,(u,9) is easily calculated by the method already used for x(w,0). 


ui 


2) 


We have 











») 
* “\ J 2 
Xai (u, 0) = _ ==. 
3 C g ‘ i - er. da ) 
x,0(u,0) = — ¥-— 208 | K,,(AL)fA)— }, 
7 Cu yi A ; 
2 2 ee ” cos 7g+cos iz 
= =--— | y242-L2 sec hare steraice beatae 
4271 } ix ~ “cos 7q—cos 47 
'(4g—2)0(3¢+3 
A34 a) (29 3) ag, eres 
I (sq+ 6)! (39 ‘2 é) 
l 
mses. ~sageiial a * = 3n—2 > 
a. * 2, Gal u) . (u > 1), 
= as y I 4 1 . 44 2 i 3n < 
ae: (6, +3¢,u+d, u®)(—uy>", (u <1), (5.3) 


where the coefhcients are given in (3.10). These are alternating series, 
which converge at uw = 1, but the convergence is improved by introducing 
terms in log(1+.u-*) and log(1+u°), as before. 


6. 'THE DRAG 
In view of the singularity at the edge £, which makes our approximate 
solution locally invalid, the intake drag (as defined in § 1) must be calculated 
from the momentum flux across some suitable curve, such as BGO, in 
figure 3. (In this figure, O,, O., /,, J, represent points which are actually 
at infinity.) The general drag integral of transonic theory is derived in 


Appendix III; in the present case it becomes 


-Cyp =s > = | {2uv dx + (v7 + 3u3—1[u]®)dy}, (6.1) 
O- o- Log l 2 / BGO, : ; ‘ 
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where [u] denotes the velocity jump across the shock wave (downstream 
value minus upstream value), and is a function of y only. The term in [wu] 
accounts for the change in entropy and stagnation pressure across the shock, 
and is to be included only at points on BGO, behind the shock. The 
integral (6.1) is invariant under the changes of the path (provided that 
the termini remain at some point on the outer intake wall and at O,), and 
the path may therefore be deformed to BE,E,/,],0,O,. Then the 
integral around the small contour E, £, gives the edge force, and that 
along J,/, gives the pre-entry drag: the other parts of the path make 
no contribution. 








en a 
rd 
0, | y=0 ____ 6 
| E, i, 
| Pt [l, 
O, 


Figure 3. 


Now near E the flow is represented asymptotically by a similarity 
solution, y = u-"fn(«), and the corresponding drag integral is 


Dy, p, ~ O(u-"*32) ~ O(y2n-B2n+9) (6.2) 


(details of this expression are given in Appendix III). 

It follows that (i) the far-field solution y, must be removed in the 
neighbourhood of the edge (x, v0; u,v-+ «) because it would lead to 
an infinite force, and (ii) the true solution yields zero edge force. For x, is 
the similarity solution with n = 1, so that the corresponding edge-force 
integral would — «© as u-» ©. (Since the geometry behind the limiting 
Mach wave may be varied without affecting the flow upstream, there is 
no general possibility of a cancellation between terms tending to infinity 
in the edge-force integral.) On the other hand the true edge flow is 
represented by the similarity solution with n = 4, so that the true edge- 
force integral +0 as u-> ~. 

The pre-entry drag is easily calculated from (6.1). On J, J, this 


reduces to 


= 243 dy = 3. (6.3) 


It remains to find the cowl drag of the slender-wedge intakes discussed 
in $1. Let the slope of the wedge face be «6, and let its length be /L(61)#%, 
where / is O(1). For «0 the pressure changes induced by the wedge 
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are O(e5), and their contribution to the drag is therefore small compared 
with that of the flat-plate pressures. Hence, by (6.1) 

1 l 

—C ~ —2e | udx, 


62 ~ Pow! 0 


where u refers to values on the outside wall of the flat-plate intake. 
Numerical values are given in the next section. 
7. RESULTS 
Pressure distributions along the centre-line, y = —1, along the inner 


wall y = 0—, and along the outer wall, y = 0+, are shown in figure 4. 


rv 


OUTER WALL y= 0+ 
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Figure 4. 


This picture is supplemented by the following asymptotic approximations 
(which, for the external wall, are based on (5.1) rather than on the 
approximation (5.2)) 


u(x, —1) ~ —0-664( -~x)-'2 for x> — 0, 
u(x,0+) ~ 0:704x-}? for x—> ©, 
u(x,0+) ~ 0°-713x-"5 for x0, 
u(x,0—-) ~ —0°712x%-)° for x0, 


u(x,0—)~ —1-—0-147e-"* for x> &, 
u(x, —1)~ —14+0-147e-" for x> 2 
It is evident that the external decay of disturbances before and after the 
inlet station is slow (although it must be remembered that x is a ‘stretched ’ 


coordinate), and that the edge singularity is very weak. It is also striking 
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that over 90°, of the compression of the entering fluid takes place ahead 
of the inlet station, and that the decay of disturbances within the duct 
is extremely rapid, so that one-dimensional theory provides a very good 
approximation to the flow within the duct. 

As a basis for comparison, we quote some results for the present intake 
in an incompressible stream. In this case the reduced variables are given by 


§=—°_—, X=Lx, Y=Ly, U'=5u(x,y), V' = dv(x,y), 


and the small-disturbance solution is 


‘ 1 1 u—1v 
x+iy=- — +log ———— ], 
: a \u—I1v u—iv+1 





so that 


u(x, —1) ~ 0:318x-1 for x>— &, 
u(x,0+) ~ 0-318x— for x—> 0, 
u(x,0+)~ +0-399x-12 for x 0, 
u(x,0—) ~ —1-—0-368e-" for x+ ow, 


u(x, —1) ~ —1+0-368e-"* for x> 0. 
Within the framework of reduced coordinates, the external disturbances 
of the incompressible flow are therefore more concentrated near the 
inlet station than are those of the transonic flow: the extent to which 
this effect is reversed, when physical coordinates are used, depends upon 
the magnitude of 6. 
The pre-entry drag of the intake in transonic flow is 


b.. ae 
eee re: 


Y 


Cp 


” spy U5 21 3 
Slender-wedge cowls, whose wedge face has slope «é, (« < 1), and length 
IL(51)'3, have the following drag: 








The decrease in drag with increasing frontal area is therefore quite rapid, 
at least initially. 


I am indebted to Professor P. Germain and to Professor J. D. Cole 
for a number of most helpful discussions, and to Mrs A. ‘Tingley for help 
with the numerical computations. ‘This work was completed while I was 
on leave of absence at the Guggenheim Aeronautical Laboratory, California 
Institute of Technology, and was partly sponsored by the U.S. Air Force 
Office of Scientific Research, under Contract AF-18(600)-383. 
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APPENDIX I. ELEMENTARY SOLUTIONS OF THE TRICOMI EQUATION 
By separation of variables one finds the following elementary solutions 
of equation (2.5): 


i 


e= "03K, (AC) = 





. 
e Meh 3L J 3(Az) +d _13(Az)], | 





= 
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‘ 
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Figure 5, 


Alternatively, if one seeks similarity solutions of form xy = (~" fn(x), 


where x = 77/@ = 97? 4u8, these turn out to be, in a neighbourhood of 
the characteristic 7 = ¢ (see figure 5), 
i l l 5 
C"Fl-m, ~=m+=3; m+>+; 1-—a}, 
Zz 2 3 18) 
and oe (OZ) 


— | | 1 1 .. 
(x—1) '6F( 5m +e, 5+ 55 “m+ l s), 


with corresponding hypergeometric functions applying elsewhere in the 
field. ‘The singulari 


= 


+ 
L 


ies of the hypergeometric equation, x =0, ~%, 1, 
correspond, respectively, to the line v7 = 0, to the sonic line, and to the 
important characteristics x = 1. As Guderley (1948) has demonstrated, 
one of these two characteristics is generally the limiting Mach wave, which 
divides both the physical and the hodograph fields into (i) the subsonic 
region, together with an initial part of the supersonic region which influences 


the subsonic flow, and (ii) that part of the supersonic region which does 
not influence the upstream field. 
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Most analytic solutions of problems in transonic flow have been obtained 
by superposition of the solutions (I.1) and (1.2), but the relatively simple 
relation between these forms does not appear to have been emphasized 
in the literature. With a view to the present application, we discuss this 
relation here. Consider 


eenad i aie - y 2% Ge ar a 
x= 09 |" e*TgA)K QS) +A ODS, (1.3) 
where g and h are assumed to be algebraic rather than exponential. Since 
forA>« 





~ \at 
08) ~ (gE) Meae 
1 —Av—) 
(27a)! 2 ’ 
the first part of (I.3) converges for A(v+) > 0 (that is in the regions I, 
II, III, IV of the hodograph plane, figure 5), whereas the second converges 
for A(v—L) > 0 (thatis, onlyinI] and II). Thus we expect the Aj),-solution 
to be associated with only the first solution in (1.2), which is analytic on 
v= ¢. Infact by standard integrals, and by properties of the hypergeometric 
functions (Erdélyi 1954, 1953), one finds that 


e-**Ty (AL) ~ 


g13 ] 0 e-Avym—23K (AL) dA 


“UV 





m2 P T(m+8). : 1 1 5 
= (m) e 3) C-mf (3 m, 5 M+-=—~; M+; 1-2), 
Ve yo J 


2rtU3 (m+ 3) ~ 
(m>0O). (1.4) 
A similar result for the J,,-solution involves both the hypergeometric 
functions in (1.2). 
To obtain a further relation we introduce the Mellin transform. The 
definition integral of this transform is 


B(p) = | g(t) di, (A< Ap <B), (1.5) 
“0 
which we also write 4(p) = ¢(£); the inversion integral is 
1 rc+i@ - 


(-»g(p)dp, (A<c<B); (1.6) 





$(f) = = 


and the composition product is 


/c-ia 


b,(p+a)bo(—p+5) = fa [ Ae*¥-14,(AL)p4(A) aA (1.7) 


~0 
(provided that ¢,(p+a) and ¢,(—p+6) have a common strip of conver- 
gence). Now let 
eal? (A) = G(A, 2). 
Then G(p+4,«) is given by (1.4), with m=p, v = fal2, €=1. Now 
if in (1.3) we set h(A) = 0, and take the Mellin transform of both sides, 
we have, by (1.7), 


X(p, x) = G(p +4, «)8(—p—4), 
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provided that g(q) converges somewhere in #q < —4. Hence 





x(¢,%) = ¢r° [ e- Ky3(Ag)g(A) . 
1 phe ,_ eri ie+ 
Oats eia” = 2? 8D (pt 8) 


fl l 1 5 : 
F(5p. 5P+ =; -4 = 1—2)a(—p— 3) dp, (1.8) 
b ee 3 6 


and we have transformed a superposition of the solutions (I.1) into a super- 
position of the solutions (1.2), the weighting function of the first form 
being replaced by its Mellin transform in the second. The complex integral 
can often be expressed as a series by contour integration. 


APPENDIX II. UNIQUENESS OF yp. 

In this appendix the problem for y, is transformed into a problem of 
the standard Tricomi type by means of a certain inversion which leaves 
the equation invariant (see, for example, Germain & Bader (1952)). The 
(z(¢),v)-plane is mapped on to an (a(y), 8)-plane by the transformation 








Vi(v+1)?4+37} Vi(etlP-@l=re= a E: = = Z =e 
PVE TEL VT 
v+1 l si 3 8 
\ (a - 1)? dt V/ 4 (3 ‘a 2 — (2) \/ (8? a?) \ (8? =) 
y(r,t) = p'3n(p, 7). 
he line v = 0, s real, maps on to the semi-circle p = 7, or 27 + (B— 3)? =}, 
0 <7< 1; the quarter-circler— ©,0 <t < 1,mapsontop-0,0<¢7< 1; 
the characteristic at infinity v+¢€ =k, k-- ©, maps on to B—y=0,y < }; 
and the characteristic v—¢ = 0 maps onto 8+y = 1,y < 3. The problem 
for Voy namely 
1—?? 4 t 
(1) 4 —— Vt = (>, — - v) 0, 
r 3) ? 
(11) 4 —l1 on v=0,4>220 
= () on v=U,2z> ey 


(iii) y ~ O(r-°3) for z real, r—> oo, 


~ O(v*3) for real, v> a, 


then becomes 
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The uniqueness of such a problem, in which data are prescribed on an arc 
in the elliptic half-plane and on an adjacent characteristic in the hyperbolic 
half-plane, was established by Tricomi (1923). 


APPENDIX II]. THE DRAG INTEGRAL OF TRANSONIC THEORY 
In this appendix we first develop the general drag integral of transonic 
theory: this involves no more labour than the derivation of the particular 
form which is of interest here. Let ¢ be the perturbation velocity potential 
of a three-dimensional transonic flow, such that the total velocity at any 
point is U,(1+V¢). The governing differential equation is 


(1—M2-T'$y)byy+byy +bzz = 0. (111.1) 


Consider the drag on a portion A of some solid boundary, and let S be 
a surface in the fluid such that S+A form a closed surface. Then by the 
equations of mass and momentum 


D=({{ Ny(P-P,) 4A, (11.2 a) 


TT] 
| 


-— || [Ny,(P—P)+pU3¢y) +pU24,(N. V4)] dS, (III.2b) 


where N is the normal outward from S+4A, N_,, is its X-component, 
and P and p denote the pressure and density. If P—P, is bounded on A, 
(111.2a) and (I1II.2b) are wholly equivalent, but if there is a singularity 
on A, an approximate solution will not be physically valid in its neighbour- 
hood. We assume, however, that the solution is valid in the field away 
from the singularity: then the drag is defined by (III.2b). 

By expansion of the usual equations of inviscid, homenergic flow, we 
find that 


(Vd)? + M242. + M24,(Vd)2— }(2—y) Mig + 





+1(y+1)Mt Y [dy]? + Of(Vd)"%, (IL.3 a) 


) > 1 9 > oY 2) a 

s = 1— M?4,-—}M?(Vd)?+ }(2-—y) Mid? + Of((V¢d)3}.  (IIL.3 b) 

Po 
Here [¢,] denotes the velocity jump across a shock (downstream value 
minus upstream value) and the term in [¢,]* accounts for the change 
in stagnation pressure across any shock waves upstream of the point in 
question*. ‘To our order of accuracy [¢y] is a function of Y and Z only. 
If these expressions are now inserted in the drag integral (III.2b), there 
results 

D_ _ ([ noo ((V)? + M242 —2 (2—y)M463 —Uy + IME [be P 
15, U2 a } v ( (VG)? + MF 6% — 3 (2—-y) M5 4 — 3 (vy + 1) ME D [Px]*) - 
sre ~ @ “/8 
—2N.V¢(¢dy— M? 4%)} dS. (111.4) 


* The need for this term was brought to my attention by Professor Germain and 
Professor Cole. 
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This is the drag integral of second-order compressible flow theory. For 
transonic flow we may make the further approximations 
3 243 
M3 by as M; x 
(N.V¢)6% ~ Niydh, 


to obtain 


bpy U5, = 1) Nay (V9)? + Mody + PO) — OP 2 [bx }?)-2N. Vb oy} dS. 


0 Jd 

(111.5) 
It is readily verified, by means of Gauss’s theorem, the differential equation 
(III.1), and the conservation equations tor a shock wave, that this drag 
integral is invariant under changes of the surface S (provided that S+A 
remains a closed surface and that there are no singularities between different 





surfaces S). 
For plane flows with M, = 1, (III.5) reduces to 


= 1 D . mee . 

= Cp = 3 oy = | uv dx + (v? + 8 — 3D [u]*) ay}, (111.6) 

fe) (a) spo l eB JC 
where the integral is taken counter-clockwise in the physical plane, and 
where the reduced variables of (2.2) have been introduced. 

Next, following Germain (1957), we consider the form of the drag 
integral when the flow may be represented on C by a similarity solution 
ot the form considered in (1.2); that is, when 

x = u-"h(a), (n= 3m), 


so that A is some solution of the hypergeometric equation 


a(1—a)h''+[$—(3 + 5n)a]h’—1n(1+n)h = 0. (111.7) 
Equation (III.6) now becomes 
a C, = | - (u-" +120 ( (8m + 2a? — (20 + 3)at*)h’ + $(n? + n)a*h] du+ 


tu t32[( — 8432 + 2al2)h” — (40 + 2 )al? — a1? )h’] da} — 


“ | Y [u]* dy, (II1.8) 


where C’ is the image of C in the hodograph plane. Now the invariance 
property of the drag integral ensures that the integrand in (III.8) is an 
exact differential, and this may also be verified by means of (III.7). Hence 
it follows from the du-term that 


n+3/2 
ae = | +2) (2n-+ 3)al2\H + $(n-+ mah) oa 


3 
Cc’ 


~!/ ys fujeay. (19) 


6Jo 


In general C and C’ cross shock waves, and the integrated term in (III.9) 
is then only piecewise continuous: its discontinuities must be taken into 
account when (III.9) is evaluated. 
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Germain (1957) has pointed out that if C is a large contour in the physical 
plane with termini on y = 0, 7 = 0, the appropriate solution is that with 
n= 1, so that the integrated term in (III.9) vanishes in the limit. The 
drag is then seen to be entirely due to the entropy gain across shock waves. 


NOTE ADDED IN PROOF 


It is clear from (111.8) and (III.9) that an edge force is possible, in 
transonic flow with .V/, = 1, only if 2 = 3. Such a tlow has been studied 
by Nonweiler (1958), who finds that the solution 


¥ = const. w**(1—a)-** (w > 0) 

corresponds to the flow past the body Y o Y*°, the velocity on this body 
being sonic for Y > 0. Nonweiler calculates the edge force by integrating 
surface pressure on the related bodies Y « X’ and lettingv —~- ?. Thesame 
result is readily obtained from (III.8), in which the du-term and the shock 
term vanish. ‘This shows that the only assumption implicit in Nonweiler’s 
result is that the velocity in the field away from the singularity is given correctly 


to the lowest order. 
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A note on the compression of air through repeated 
shock waves 


By J. L. STOLLERY and D. J. MAULL 


Department of Aeronautics, Imperial College, University of London 
(Received 10 May 1958) 


SUMMARY 

The results of calculations of the compression of air by 
repeated shock waves are compared with the perfeci-gas values 
given by Evans & Evans (1956). ‘The comparison emphasizes 
the increasing divergence of real from perfect-gas results as 
shock strengths are raised. ‘lhe equations relating conditions 
across a shock wave are obtained in a convenient form for 
solution using a Mollier diagram. 
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Figure 1. The displacement vs time diagram. 








Figure 1 shows the problem considered. A piston instantaneously 
accelerated to a velocity u,, which is subsequently held constant, generates 
a shock wave. ‘lhe wave is reflected from the closed end and the piston 
face. ‘The air pressure, density and temperature rise non-isentropically 
through the initial and reflected shocks, soon reaching values at which the 
excitation of vibrational modes, dissociation, electronic excitation and 
finally ionization become important. ‘These real-gas effects are taken into 
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account here by using the Mollier diagram of Feldman (1957) prepared 
from the real-air tables of Hilsenrath & Beckett (1956). All viscous, heat 
transfer and relaxation effects are neglected and the flow is assumed to be 
one-dimensional. 

The notation is as follows. 

n An integer, referring to a flow region in the time—displacement 
diagram (figure 1). Thus H,,, p,, p, and T,, are the enthalpy, 
density, pressure and temperature in the mth region. The nth 
shock wave moves into the (m—1)th region and has the mth 
region behind it. 

U,, Velocity of the mth shock wave. 

u, Flow velocity in the nth region (u, =u, for n odd, u, = 0 
for n even). 

q Velocity relative to the shock wave. 

M, Shock strength = (initial shock velocity)/(speed of sound in 
region 4) = U,/ap. 

z A pressure of one atmosphere. 


The equations for a stationary shock wave are: 


energy H, + 3q2 = H, + 442, (1) 
momentum Po + po F5 = pi +P, qi; (2) 
continuity Pode = fr M- (3) 


These may be used for a moving shock wave by superposition of velocities, 
giving for the initial shock wave 

qo = Ui, q = Uy. 
Since u, = u,,, the above three equations yield 


H, ~ Hy = 14;(2=P, Pi-Po= uj( 2f0.) U,=;—2—. (4) 
P1— Po Pi— Po on (Po Px) 


In the same way the following equations for the first reflected shock 
are obtained: 
Ri a 1u;(B 22s), tin ¥5( p2Pr ), ——. a 
P2— Py P2— Py (p2/p,)— 1 
Similar formulae may be obtained for further shock reflections. ‘The 
general forms of the equations are: 





n+1—Pn 
of Pn+1P 8 ™ 
Pnii—Pn = uj( Pett), (7) 
Pn+1— Pn 
U,, = a (+ for n odd, — for meven). (8) 
= (p,, 1 P,) 


On the Mollier diagram, figure 2, enthalpy is plotted against entropy 
for lines of constant pressure, lines of constant density and lines of constant 


2T2 








652 J. L. Stollery and D. F. Maull 


temperature. ‘hus any two of these properties defines a point on the 
diagram. Given the initial conditions in region 0 and a value for M, 
(or piston speed), two values for p,; are chosen. ‘The two corresponding 
values of H, calculated from (4) are plotted on the Mollier diagram, and 
the points joined by a straight line. ‘The two values for the pressure p, 
from (+) are also plotted, and the line linking these points will cut that 
joining the enthalpy values if the chosen densities are reasonable. ‘The 
point of intersection gives the approximate values for conditions in 
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Figure 2. Mollier diagram showing method of solution. 


region 1. ‘The approximation may be improved by choosing two new 
density values just either side of the first intersection point and repeating 
the construction. With a little experience this becomes unnecessary. 
Now, knowing conditions in region 1, those in region 2 follow from (5) on 


using the same graphical construction. Similarly the method can be 


used for the other regions. The process is then repeated for further values 
of M,. A comparison with Feldman’s electronically-computed values 


The values for region 3 


2 showed a maximum variation of 3°,. 


for region 2 s 
are thought to be accurate to 5°,, 


The results for an initial temperature 7, of 290° K are shown in 
Conditions after one, two and three shocks are plotted 


figures 3, + and 5. 
against initial shock strength and compared with the perfect-gas solutions. 
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After two reflections only (m = 3), the perfect-gas solutions diverge from 
the air values for initial shock strengths greater than two. ‘Table 1 presents 
a comparison at various piston speeds. 

Thus pe-zfect-gas performance calculations for shock tubes of the 
multiple reflection type are invalidated and real-gas tables must be used. 
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Figure 3. Pressure vs shock strength. 7 and 0-017 refer to the values of po 
(7 = 1atmos.). The curve for py) = 0-17 lies between the other two curves 
for air and has been omitted for the sake of clarity. 





Piston Perfect gas Air (pp = 0-017) 
speed 


(ft./sec) P3/Po ts (Le P3/Po Ps/Po T3 To P3/Po 











2 000 90 5-1) 17°5 90 4:8 20-0 
4000 | 595 | 15-3|37:0|] 600] 10-6 | 57-0 
6 000 | 1500 | 31:0} 46-0} 1900 | 16-8 | 107-0 
8000 | 2700 | 53-4| 51-0] 5300 | 23-0 | 158-0 
10000 | 4300 | 83-0) 52-5|11300 | 28-2 | 222-0 


12000 | 6100 | 118-0 | 53-0 | 18900 | 33-0 | 292-0 


| | | 
Table 1. 
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Figure 4. Density vs shock strength. 



























































300777 ra a a | - fi =Sr 

ies, tae ae | | 

——+— PERFECT GAS /o 70 | 
250 5 Sane Y- iii lacus mcrae mains 
F A | | | 
/ at, | | 
200}——_—__+——}—— ++ —_444 7 | — 7 — \ ——_—_—— 4 
vv a | | | | 

. [07 | | | | 
i5O + T —F T 1 
Pt. | | 

“74 rf | 

a } | | | } 
100 ——— 4 —___  —___4—__— — SS = eee: | 
“A | o.o'8 | | | 

Z n=3 — i 

ZL Pe. eae Oe i 
” _ ZY ; nae = GE ae “ae jen i i paaeyl 
ZA bap | 
a — 
et —{——— {Wormleom— | — —t— na | 

re) 
r) 2 4 6 8 10 12 14 16 18 20 


Figure 5. Temperature vs shock strength. 
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REVIEWS 


History of Hydraulics, by HuNrer Rouse and Simon INce. Iowa 
City: lowa Institute of Hydraulic Research, State University of 
Iowa, 1957. 269 pp. $5.00. 


In these days of high blood pressure and rapid advancement in science 
and technology, it is refreshing to find a book which halts the headlong 
rush and contemplates the past. Here is a harmless tranquiliser to which 
all those interested in fluid mechanics should become addicted. ‘The 
historical development of hydraulics since early antiquity is briefly but 
expertly recounted in a way which enables the reader to see the progress 
of the subject in its proper perspective. The aim has been to review the 
advances made in formulating the basic principles of fluid motion, together 
with those relevant aspects of experimental research, engineering design, 
related scientific progress and world history. Hydraulicists will be aware 
that much of the material for the book has already been published as a 
series of supplements to La Houille Blanche which has regularly featured 
historical and biographical notes. 

The distinguished senior author, who is prominent in the advancement 
of engineering fluid mechanics, has also, most fortunately, been interested 
in the historical and philosophical progress of his subject. In 1950 he 
therefore encouraged a research student at Iowa, Simon Ince, to prepare 
a thesis on ‘“‘ A history of hydraulics at the end of the 18th century ” 
which was submitted in 1952. When, some years later as a Fulbright 
Scholar, Professor Rouse spent a year at the University of Grenoble, 
it was an excellent opportunity to expand the original thesis. He was 
greatly aided by the unrivalled library facilities which had been built up 
in the neighbouring Laboratoire Dauphinoise d’Hydraulique under the 
genius of another eminent hydraulicist, Pierre Danel. With such a fusion 
of backgrounds it is little wonder that the book flourished. 

Inevitably, because water is such an important element in life, 
‘Hydraulics’ must be one of the oldest subjects studied by man. 
Initially, its study was presumably empirical, but later it proved a 
fruitful field for theoreticians. One of the fascinating aspects of the 
book is its revelation of the close link in the past between the mathematician, 
the scientist and the hydraulic engineer, who in the 18th century 
were usually one and the same person. Nowadays in this era of 
specialization the field has broadened into what is fashionably termed 
fluid mechanics and the pursuit divided between aero- and hydro- 
dynamicists, engineers, meteorologists and a host of others who sometimes 
delude themselves that their own subject is ‘ different’. It is, in fact, 
often difficult to see the wood for the trees, and for this reason the book 
is welcome because of the clear view that it presents. 

The story begins in the days, difficult for historians, when ancient 
civilizations were developing the middle east. Little is recorded, but it 
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is interesting, for instance, to note that the siphon dates at least from those 
times. ‘The later periods of Greek and Roman culture are fortunately 
better documented in classical literature, from which numerous quotations 
are taken. ‘The theoretical contributions to hydrostatics by Archimedes 
and the philosophy of Aristotle and others are famous but it may not 
always be appreciated that the Greeks were also competent hydraulic 
engineers, their skills as canal and tunnel builders probably being inherited 
from the middle east. On the mechanical side, Ctesibus is accredited 
with the first recorded positive displacement pump, and the steam-jet 
turbine described by Hero is well known. ‘The Romans were outstanding 
as water handlers and have left behind many hydraulic structures as 
evidence of their engineering skill. Although illustrations from the 
middle ages depict the use of hydraulic machines in mills and forges, no 
outstanding figures in hydromechanics seem to have emerged. Examples 
are quoted of a few gropings towards mechanics made by Buridan of the 
University of Paris, Heytesbury of Oxford and others, but little is known 
of scientific progress during this period. In addition to the references 
quoted by Rouse and Ince the translated modern edition of De Re Metallica 
by Georgius Agricola might also have been mentioned as a valuable source 
of illustrations of hydraulic machinery used in mining during the middle 
ages. However, not until Leonardo da Vinci and the Renaissance do the 
individual contributions of hydraulicists become apparent. The drawings 
of hydraulic phenomena made by da Vinci and reproduced in the book 
show da Vinci’s acute powers of observation, and extracts from his writings 
reveal a remarkable understanding of physical laws such as those of force 
and work. Others such as Copernicus, Stevin and Galileo also contributed 
to the progress of mechanics and the experimental method during this 
early period. From here on the story moves more rapidly through chapters 
on the naissance of the experimental method, post-renaissance hydraulics, 
17th-century mathematics and mechanics, to the 18th century when 
theoretical hydrodynamics really began and remarkable advances were made. 
The 16th, 17th and 18th centuries may be the most fascinating periods for 
many readers, partly because they cover the birth and development of 
hydrodynamics as we know it and partly because we are less familiar 
with these times than with the more recent developments in the 19th 
and 20th centuries. ‘The versatility of many of the scientists is remarkable. 
For instance, Newton is said to have discovered the contraction of jets in 
comparing the theoretical and measured flow through orifices and Euler 
was the first to consider cavitation in connection with the design of water 
turbines. ‘The progress made in the 18th century is truly remarkable 
and it makes a fascinating story. ‘The history is completed by chapters 
on the 19th century, the rise of fluid mechanics in the 20th century and 
an appraisal of the present situation at mid-century. Somehow the 
advances made in the 19th and 20th centuries lack the romance and glamour 
of the earlier golden age but the practical contributions, particularly to 
engineering, are nevertheless very great. 
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In the 18th century the distinguished tradition of the French Ecoles 
was clearly established by the remarkable numbers of very eminent 
hydraulicists which they produced. Because of the easier availability of 
French records for this period, there may be a tendency to neglect the 
work of other countries during this time. Balance is always difficult to 
preserve, particularly in such a concise history, and although one may 
have doubts here and there, the authors have maintained a good general 
level of equilibrium. 

Most of the important figures seem to have been mentioned, although 
a few distinguished names may have been omitted because they did not 
come within the particular scope of the book. Perhaps some mention 
should have been made of such names as Stodola, Schoklitsch and a few 
others. 

This book provides a good many details which present-day scholars 
miss by not having easy access to original or translated versions of the 
great technical works. How much is lost by reference to second-hand 
and often incorrect versions resulting from successive generations of 
textbooks! ‘‘ One of the ironies of history is the fact that Froude’s 
name has since become inseparably associated with a law of similarity 
and a non-dimensional number, the first of which he did not originate 
and the second of which he never even used.” ‘The authors have rendered 
a service by pointing out such anomalies and by quoting direct from old 


” 


works whenever possible. 

The standard of printing is high, and there are many illustrations, 
including reproductions of portraits of the outstanding figures. ‘The book 
is written in a concise and easily readable style and is remarkably complete 
and well-documented for a brief survey. In fact, it whets the appetite 
for a larger version and it is to be hoped that the authors may consider 
writing a more detailed history on the same lines. ‘To all those wandering 
in the maze of modern technology, this history of hydraulics is warmly 
commended. 

S. P. Hutton 


Fundamentals of Gas Dynamics, edited by H. W. Emmons (Volume 
III of High Speed Aerodynamics and Jet Propulsion). Princeton 
University Press, 1958. 749 pp. $20.00. 


Here is a book written by some of the world’s experts on gas dynamics. 
The result, as one would expect, is a book of great value which will be used 
continually by everyone interested in the subject. It must also be said, 
however, that this reviewer found it disappointing in a number of ways. 
This should not be unduly emphasized, and if a relatively large amount of 
space is devoted here to points of adverse criticism, it is with the under- 
standing that the reader does not need or want to be told many times that 
articles by people like ‘I’. von Karman and G. I. Taylor have much to offer! 
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The eight main sections are: 

A. The equations of gas dynamics (Tsien). 

B. One-dimensional treatment of steady gas dynamics (Crocco). 
One-dimensional treatment of unsteady gas dynamics (Kantrowitz). 


C 
D. The basic theory of gas dynamic discontinuities (Hayes). 
E. Shock wave interactions (Polachek and Seeger). 
F. Condensation phenomena in high speed flows (Stever). 
G. Gas dynamics of combustion and detonation (von Karman, Emmons, 
Tankin, Taylor). 

H. Flow of rarefied gases (Schaaf and Chambre). 

One of the main criticisms is that these sections remain essentially 
separate, each one having a different purpose and scope, and a different 
assumption as to the type of reader. ‘The material is not combined into a 
connected comprehensive account. ‘There are practically no cross- 
references between sections, and there is considerable duplication. All 
this could have been just a minor irritation but it turns cut to be more. 
For as a consequence, the value of some articles is reduced, matter of 
interest is left out and some basic questions emerge confused. Let us 
consider the articles in turn. 

A first section of the standard type on the equations of motion is 
appropriate if the book is planned as a complete account that a newcomer 
would read. But the later sections (with the possible exception of C) are 
not of that type. Thus A loses much of its value. Most of it is in fact 
not used at all in the rest of the book; each section rederives the few 
equations needed from first principles. ‘There would be a good case for 
omitting the section altogether (since there is no shortage of references) 
and replacing it by a careful review of the results of kinetic theory, the 
physical chemistry of dissociation, etc., topics which have become so 
important and yet are not familiar to many gas dynamicists. (It is true 
that Vol. I of the series is Thermodynamics and Properties of Matter but 
that is another 700 pages!). 

Although a detailed discussion of the contents of each section would be 
too long, it must be noted that some of the derivations in A are misleading 
and some appear to be wrong. For example, in deriving the energy 
equation (pages 10 and 11), the work done by the stresses 7,; per unit time 
per unit volume is calculated to be z,; du, 0x; instead of O(7;;u;),Ax,;; the 
resulting equation is correct because the increase of kinetic energy is also 
omitted from this energy equation. In the treatment of Helmholtz’s 
equation for the vorticity Q (pages 30 and 31), the term (Q. V)u is 
interpreted as a ‘ bending’ of vortex lines and the term — 92(V.u) as 
the ‘stretching’. This is wrong, of course; the former represents 
the stretching as well and the latter would be absent altogether in 
incompressible flow. The so-called ‘intrinsic coordinates’ m, s, the 
distances along and perpendicular to streamlines, are used on pages 
53, 61 and 63. ‘These variables cannot be introduced in the large and 
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0/On, 0/ds are not derivatives in the usual sense; they stand for h>! d/dq, 
hz 0/08, where « and £ are generalized coordinates defining the streamlines 
and orthogonals, and h, dx and h, df are the line elements; for instance 
0/ds(0f on) # O On(Of ds). A final example: the equations of motion are 
deduced by considering conservation of mass, etc., for a parallelepiped, 
and are then integrated over an arbitrary volume, the results being interpreted 
as conservation laws. Why not treat an arbitrary volume in the first place ? 

Section B is a complete book in itself. It is nearly 300 pages long, 
out of about 750 for the whole book, and is completely independent of the 
rest. All aspects of one-dimensional steady flow are developed and 
applied in considerable detail. ‘This must surely be the most thorough 
treatment that exists. Particularly impressive is the amount of useful 
information that can be got out of rather a small number of basic ideas, 
supplemented by careful physical arguments. There are many practical 
applications, some of which are far from one-dimensional at first sight. 
Like most of the other articles, great care is taken to be realistic about the 
thermodynamics; after reading this book, it will be almost embarassing 
to assume p o ev, ‘This section is perhaps more suitable for reference 
as needed than for continuous reading; reading it from beginning to end 
required some persistence. 

In contrast to B, section C is relatively short and not compendious, 
with a general reference to Supersonic Flow and Shock Waves by Courant 
& Friedrichs for full mathematical details. Here some of the space wasted 
elsewhere in duplication could have been used. For instance, Riemann’s 
classical solution for the interaction of two simple waves and the hodograph 
transformation are mentioned only in a footnote. On the other hand, 
there are some physical arguments about simple waves, with or without 
weak shock waves, and Professor Kantrowitz’s original work on the 
stability of shock waves in steady channel flows, which do not appear in 
Courant & Friedrichs. It should be noted that this particular section was 
written in 1952 and as a consequence is a little out of date now. 

Part of D is concerned with the Navier-Stokes structure of shock 
waves. ‘This is very well done and indeed some of the original ideas, 
such as the precise use of Burgers’s equation, have less impact than they 
deserve because they were introduced independently, but published 
earlier, by Lighthill in Surveys in Mechanics (the G. I. ‘Taylor 70th 
Anniversary Volume). It is disappointing, however, to find that the 
kinetic theory approaches to the question are not also included in some 
detail. ‘The author claims correctly that their foundations are a little 
shaky and their results seem to be less successful than the Navier-Stokes 
treatment, but these are not sufficient reasons unless space is scarce. 
But space is used elsewhere on much less respectable theories, and section 
B is certainly not limited in this respect. After all, the beginner can find 
the classical treatment of shock structure in many books, and the more 
experienced reader would welcome an authoritative account of some of 
the more recent attempts just because they do pose considerable difficulties. 
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The remainder of D concerns ‘ discontinuities ’ in which reactions and 
change of state occur, examples being condensation shocks, flame fronts 
and detonations. ‘This has to be considered in conjunction with sections 
F and G where these topics are studied in detail. Now, F and G make no 
serious reference to the longish discussion in D; instead there is 
appreciable overlap. For example, the classification of combustion pheno- 
mena into weak and strong deflagrations and detonations is given in D3, 
G1 and G10, in each case as if for the first time. Both D3 and G10 include 
detailed discussions of the properties of the Hugoniot curve which are 
essentially the same. ‘The equations for a plane stationary flame front 
are given in both G1 and G2 (which are by different authors) in different 
form and notation without explaining why this was necessary. In a 
difficult and only partially understood subject, the lack of a unified treat- 


4 


ment is very confusing. Moreover, what is the relation of these articles 
on combustion to Vol. II of the series which is an equally large volume 
devoted entirely to combustion ? 

‘The interesting question of the existence of strong deflagrations and 
weak detonations is clouded by some of the statements. After referring 
to Friedrichs extensively in the arguments suggesting that strong 
deflagrations cannot exist, the final paragraph of D6 begins: ‘* With 
regard to the weak detonation, which has not been observed (again to the 
writer's knowledge) in explosion experiments, there is no theoretical 
evidence for non-existence.”” Yet such evidence is given in Courant & 
Friedrichs and on p. 231 of their book there is the italicized statement: 
‘Consequently, weak detonations are impossible.” I am open to argu- 
ments that their reasoning is unsound but not that it is not there. 

Of course, in spite of the above criticism, F and G are bulging with 
information, and one can really get some idea of the present state of 
these complicated fields. There is again the feeling of reality which 
is an outstanding feature of the book. For example, in G10 we are 
led through the details of calculating the Hugoniot curve for the mixture 
H,+ $0,+2-5N,, first without and then with dissociation, before going 
on to discuss the propagation of detonation waves. At every stage there 
are checks with whatever experiments are available. 

Section E requires little comment. It is a complete account of the 
present knowledge of the interaction of plane shock waves with other 
waves, solid walls and interfaces. ‘The methods of deriving theoretical 
formulae, the subsequent computations, the comparison with experiment 
(and lack of agreement for Mach reflection of weak shocks) are all set out 


exactly as they should be. 

The final section H stands out as a piece of careful exposition. All the 
ideas and results are presented clearly and concisely; it makes interesting 
and informative reading. It even remedies earlier deficiencies in part. 
For example, the section on slip flow contains an excellent review of the 
derivation of the Navier-Stokes, Burnett and Thirteen-moment equations 
from the Boltzmann equation, with a sober consideration of the validity 
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and degree of success of the various approximations involved. ‘This is of 
general value for the book, quite apart from its relevance to slip flow. Some 
brief comments on shock structure also appear here. Odd though it may 
seem, section H was added on to this volume as an afterthought; it was 
originally prepared for Vol. IV. 

After a classification of the various flow regimes, with the significant 
relative magnitudes of Mach number and Reynolds number, section H is 
divided into three main sections; free molecule flow, slip flow, and experi- 
mental results in slip flow and transition regimes. ‘The first of these shows 
how theories of the heat transfer and the forces, for bodies moving at the 
appropriate high speeds, may be developed assuming a Maxwellian 
distribution about the main stream velocity for the molecules incident on 
the body surface and using empirical coefficients to specify the interaction 
of the incident molecules with the surface. ‘Typical cases are worked out 
and detailed references are given for further work. In the chapter on 
slip flow, the review noted above is given and the conclusion reached that 
the Navier-Stokes equations, together with the slip velocity and tem- 
perature jump boundary conditions, are at present more successful than 
any of the intended improvements. ‘lhe final chapter sets out the 
comparison of theoretical results worked out in this way with experiment. 

G. B. WHITHAM 
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